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Preface 


The  Air  Force  has  seen  a  growing  number  of  weapon 
systems  that  employ  some  type  of  digital  signal  processing 
(DSP).  Many  analog  circuits  are  being  replaced  with 
digital  equivalent  circuits.  AFIT  supports  this  trend 
with  courses,  laboratory  projects,  and  theses  designed  to 
give  students  sufficient  background  in  this  area  to  satisfy 
Air  Force  needs.  The  computer  facilities  in  the  AFIT 
engineering  laboratory  are  an  important  part  of  this 
support . 

One  of  the  computer  systems  has  recently  had  more 
hardware  added  to  it,  making  it  potentially  a  very  valuable 
resource  in  support  of  DSP.  This  report  describes  what 
was  accomplished  with  the  system. 

Thanks  to  Professor  Gary  Lamont,  advisor,  who  proposed 
the  topic  and  helped  limit  its  scope.  His  interest,  en¬ 
couragement,  and  suggestions  were  invaluable.  Thanks  to 
John  Bankovskis,  currently  assigned  to  the  Optical  Sensors 
Lab,  who  volunteered  his  time  and  expertise  to  assist  in 
the  generation  of  the  operating  system.  Thanks  also  to 
Maj  Lillie  and  Capt  Kizer,  my  thesis-committee  members. 


Wayne  Todd 


t 


Table  of  Contents 


Preface  .  ii 

List  of  Figures  . vi 

Abstract  .  vii 

Chapter  1  Introduction  .  1 

1 . 1  Background  .  1 

1.2  Requirements  .  1 

1 . 3  Scope  .  2 

1.4  Approach  .  2 

1.5  Overview  of  Thesis  .  3 

Chapter  2  Requirements  .  5 

2.1  Operating  System  .  5 

2.2  Data  Acquisition  .  7 

2.3  Data  Storage  .  8 

2.4  Computational  Tasks  .  9 

2.5  Output  Display  . 10 

2.6  User  Interaction  . 11 

2.6.1  The  User  . 12 

2.6.2  Response  Time  . 13 

2.6.3  Language  . 14 

2.6.4  Input  . 15 

2.6.5  User  Assistance  . 16 


Chapter  3  Environment  .  17 

3.1  Operating  System  . 17 

3.1.1  Generation  of  RTE-III  . 19 

3.1.2  Capabilities  and  Limitations  . 21 

3.2  Data  Acquisition  . 22 

3.2.1  A/D  Converter . 22 

3.2.2  Software  Support  for  A/D  Converter  . 23 

3 . 2 . 2 . 1  Pol 1 ing  . 23 

3. 2. 2. 2  Interrupts . 27 

3.2.2. 3  Direct  Memory  Access  . 28 

3. 2. 2. 4  I/O  Driver . 29 

3.3  Data  Storage . 30 

3.3.1  Disk  . 30 

3.3.2  Cassette  . 31 

3.4  Output  Display  . 32 

3.4.1  Conventional  CRT  Terminal  and  Line  Printer  ....32 

3.4.2  Graphics  Terminal/Printer  . 33 

3.4.3  D/A  Converter  . 33 

3.4.4  Software  Support  for  D/A  Converter . 34 


i  ii 


Chapter  4  Computational  Tasks  .  35 

4.1  FFT/IFFT  .  36 

4.2  FIR  Filter  Design  .  37 

4.3  Autocorrelation  and  Covariance  .  38 

4.4  HR  Filter  Design .  38 

4.5  Real-Time  Processing  .  39 

Chapter  5  User  Interaction  .  42 

5.1  Graphics  Terminal  Control  .  42 

5.2  Segmented  Program  .  43 

5.3  Executive  Programs  .  44 

5.4  Recommended  Software  Structure  .  45 

Chapter  6  Recommendations  .  46 

6.1  Additional  Hardware  .  46 

6.1.1  Magnetic  Tape  Drive .  46 

6.1.2  Faster  A/D  Converter .  47 

6.1.3  HP  Fourier  Analyzer  .  48 

6.1.4  Graphics  Printer  .  48 

6.2  DMA  for  D/A  Converter  .  48 

6.3  Additional  Applications  Programs  .  49 

6.4  Real-Time  Processing .  49 

6.5  Operating  System  .  49 

Bibliography  .  51 

Appendix  A:  RTE-III  Generation  Answer  File  .  54 

Appendix  B:  Using  SWTCH  and  DSKUP  .  59 

Appendix  C:  HP  Users  at  WPAFB  .  61 

Appendix  D:  Computational  Tasks  Listings  .  62 

1.  FAFT4  Source  Listings  .  63 

2.  FAFT4  Original  Output  .  68 

3.  FAFT4  Output  in  Auto  Plot  Format  .  70 

4.  FIRF2  Source  Listing  .  71 

5.  FIRF2  Output  .  85 

6.  COVAU  Source  Listing  .  87 

7.  COVAU  Output  .  92 

8.  FWIIR  Source  Listing  .  93 

9.  FWIIR  Input  . 106 

10.  FWIIR  Output  . 107 

Appendix  E:  A/D  Converter  Software  Listings . 119 

1 .  Pol  1  ing  . 120 

2.  DMA  . 127 

3.  Dummy  Drivers  . 130 

Appendix  F:  System  User's  Guide  .  132 


Appendix  G:  Disk  Storage  Programs . 135 

1 .  BINRY  Utility  . 136 

2.  FMP  Calls  . 139 

Appendix  H:  Auto  Plot  Program  . 141 

Appendix  I:  D/A  Converter  Software  Listings  .  144 

Appendix  J:  Segmented  Program . 152 

Appendix  K:  Executive  Programs . 155 

1.  Program  Schedule  EXEC  Call . 156 

2.  MESSS  Utility . 159 

Appendix  L:  Recommended  Software  Structure  .  161 

Vi  ta . 165 


v 


List  of  Figures 


Block  Diagram  of  Configured  Hardware . 

Connections  for  A/D  Converter  and  Interface 

Oscilloscope  Display  of  Status  Bit . 

Plot  of  Sample  Time  Measurements . 

Flowchart  of  Real-Time  Processing . 


Abstract 


This  project  involved  the  integration  of  hardware 
and  software  to  develop  a  general  purpose  digital  signal 
processing  (DSP)  laboratory  based  on  a  Hewlett  Packard  (HP) 
21MX-type  computer.  The  hardware  included  an  HP  2108A 
computer,  HP  7906  disk  drive,  HP  26/8A  graphics  terminal, 
Texas  Instruments  Silent  700  printing  terminal,  Remex 
optical  paper  tape  reader,  two  Analog  Devices  ADC-12QZ 
analog- to-digital  converters,  and  an  HP  12555B  digital-to- 
analog  converter. 

The  operating  system  generated  for  the  laboratory 
was  RTE-III.  It  is  the  latest  version  of  RTE  that  can 
be  used  with  the  currently  available  hardware.  Software 
support  for  I/O  devices  not  supported  by  RTE-III  was 
developed  and  implemented. 

Four  DSP  applications  programs  were  implemented  and 
tested.-- The  first  program  performs  both  a  fast  Fourier 

N 

transform  ( FFT )  and  an  inverse  fast  Fourier  transform 
(IFFT)  of  a  complex  function.  The  FFT  and  IFFT  are 
computationally  efficient  methods  for  time-to-f requency 
and  f requency-to- time  transformations,  respectively.  The 
second  program  designs  a  linear  phase  finite  impulse 
response  digital  filter.  The  third  program  uses  auto¬ 
correlation  and  covariance  methods  of  linear  prediction 
analysis.  The  fourth  program  designs  a  finite  word-length 
infinite  impluse  response  digital  filter. 


Chapter  1  Introduction 


This  chapter  begins  with  background  information  on 
the  motivation  supporting  this  study.  This  is  followed 
by  a  brief  look  at  the  requirements  for  a  general  purpose 
digital  signal  processing  (DSP)  laboratory,  the  extent  to 
which  these  requirements  have  been  met,  and  the  approach 
taken.  The  chapter  concludes  with  an  introductory 
description  of  the  thesis  project. 

1 . 1  Background 

The  AFIT  engineering  laboratory  has  had  an  HP  21MX 
computer  system  for  several  years.  It  has  not  been  used 
extensively  by  either  students  or  faculty  because  its 
limited  peripheral  devices  and  limited  software  support 
made  its  use  difficult.  Recently,  additional  hardware 
(e.g.,  disk  drive,  graphics  terminal,  memory)  has  been 
added  to  the  system  which  would  potentially  make  it  a 
very  valuable  asset  to  AFIT  engineering  students  and 
faculty.  These  additional  capabilities  have  not  previous¬ 
ly  been  exploited.  Also,  professor  Gary  Lamont,  who  is 
responsible  for  the  system's  use,  has  designated  it  to 
function  as  a  DSP  laboratory.  DSP  is  an  important  part 
of  the  engineering  curriculum,  and  additional  support 
facilities  are  needed  to  meet  the  demands  of  a  growing 
number  of  students. 

1.2  Requirements 

The  basic  requirement  of  the  laboratory  is  that  it 
be  user-oriented  and  provide  hardware  and  software  to 


support  most  common  DSP  applications.  It  should  also 
be  expandable  to  allow  the  integration  of  additional  hard¬ 
ware  and  software.  The  basic  requirement  has  been  divided 
into  six  requirements  for  the  purpose  of  discussion  and 
implementation.  These  are:  1)  the  operating  system, 

2)  data  acquisition,  3)  data  storage,  4)  computational 
tasks,  5)  output  display,  and  6)  user  interaction. 

1 . 3  Scope 

This  project  investigated  and  demonstrated  most  of 
the  system's  capabilities  as  they  existed  at  the  beginning, 
and  as  they  were  added  during  the  project.  All  the  hard¬ 
ware  had  been  specified  and  ordered  previously.  Software 
has  been  provided  (including  an  operating  system)  where 
necessary  in  support  of  the  hardware.  Also,  several 
DSP  applications  programs  were  taken  from  the  literature 
and  modified  to  run  on  the  HP  system.  Investigation  of 
the  present  capabilities  led  to  conclusions  about  what 
additions  should  be  made  to  meet  the  ultimate  design 
goal.  The  recommended  additions  are  presented  as  either 
essential  or  desireable. 

1.4  Approach 

The  first  task  was  to  investigate  the  requirements 
for  a  general  purpose  DSP  laboratory.  This  involved 
both  a  literature  search  and  personal  conversations  with 
several  people  currently  using  DSP.  The  next  task  was 
to  generate  an  operating  system  (OS).  The  OS  provides 


software  support  for  most  of  the  I/O  devices  used.  It 
also  has  program  development  and  testing  facilities  that 
were  necessary  later  in  the  project.  The  OS  was  generated 
on  another  HP  computer  system  because  a  disk  loader  ROM, 
which  is  necessary  to  load  the  OS  into  main  memory,  was 
rot  available  early  in  the  project. 

Next,  the  hardware  was  configured.  Most  of  it  was 
available  from  the  start,  but  some  (e.g.,  the  disk 
loader  ROM)  arrived  later.  RTE-III  supports  most  of 
the  hardware,  but  software  support  had  to  be  provided  for 
the  A/D  converters  and  the  D/A  converter.  Four  methods 
of  software  support  for  the  converters  were  investigated: 
polling,  interrupts,  direct  memory  access,  and  system 
I/O  drivers.  This  development  followed  the  hardware 
integration . 

Finally,  several  DSP  applications  programs  were 
implemented  and  tested. 

1.5  Overview  of  Thesis 

Chapter  2  covers  the  requirements  that  have  been 
identified  for  the  laboratory.  The  requirements  are 
a  combination  of  observing  similar  systems  and  the 
author's  own  judgement.  Chapter  3  discusses  the  system's 
hardware  and  support  software,  which  together  comprise 
the  laboratory's  environment.  The  main  subject  areas 
covered  are  the  operating  system,  data  acquisition, 
data  storage,  and  the  output  display.  Chapter  4  is  a 
discussion  of  the  DSP  applications  programs  that  have 
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been  implemented  and  tested.  The  four  programs  presented 
were  taken  from  a  textbook  (Ref  7)  and  modified  to  run 
on  the  HP  system.  Also  presented,  is  a  possible  way  to 
use  the  system  for  real-time  processing.  Chapter  5  is  a 
discussion  of  the  issues  which  have  been  addressed  per¬ 
taining  to  user  interaction.  Chapter  6  lists  and  justifies 
numerous  recommendations  which  would  allow  the  DSP 
laboratory  to  meet  all  requirements. 
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Chapter  2  Requirements 

A  literature  search  did  not  lead  to  any  book,  per¬ 
iodical,  or  other  document  which  specified  and  justified 
either  the  hardware  or  software  requirements  for  a  general 
purpose  DSP  laboratory.  However,  existing  hardware  and 
software  implementations  were  studied  (Ref  4:433-441,  7, 
10,  12,  18,  23,  25:655-657,  31,  32,  33,  Appendix  C)  in 
an  attempt  to  determine  the  most  common  requirements, 
and  survey  possible  capabilities.  Six  basic  requirement 
areas  have  been  identified:  the  operating  system,  data 
acquisition,  data  storage,  computational  tasks,  output 
display,  and  user  interaction.  Even  though  each  of  these 
areas  is  discussed  in  a  separate  section  of  this  chapter, 
they  are  closely  interrelated  and  must  be  considered 
together  when  implementing  a  DSP  laboratory.  Note  that 
for  brevity,  in  this  context,  the  term  operating  system 
includes  the  assembler,  interactive  editor,  compilers, 
and  system  utilities. 

The  overall  objective  is  to  design  a  user-oriented 
system  so  that  the  DSP  system  user  does  not  require  a 
computer  background.  This  objective  is  to  be  fullfilled 
in  both  an  educational  and  research  environment. 

2.1  Operating  System 

Since  the  user  of  this  DSP  laboratory  is  assumed 
to  have  minimum  knowledge  of  the  specifics  of  computer 
operations,  automated  resource  management  is  essential. 

The  operating  system  (OS)  provides  this.  In  effect, 
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the  details  of  resource  management  should  be  "user- 
transparent"  (Ref  20).  These  resources  include  the  CPU, 
memory,  peripheral  devices,  and  files.  Because  the  OS 
manages  these  resources,  it  is  an  important  part  of  user 
interaction,  discussed  separately  in  section  2.6. 

An  OS  can  be  classified  as  a  real-time  OS,  batch  OS, 
time-sharing  OS,  or  multipurpose  OS  (Ref  19:61-99).  A 
real-time  OS  is  characterized  as  being  able  to  perform 
a  sequence  of  processing  tasks  (e.g.,  real-time  data 
acquisition,  computations,  and  display)  within  rigid  time 
constrain ts .A  batch  OS  is  capable  of  executing  more  than 
one  job  without  the  need  for  operator  intervention.  The 
concept  of  time-sharing  refers  to  distributed  computer 
power  among  a  set  of  users  at  interactive  terminals. 

A  multipurpose  OS  will  support  the  capabilities  of  a  com¬ 
bination  of  at  least  two  of  the  other  three  types  of  OSs. 

The  DSP  laboratory  requires  a  multipurpose  OS.  It 
is  intended  to  have  the  capability  for  real-time  pro¬ 
cessing  as  described  in  section  4.5.  Also,  many  DSP 
programs  have  lengthy  (one  or  more  hours)  execution  times. 
Therefore,  a  multiprogramming  batch  capability  is  also 
required,  so  that  the  user  can  initiate  actions  while 
one  or  more  programs  are  executing.  The  only  aspect  of 
time-sharing  required  is  the  use  of  an  interactive 
terminal.  This  is  because  the  system  is  being  designed 
such  that  there  can  only  be  one  user  at  a  time.  An 
interactive  terminal  is  more  user-responsive  than  a 
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batch  system  because  it  allows  on-line  inputs  to  be 
entered  by  the  user.  Finally,  as  an  additional  aid  in 
user-responsiveness,  the  user  should  be  able  to  initiate 
batch  programs  interactively. 

2.2  Data  Acquisition 

Data  can  be  acquired  through  a  program  that  generates 
data,  from  a  data  file,  or  from  an  external  source.  In 
the  case  of  an  external  source,  the  phenomena  is  measured 
by  some  device  and  is  represented  as  an  analog  signal. 

An  analog-to-digital  (A/D)  converter  will  convert  the 
analog  signal  (samples  taken  in  real-time  and  equally 
spaced  in  time)  into  a  form  (binary  numbers)  suitable 
for  digital  computer  processing. 

Several  problem  areas  exist  with  regard  to  repre¬ 
senting  analog  data  in  digital  form.  The  first  is  the 
accuracy  and  precision  of  the  A/D  converter  The  pre¬ 
cision  is  determined  by  the  number  of  bits  which  repre¬ 
sent  each  data  point.  The  accuracy  concerns  the  ability 
of  the  A/D  converter  to  provide  correct  data  within  the 
precision  limit.  Two  other  common  problem  areas  are 
aliasing  and  leakage  (Ref  2). 

Aliasing  refers  to  the  fact  that  high-frequency 
components  of  a  time  function  can  impersonate  low  fre¬ 
quencies  if  the  sampling  rate  is  too  low.  To  effectively 
deal  with  aliasing,  the  sampling  frequency  must  be  at 
least  twice  that  of  the  highest  frequency/  W  present  in 
the  signal  being  analyzed.  However,  increasing  the 


sampling  frequency  results  in  more  data  points  and  con¬ 
sequently,  more  computations.  For  band  limited  signals, 
Chen  (Ref  4:82-83)  recommends  a  sampling  frequency  between 
2.0W  and  2.5W. 

The  problem  of  leakage  arises  from  sampling  a  sig¬ 
nal  for  a  finite  period  and  neglecting  what  happened 
before  or  after  this  period.  This  is  equivalent  to  multi¬ 
plying  the  signal  by  a  rectangular  data  window,  which 
results  in  the  leakage  of  energy  from  one  discrete  fre¬ 
quency  into  adjacent  frequencies.  Various  types  of  data 
windows  can  be  used  to  reduce  leakage  (Ref  24). 

The  user  should  be  able  to  interactively  specify 
the  sampling  rate,  the  number  of  samples  per  data  frame, 
the  type  of  data  window  (if  any),  and  the  data  source 
(e.g.,  an  I/O  channel  number).  By  changing  the  first 
three  of  these  four  variables  and  analyzing  the  same 
data  several  times,  error  analysis  can  be  performed 
with  respect  to  aliasing  and  leakage.  However,  except 
for  periodic  calibration,  the  user  has  to  assume  that 
the  data  supplied  by  the  A/D  converter  is  correct. 

2.3  Data  Storage 

Data  acquired  from  an  external  source  or  generated 
by  a  program  is  initially  stored  in  main  memory.  In 
the  case  of  the  external  source,  it  is  because  main 
memory  is  faster  than  secondary  memory.  The  data  can 
then  be  sent  to  a  display  device,  be  used  in  computations, 
or  stored  on  a  mass  storage  device,  such  as  a  disk. 


The  capability  for  storing  data  is  especially  important 
because  some  data  is  relatively  expensive  or  difficult 
to  obtain.  For  example,  it  may  be  expensive  to  turn  on 
a  rocket  engine  for  the  purpose  of  obtaining  sensor 
data.  This  kind  of  test  may  also  be  difficult  to  set  up. 

Another  reason  for  storing  data  is  that  it  may  be 
desired  to  use  the  same  data  more  than  once.  For 
example,  it  can  be  recalled  and  processed  by  one  algor¬ 
ithm  or  several  algorithms,  and  results  compared. 

2.4  Computational  Tasks 

The  computational  tasks,  or  algorithms  presented 
here  are  for  a  general  purpose  DSP  laboratory.  However, 
specific  applications  may  require  additional  types  of 
computation,  or  analysis  of  several  methods  of  accomp¬ 
lishing  the  same  task.  Therefore,  the  overall  software 
structure  should  be  modular,  so  that  additional  modules 
can  be  added  without  the  need  to  alter  the  existing 
software.  A  module  is  a  separate  entity  of  instructions, 
which  by  itself  performs  a  predefined  task,  computational 
or  otherwise.  Computational  tasks  required  for  a 
general  purpose  DSP  laboratory  include  the  fast  Fourier 
transform  (FFT),  inverse  FFT  (IFFT),  transfer  function, 
coherence  function,  auto  correlation,  cross  correlation, 
convolution,  waveform  averaging,  various  statistical 
computations,  digital  filter  design,  different  types  of 
data  windows,  and  power  spectral  density  (Ref  32). 

The  system  should  have  the  capability  of  performing 
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FFT 


these  tasks  in  real-time,  if  appropriate  (e.g., 

I FFT ,  correlation,  etc.).  To  insure  this,  the  CPU, 
memory,  and  executing  programs  must  be  fast  enough  to 
allow  real-time  computations  without  the  loss  of  data 
(i.e.,  losing  data  by  replacing  old  data  with  new  data 
before  the  old  data  has  been  used).  In  this  case, 
computations  must  be  performed  in  parallel  with  real-time 
data  acquisition  and  real-time  display. 

2.5  Output  Display 

A  general  purpose  DSP  laboratory  requires  the  cap¬ 
ability  for  displaying  data  and  results  from  the  afore¬ 
mentioned  computational  tasks.  Also,  a  display  device 
(e.g.,  a  graphics  terminal)  is  an  essential  communica¬ 
tion  medium  between  the  user  and  the  computer  in  an 
interactive  environment.  However,  a  discussion  of  this 
second  aspect  of  the  output  display  is  reserved  for 
section  2.6. 

Requirements  for  the  output  display  must  be  specified 
both  in  terms  of  the  types  of  devices  and  the  format 
of  the  display.  Several  criteria  exist  for  selecting 
the  display  devices.  First,  the  system  should  provide 
the  user  with  a  hard  copy  of  any  useful  display,  for 
documentation  purposes.  This  might  include  a  line 
printer  for  listing  tables  of  values,  or  a  camera  for 
taking  a  picture  of  an  oscilloscope  display.  A  soft 
copy  device  (CRT)  is  necessary  for  some  types  of  real¬ 
time  display  where  the  display  is  continuously  replaced 
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with  a  new  display.  It  is  also  important  to  consider  whet¬ 
her  the  soft  copy  device  is  a  storage  or  non-storage  type. 

A  storage  device  does  not  require  continuous  input  from  the 
computer  for  a  display,  whereas  a  non-storage  device  does. 

Rather  than  stating  that  one  or  more  particular 
devices  be  required,  the  combined  capability  of  a  par¬ 
ticular  group  of  devices  should  be  considered  in  light  of 
the  relevant  issues  presented.  In  addition,  some  applica¬ 
tions  may  require  an  output  which  is  not  a  visual  display. 
For  example,  an  audio  amplifier  and  speaker  are  needed 
for  speech  synthesis.  A  general  purpose  DSP  laboratory 
will  not  be  able  to  support  all  applications,  but  it 
should  be  flexible  enough  to  allow  additional  output 
devices  to  be  added. 

The  information  displayed  should  be  in  a  format 
which  is  convenient  and  comprehensible  to  the  user. 
Excluding  the  type  of  display  device,  the  type  of  output 
format  should  be  governed  primarily  by  the  information 
the  user  hopes  to  extract  from  the  output.  For  example, 
a  graph  generally  provides  a  good  panoramic  view  of  a 
set  of  data,  whereas  a  list  of  values  can  provide  more 
information  (e.g.,  number  of  significant  digits)  about 
specific  data  points.  The  user  should  be  able  to 
interactively  select  from  various  formats  to  suit  his 
particular  needs. 

2.6  User  Interaction 

Generally  speaking,  the  computer  cannot  serve  as  a 
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useful  tool  unless  a  good  system  of  communication  has 
been  established  between  the  user  and  the  computer.  The 
purpose  of  this  section  is  to  present  the  major  issues 
involved  in  designing  a  good  communication  link.  Both 
psychological  and  technical  aspects  should  be  considered, 
as  with  any  other  human  engineering  problem  (Ref  9:2). 
Basically,  the  computer  must  be  adapted  to  the  user. 

The  OS,  as  previously  discussed,  is  an  important  first 
step  in  providing  a  degree  of  transparency  for  the  user. 
However,  this  only  applied  to  system  resources.  A 
higher  level  of  transparency  is  required  for  application 
programs  (e.g.,  DSP)  so  that  they  will  also  be  user- 
oriented.  Higher  level  just  means  that  applications  pro¬ 
grams  are  scheduled  by  the  OS. 

The  issues  concerning  user  interaction  presented 
are:  the  user,  response  time,  language,  input,  and 

user  assistance.  Even  a  system  which  addresses  all  of 
these  issues  with  some  degree  of  success  should  make 
provision  for  user  feedback  and  subsequent  changes 
(Ref  22:61).  The  system  designer  should  either  take  the 
responsibility  for  this  or  provide  sufficient  documen¬ 
tation  for  the  system  manager  to  take  appropriate  action. 
With  adherence  to  these  criteria,  the  ultimate  goal  of 
user  satisfaction  can  be  attained. 

2.6.1  The  User 

In  order  to  design  a  DSP  laboratory  for  user  sat¬ 
isfaction,  the  user  must  be  defined.  As  stated  previously. 
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the  users  will  be  students  and  faculty  members  at  AFIT. 

This  presents  the  problem  of  designing  a  system  for  a  group 
of  users  with  diverse  backgrounds  (Ref  29:17).  One 
requirement  is  that  the  user  should  not  need  a  computer 
background.  Users  will  also  have  varying  objectives 
when  using  the  laboratory.  The  spectrum  of  objectives 
can  be  defined  as  ranging  from  simply  entering  commands 
and  accompanying  parameters  for  performing  DSP  experi¬ 
ments,  to  modification  and  upgrade  of  the  system's  capa¬ 
bilities.  The  issues  of  user  interaction  which  are 
addressed  here  pertain  to  users  with  the  former  objective. 

2.6.2  Response  Time 

User  satisfaction  regarding  response  time  after  a 
command  has  been  entered  depends  both  on  the  user's 
expectations  and  the  system's  capabilities  (Ref  22:73). 
Expectations  can  be  modified  initially  by  providing 
approximate  execution  time  values  in  the  documentation. 

Of  course,  expectations  will  also  change  through  ex¬ 
perience  in  using  the  laboratory.  This  section  will 
focus  only  on  design  issues  pertinent  to  minimizing 
response  time.  Response  time  in  this  context  is  defined 
as  that  time  between  the  user  initiating  an  action 
(e.g.,  retrieval  of  a  data  file,  computations )  to  the 
system's  completion  of  that  action  (if  possible),  and 
notifying  the  user  that  it  is  completed  or  that  it 
cannot  be  completed. 

As  stated  previously,  the  system  is  intended  to 
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support  only  one  interactive  user  at  a  time.  Therefore, 
the  problems  associated  with  minimizing  response  time 
for  multiple  users  in  a  time-sharing  environment  are 
not  discussed. 

With  this  constraint,  the  focus  of  attention  is  on 
speed  of  execution.  The  speed  depends  on  the  hardware, 
software  (including  the  OS),  and  the  algorithms  employed. 
Issues  pertinent  to  execution  speed  are  far  too  exten¬ 
sive  to  adequately  cover  in  this  report.  The  interested 
reader  is  referred  to  the  literature. 

2.6.3  Language 

The  language  that  the  user  uses  to  communicate  with 
the  system  for  the  purpose  of  initiating  DSP  operations 
is  called  the  command  language  (CL).  The  CL  should 
resemble,  as  much  as  possible,  human  communication 
(Ref  9:27).  This  is  because  the  design  process  has 
the  goal  of  adapting  the  system  to  humans,  rather  than 
the  reverse.  Therefore,  commands  should  be  English 
words  instead  of  codes.  For  example,  the  command 
"REPLACE"  is  much  clearer  than  "RP".  However,  once 
the  user  gains  familiarity  v/ith  the  system  it  is  con¬ 
venient  to  have  the  option  of  entering  a  shortened  ver¬ 
sion  of  each  command  to  reduce  the  number  of  input 
characters  required. 

The  other  half  of  designing  for  human  communication 
involves  the  computer's  responses  to  user  requests.  Both 
requests  for  input  by  the  computer  and  notification  of 


the  system's  status  should  be  in  a  readable  English 
format . 

2.6.4  Input 

The  discussion  of  user  input  in  this  section 
assumes  an  interactive  environment,  although  most  of 
the  requirements  also  apply  to  batch.  These  require¬ 
ments  are:  user  interruption  of  a  requested  action, 
interception  of  user  input  errors  by  the  system,  correc¬ 
tion  of  input  errors  by  the  user,  provision  of  default 
values,  logical  sequencing  of  inputs,  construction  of 
a  sequence  of  user  requests,  and  useful  error  messages 
(Ref  17). 

Once  the  user  has  caused  the  system  to  begin  some 
action,  he  should  have  the  means  of  interrupting  the 
process.  This  feature  is  useful  for  recovering  from 
user  mistakes  which  are  not  unacceptable  errors  to  the 
system.  Those  errors  which  are  invalid  (e.g.,  exceed 
parameters,  illegal  commands)  should  be  in  intercepted 
by  the  system  before  execution  is  attempted.  The  user 
should  be  notified  at  that  point  of  the  specific 
error  and  be  given  options  for  further  action.  One  of 
the  options  should  allow  the  user  to  correct  the 
erroneous  input.  To  aid  in  avoiding  errors  in  the  first 
place,  default  values  should  be  provided,  as  appropriate. 
Another  way  to  avoid  input  errors  is  for  the  system 
to  request  the  input  in  a  logical  order.  Finally, 
the  user  should  be  able  to  construct  a  sequence  of 
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commands  (commonly  called  a  MACRO)  so  that  he  doesn't 
have  to  wait  for  each  action  to  complete  prior  to 
entering  the  next  command. 

2.6.5  User  Assistance 

Both  on-line  and  off-line  assistance  are  necessary 
for  a  user-responsive  system.  On-line  assistance  includes 
system  prompts  for  user  input,  lists  of  options,  and 
summary  information  for  each  valid  command.  Prompts  are 
required  so  that  the  user  will  always  be  aware  of  the 
nature  of  the  input  requested.  A  list  of  options  should 
always  follow  an  error  message  (Ref  22:112),  and  may 
also  be  used  with  prompts.  Summary  information,  often 
called  HELP  files,  should  give  the  user  a  brief  descrip¬ 
tion  of  each  valid  command  and  the  limits  for  each  of 
its  parameters  (Ref  17). 

Off-line  assistance  refers  to  the  documentation 
provided  for  the  DSP  laboratory.  A  system  reference 
manual  should  be  included  which  contains  comprehensive 
information  regarding  all  of  the  system's  capabilities 
not  covered  in  the  other  system  documents.  A  training 
manual  should  also  be  included.  It  should  teach  the 
user  how  to  properly  use  every  available  command  through 
the  use  of  examples  and  exercises  (Ref  22:56). 
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Chapter  3  Environment 

The  laboratory  environment  can  be  divided  into 
four  separate  categories:  the  operating  system  (OS), 
data  acquisition,  data  storage,  and  output  display.  These 
four  together  provide  the  capabilities  and  impose  the 
ultimate  limitations  under  which  the  user  must  work. 

At  the  beginning  of  this  project,  the  specific 
hardware  environment  included  an  HP  2108A  computer,  HP 
7906  disk  drive,  HP  2648A  graphics  terminal,  Remex 
high  speed  paper  tape  reader,  teletype  terminal  with 
paper  tape  punch/reader,  two  A/D  and  one  D/A  converter, 
memory  boards,  and  various  I/O  interfaces. 

The  graphics  terminal  was  not  initially  configured 
because  the  interface  board  for  it  had  not  yet  arrived. 

The  disk  drive  couldn't  be  used  at  first  because  a 
loader  ROM  had  not  arrived.  Without  the  disk  drive, 
the  operating  system  couldn't  be  generated  for  the  de¬ 
sired  configuration.  However,  all  the  hardware  listed 
above  (except  the  teletype)  is  now  configured,  with 
the  arrival  of  the  required  components  (Fig.  1).  The 
teletype  has  been  replaced  with  a  Texas  Instruments  (TI) 
Silent  700  printing  terminal.  The  OS  has  been  generated 
to  support  all  but  the  A/D  and  D/A  converters.  RTE-III 
does  not  include  system  drivers  for  these  devices. 

3.1  Operating  System 

RTE-III  was  chosen  as  the  OS  because  it  is  the 
latest  version  of  RTE  that  can  operate  with  the  currently 
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Fig.  1  Block  Diagram  of  Configured  Hardware 


available  support  hardware  ( RTE  stands  for  Real  Time 
Executive)  (Ref  26:1-1).  The  unique  hardware  required 
for  RTE- III  is  the  Dual  Channel  Port  Controller  (DCPC- 
used  for  direct  memory  access),  dynamic  memory  mapping, 
memory  protect,  and  a  minimum  of  32K  words  of  memory. 

First,  the  generation  of  RTE-III  is  discussed. 

This  should  assist  future  users  v/ho  need  to  generate 
a  new  OS.  Next,  the  capabilities  and  limitations  of 
RTE-III  are  discussed  with  respect  to  the  DSP  applica¬ 
tion  . 

3.1.1  Generation  of  RTE-III 

The  procedures  for  generating  a  new  OS  begin  by 
running  a  program  called  "RT3GN".  The  program,  as  it 
is  running,  will  ask  for  various  parameters  to  be 
entered  from  the  system  console.  Since  over  a  hundred 
responses  are  required,  the  preferred  method  of  entering 
the  parameters  is  to  create  an  answer  file.  See  Appendix 
A  for  a  listing  of  the  answer  file  used  in  the  most 
recent  generation. 

Generating  an  OS  takes  about  two  hours  if  every 
entry  in  the  answer  file  is  valid.  Unfortunately, 
the  computer  usually  halts  before  the  generation  is 
complete.  The  point  at  which  it  halts  appears  to  be 
random.  Why  this  anomaly  occurs  is  presently  unknown. 

The  current  version  of  the  OS  required  50  to  60  attempts 
before  it  finally  completed  (using  the  uame  answer 
file  every  time). 
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Once  the  OS  has  been  generated,  a  program  called 

“ SWTCH "  is  run  to  switch  from  the  current  OS  to  the 

new  OS.  These  procedures  are  summarized  below  and  differ 

somewhat  from  those  in  the  Operator's  Manual  (Ref  26). 

For  example,  the  Operator's  Manual  refers  to  the  program 

RTGEN  instead  of  RT3GN  ( RTGEN  may  have  been  an  earlier 

version  of  RT3GN).  Also,  the  Operator's  Manual  does 

not  discuss  the  use  of  an  answer  file. 

*RUN , FMGR 
: RP , RT3G1 
: RP , RT3G2 
: RP , RT3G3 
: RP , RT3G4 
: RP , RT3G5 
: RP , RT3G6 
: RP , RT3G7 
: RP , RT3GN 
:  EX 

*RUN , RT3GN 

LIST  FILE? 

L81914 

ECHO? 

TR, A81914 

(The  OS  and  listing  file  are  generated) 


RT3GN  FINISHED 
* RUN, SWITCH 

The  procedure  for  using  SWTCH  is  summarized  in  Appendix 


Once  the  OS  was  executing,  the  first  task  was  to 
create  a  backup  disk  so  that  the  OS  could  be  recovered 
in  the  event  of  a  disk  failure.  An  off-line  utility 
called  DSKUP  was  used  to  copy  the  lower  platter  (where 
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the  OS  resides)  to  the  upper  platter  (which  is  removable). 
DSKUP  was  loaded  from  the  right  cassette  drive  of  the 
graphics  terminal  (See  Appendix  F).  The  procedure  for 
using  DSKUP  is  shown  in  Appendix  B. 

One  error  has  been  discovered  in  the  most  recent 
answer  file.  Three  memory  partitions  were  specified 
as  18K,  instead  of  the  maximum  of  17K,  as  reported  by 
RT3GN  in  the  list  file.  This  didn't  cause  any  genera¬ 
tion  errors.  These  partitions  are  assumed  to  be  17K 
each,  so  that  3K  words  of  physical  memory  are  unavailable 
for  use. 

3.1.2  Capabilities  and  Limitations 

RTE-III  provides  useful  software  support  for  pro¬ 
gram  development  and  testing.  Included  are  an  inter¬ 
active  editor,  an  assembler,  and  compilers  for  BASIC, 

ALGOL,  FORTRAN,  and  FORTRAN  IV.  The  only  compiler  that 
has  been  used  is  FORTRAN  IV  because  DSP  programs  taken 
from  the  literature  were  written  in  this  language.  How¬ 
ever,  two  limitations  should  be  noted  in  light  of  the 
DSP  application.  The  primary  limitation  is  the  size 
of  logical  memory  (See  Ref  26  for  a  discussion  of  physical 
memory  vs.  logical  memory).  DSP  programs  are  typically 
quite  large  (the  FIR  filter  design  program  in  Appendix 
D  required  15K  words  of  user  partitioned  memory),  so 
that  many  programs  cannot  be  used  on  this  system,  unless 
they  are  segmented.  This  is  because  the  maximum  parti¬ 
tion  size  for  user  programs  is  17K  words.  The  maximum 


partition  size  is  set  at  generation  time,  where  RT3GN 
provides  the  maximum  value. 

Another  notable  limitation  is  that  writing  an  ex¬ 
ecutive  program  which  can  reference  any  of  the  DSP  pro¬ 
grams  is  not  as  simple  as  it  might  be.  The  simplest 
way  to  write  an  executive  program  is  to  write  it  in 
FORTRAN  IV  (or  another  high  order  language)  and  load  it 
into  memory  with  all  of  its  subprograms.  This  method 
is  not  practical  for  the  HP  system  because  logical 
memory  size  is  relatively  limited.  Alternate  solutions 
are  discussed  in  chapter  five. 

3.2  Data  Acquisition 

As  stated  previously,  data  can  be  acquired  in  real¬ 
time  from  an  A/D  converter,  or  with  software  that  gen¬ 
erates  data.  This  section  focuses  on  data  obtained 
which  represents  real-time  phenomena.  Examples  of  soft¬ 
ware  generated  data  can  be  found  in  chapter  four.  Two 
related  subjects  are  presented  here:  The  A/D  hardware 
and  its  software  support. 

3.2.1  A/D  Converter 

Two  A/D  converters  were  provided  for  this  project. 
They  are  the  ADC-12QZ  model,  made  by  Analog  Devices 
(Ref  6:175-173).  Each  is  mounted  on  an  AC-1548  mounting 
card.  An  input  range  of  -lOv  to  +10v,  represented  by 
12  bits  in  twos-complement  format  is  configured.  The 
bipolar  range  was  chosen  to  allow  more  flexibility  in 
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the  type  of  signals  that  can  be  analyzed.  Twos-comple- 
ment  format  was  chosen  because  this  is  how  the  HP 
computer  represents  numbers  in  memory.  Each  A/D  conver¬ 
ter  is  connected  to  an  HP  12566B  micro-circuit  inter¬ 
face  card  (Ref  16)  as  shown  in  Fig.  2. 

3,2.2  Software  Support  for  A/D  Converter 

Four  different  kinds  of  software  support  have  been 
investigated  for  the  A/D  converter  (Ref  34).  The  first 
is  polling  (i.e.,  check  the  I/O  channel  flag  bit)  soft¬ 
ware,  whereby  the  computer  initiates  an  A/D  conversion, 
waits  in  a  loop  for  the  conversion  to  complete,  loads 
the  data  into  a  register,  and  then  stores  it  in  a 
memory  buffer.  The  second  type  is  direct  memory  access 
(DMA).  In  this  case,  the  computer  initializes  one 
channel  of  the  DCPC,  which  will  perform  the  required  oper¬ 
ations  for  data  acquisition  and  storage.  The  third 
method  is  to  use  a  system  driver.  This  method  as  only 
partially  been  implemented.  The  fourth  method  involves 
the  use  of  interrupts.  The  interrupt  method  was  found 
to  be  inappropriate  for  this  application.  Each  of 
these  methods  is  discussed  with  regard  to  its  respective 
advantages,  disadvantages,  and  implementation.  All  re¬ 
levant  program  1  istings  are  in  Appendix  E. 

3. 2. 2.1  Polling 

The  main  advantage  of  polling  software  is  that  a 
software  delay  loop  can  be  included  to  allow  the  user 
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to  interactively  specify  the  sampling  rate.  In  the  sub¬ 
routine  WAIT  1 ,  the  delay  is  increased  by  increasing  the 
variable  DELAY.  Of  course,  DELAY  is  just  an  integer 
which  will  cause  a  different  amount  of  delay  for  differ¬ 
ent  computers  (even  of  the  same  type)  and  different 
operating  conditions  (e.g.,  temperature).  In  order  to 
relate  the  delay  value  to  time,  a  test  subroutine  (AD14 
or  AD15)  was  implemented.  This  subroutine  is  nearly 
identical  to  the  original  (AD14B  or  AD15B),  so  that  the 
speed  of  both  is  approximately  the  same  in  the  data  ac¬ 
quisition  loop.  The  speed  of  the  instructions  in  this 
loop  will  determine  the  amount  of  time  from  the  completion 
of  one  conversion  to  the  beginning  of  the  next  conversion. 
This  time  -  the  instruction  execution  time  (IET)  -  and 
the  conversion  time  (CT)  can  be  measured  with  an  oscillo¬ 
scope  by  checking  the  status  bit  of  the  A/D  converter 

(Fig.  2,  3).  The  sum  of  these  two  times  is  the  sample 

time  (ST).  Several  measurements  were  taken  for  different 
values  of  DELAY,  and  are  summarized  in  Table  I.  The 
measurements  are  also  plotted  in  Fig.  4.  In  each  case, 
the  CT  was  about  33  us. 

The  theoretical  sample  time  (TST)  for  a  periodic 
signal  of  a  single  frequency  can  be  expressed  as 
TST  =  NC/(fxNP) 

where 

NC  =  number  of  cycles 

f  =  frequency  in  hertz 

NP  =  number  of  sampled  data  points 


_ _ 
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DELAY 


Fig.  4  Plotted  Sample  Time  Measurements 


Table  I 


Sample  Time  Measurements 


DELAY 

ST,  us 

0 

76 

10 

130 

30 

240 

50 

350 

80 

510 

100 

620 

150 

890 

200 

1150 

500 

2775 

Once  the  TST  is  calculated,  the  number  for  DELAY  can  be 
selected  from  Fig.  4.  Data  can  then  be  acquired  and 
plotted  using  this  value  for  DELAY  and  the  value  of  NP 
selected  above.  The  actual  NC  plotted  is  used  to  cal¬ 
culated  ST.  By  comparing  ST  with  TST,  error  in  relating 
DELAY  to  time  can  be  determined. 

The  disadvantages  of  polling  software  are  twofold. 
First,  it  is  slower  than  DMA  because  data  must  first  be 
transferred  to  a  register.  The  effect  of  this  is  to 
lower  the  maximum  frequency  that  can  be  accurately 
represented.  Secondly,  this  method  precludes  real-time 
processing,  whereby  one  data  buffer  is  being  filled 
while  another  is  being  processed  (Ref  32).  Instead, 
the  processor  fills  one  buffer  until  it  is  full. 

3. 2. 2. 2  Interrupts 

Performing  I/O  by  the  use  of  interrupts  is  inappro¬ 
priate  for  this  application  for  several  reasons.  First, 
if  a  system  driver  is  not  used,  I/O  for  the  A/D  conver¬ 
ter  can  only  be  accomplished  with  the  interrupt  system 
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turned  off.  Otherwise,  an  illegal  (not  allowed  by  the  OS) 
interrupt  error  will  occur.  Secondly,  because  of  the 
high  speed  synchronous  nature  of  the  A/D  converter,  an 
interrupt  must  be  handled  within  a  certain  amount  of  time. 
RTE-III  has  many  potential  interrupts  from  other  I/O 
channels  that  it  must  handle  by  priority,  so  there  is  no 
guarantee  that  data  will  not  be  lost  from  the  A/D  con¬ 
verter.  Thirdly,  data  which  is  acquired  must  first  be 
transferred  to  a  register  as  with  the  polling  method, 
so  that  about  the  same  software  overhead  is  incurred. 
Therefore,  the  interrupt  method  offers  no  speed  advantage 
over  the  polling  method.  If  the  interrupt  method  could 
be  used,  it  would  have  the  advantage  that  some  useful 
computations  could  be  performed  between  interrupts. 

3. 2. 2. 3  Direct  Memory  Access 

The  fastest  kind  of  data  acquisition  available  for 
this  system  is  direct  memory  access.  The  acquisition 
rate  is  almost  that  of  the  A/D  conversion  rate.  The 
entire  transfer  process  is  handled  by  hardware  (the  DCPC) 
once  it  has  been  initialized  by  the  appropriate  instruc¬ 
tions.  Therefore,  no  software  overhead  is  incurred  in 
acquiring  and  storing  data.  Additionally,  the  processor 
is  free  to  perform  useful  computations  in  parallel  with 
data  acquisition.  However,  the  processor  speed  is  reduced 
slightly  because  the  processor  cannot  reference  memory 
at  the  same  instant  of  time  that  the  DCPC  is  transferring 
data  to  or  from  memory.  Thus,  the  number  of  computations 


during  data  acquisition  can  be  increased  by  reducing 
the  number  of  memory  reference  instructions.  The  capa¬ 
bility  of  simultaneous  data  acquisition  and  processing 
permits  real-time  processing,  which  is  discussed  in 
chapter  four. 

The  only  disadvantage  of  DMA  is  that  the  sampling 
rate  cannot  be  specified  interactively  by  the  user,  with 
out  additional  hardware.  This  hardware  might  include  a 
programable  timer,  which  is  initialized  prior  to  the 
transfer  of  each  data  sample.  A  count-finished  signal 
would  then  be  substituted  for  the  status  signal  in  Fig. 

The  subroutine  DMA14  (or  DMA15)  in  Appendix  E 
demonstrates  the  software  needed  to  use  DMA.  It  is 
called  from  the  main  program  ANDG2. 

3. 2. 2. 4  I/O  Driver 

A  driver  is  a  program  designed  to  support  I/O  for 
a  particular  device.  It  is  placed  in  the  system  area 
at  generation  time,  and  gives  the  OS  total  control 
over  I/O.  Although  RTE-III  does  not  provide  drivers  for 
the  D/A  converter  card  and  the  microcircuit  interface 
card,  user  written  drivers  can  be  included  in  the  OS 
generation  (Ref  27).  Unfortunately,  in  debugging  a 
driver,  each  new  version  requires  another  generation 
of  the  OS.  This  is  a  lengthy  process  even  if  the  genera 
tion  completes  every  time  with  the  first  attempt. 

In  order  to  permit  user-written  drivers  to  be  im¬ 
plemented  and  debugged  without  generating  a  new  OS, 


29 


three  dummy  drivers  have  been  included  as  part  of  the 
latest  OS  generation.  By  themselves,  these  drivers  do 
nothing  but  return  to  their  calling  program.  However, 
by  changing  the  return  address  in  the  system  area  to 
the  beginning  address  of  a  user-supplied  driver  before 
calling  a  dummy  driver,  the  actual  driver  can  be  tested. 

A  large  buffer  area  is  included  in  the  dummy  driver  so 
that  the  actual  driver  can  eventually  be  copied  into  this 
area.  This  method  has  not  been  explored  beyond  the 
conceptual  stage. 

Note  that  the  OS  will  assign  a  DCPC  channel  for  the 
driver.  If  a  DCPC  channel  is  not  available,  the  device 
goes  into  a  "waiting  for  DCPC"  state  (Ref  27:3-16).  Be¬ 
cause  of  the  possibility  of  waiting,  drivers  are  probably 
not  suitable  for  real-time  processing. 

3.3  Data  Storage 

The  laboratory  has  two  types  of  devices  to  meet  the 
requirements  of  data  storage  from  main  memory:  the 
cassette  drives  in  the  graphics  terminal  and  the  disk 
drive.  The  merits  and  some  possible  applications  of 
each  are  presented. 

3.3.1  Disk 

The  disk  drive  can  be  used  in  two  different  ways  to 
store  data.  One  method  is  useful  for  temporary  storage 
while  a  program  is  executing.  The  other  method  can  be 
used  for  permanent  storage  and  later  retrieval.  Both 
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methods  can  be  useful  in  DSP  applications.  Each  is 
demonstrated  by  programs  listed  in  Appendix  G. 

Storing  data  temporarily  on  a  disk  involves  the  Disk 
Track  Allocation  EXEC  Call  (Ref  26:3-18)  and  the  system 
utility  B1NRY  (Ref  11:3-2).  With  this  utility,  the  exact 
area  on  the  disk  where  the  data  transfer  is  to  occur  must 
be  specified.  RTE-III  will  assign  the  tracks  to  the  program, 
and  release  the  tracks  when  the  program  terminates.  FMGR 
and  other  programs  cannot  access  data  stored  in  this 
manner  once  the  tracks  have  been  released. 

One  specific  application  of  BINRY  is  to  pass  data 
between  two  concurrently  running  programs  which  are  in 
separate  partitions.  In  addition  to  the  data,  some  kind 
of  message  must  also  be  transferred  to  signal  that  the 
data  transfer  is  complete  or  in  progress.  This  applica¬ 
tion  has  not  been  implemented. 

For  permanent  storage,  FMP  calls  are  used.  They 
allow  a  user  program  to  create  a  disk  file,  and  read 
from  or  write  into  a  disk  file.  The  data  will  remain 
accessible  to  other  programs  (such  as  FMGR)  after  the 
user  program  has  terminated.  To  transfer  data  among  a 
large  system  of  DSP  programs,  a  standard  format  for  a 
data  file  should  be  established.  This  would  make  it 
easier  for  a  program  to  find  the  desired  data  within  any 
given  data  file. 

3.3.2  Cassette 

One  way  to  store  program  generated  data  into  a  disk 
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file  is  to  first  store  it  on  a  cassette  tape.  Although 
the  transfer  rate  is  not  as  high  as  that  of  the  disk, 
data  stored  on  a  cassette  tape  can  be  accessed  by  FMGR 
without  the  need  for  a  user  program  to  create  a  disk 
file.  Before  the  data  is  sent  to  the  tape  ( LU  4  or  LU  5) 
tape  positioning  characters  should  be  sent.  After  com¬ 
pleting  the  data  transfer,  an  end-of-file  (EOF)  mark 
should  be  written  on  the  tape.  These  special  characters 
can  be  sent  by  using  the  method  described  in  chapter  five 

The  data  can  be  transferred  to  a  disk  file  with  a 
FMGR  command  as  follows: 

ST , 4 , DATA1 : : 20 

In  this  case,  the  file  is  transferred  from  the  left 
cassette  to  a  file  called  DATAl ,  which  is  created  by  FMGR 
on  cartridge  20  (upper  platter  disk)  (Ref  1:2-23). 

3.4  Output  Display 

Six  types  of  output  display  can  be  supported  with 
the  system  in  its  current  configuration  using  the  fol¬ 
lowing  devices:  ccnventional  CRT  terminal,  line  printer, 
graphics  terminal,  graphics  printer,  oscilloscope,  and 
X-Y  plotter  (the  printers  must  be  connected  directly  to 
the  graphics  terminal).  A  line  printer  or  graphics 
printer  can  be  connected  to  the  computer  if  the  RTE-III 
printer  driver  is  included  in  a  new  OS  generation. 

3.4.1  Conventional  CRT  Terminal  and  Line  Printer 

These  types  of  displays  are  best  suited  for  listing 
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tables  of  values  because  no  special  manipulation  of  the 
data  is  required.  Graphs  and  charts  can  also  be  imple¬ 
mented,  but  require  additional  instructions  to  draw  axes, 
label  axes,  scale  the  data,  etc.  The  CRT  has  the  advan¬ 
tage  of  speed  for  checking  results  quickly.  The  printer 
has  the  advantage  of  providing  a  permanent  copy  of  useful 
results . 

3.4.2  Graphics  Terminal/Printer 

The  graphics  terminal  and  graphics  printer  displays 
are  the  best  method  of  implementing  graphs  and  charts 
because  the  resolution  is  only  limited  by  the  size  of 
a  single  point.  However,  the  display  size  is  limited  to 
one  screen  for  the  HP  2648A  terminal.  A  graphics 
printer  would  simply  provide  an  exact  copy  of  the  screen 
display.  The  graphics  terminal  and  graphics  printer  can 
also  serve  the  purpose  of  a  conventional  terminal  and 
printer,  respectively.  A  graphics  printer  is  not  cur¬ 
rently  configured  because  none  is  available. 

The  HP  2648A  terminal  has  a  very  useful  feature  for 
exploiting  the  graphics  capability,  called  Auto  Plot 
(Ref  15).  The  only  requirements  for  using  it  are  to 
fill  out  the  Auto  Plot  menu,  which  specifies  the  data 
parameters,  and  then  provide  the  data  (usually  sent  from 
the  computer).  These  procedures  are  implemented  entirely 
by  the  computer  in  an  example  in  chapter  five. 

3.4.3  D/A  Converter 

The  D/A  converter  used  for  this  project  is  an  HP 
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12555A.  It  has  two  analog  outputs  corresponding  to  X 
and  Y  axes.  Each  16-bit  word  in  memory  that  is  transferred 
is  separated  into  two  bytes.  Each  byte  is  then  converted 
to  an  analog  signal,  ranging  from  0  to  10  volts.  The 
D/A  converter  also  has  a  20  msec  refresh  timer  for  con¬ 
ventional  oscilloscopes,  and  an  unblanking  signal  for 
storage-type  oscilloscopes  {Ref  13).  The  converter 
is  also  suitable  for  driving  an  X-Y  plotter. 

3.4.4  Software  Support  for  D/A  Converter 

The  program  that  has  been  implemented  to  support 
the  D/A  converter  will  work  for  either  a  conventional 
oscilloscope  or  an  X-Y  plotter  (See  Appendix  I).  The 
inputs  to  the  subroutine  DAI 6B  are  two  real  arrays, 
one  containing  the  X  values,  and  the  other  containing 
the  Y  values.  DA16B  scales  the  data  into  bytes,  for¬ 
mats  the  data  into  one  array,  and  then  outputs  it  to 
the  D/A  converter. 

The  main  differences  between  the  oscilloscope 
approach  and  the  plotter  approach  are  the  refresh  timer, 
which  is  required  for  the  conventional  oscilloscope, 
and  the  delay  loop,  which  matches  the  speed  of  the  com¬ 
puter  to  that  of  the  plotter.  Both  devices  are  cur¬ 
rently  connected  in  parallel  to  the  X  and  Y  outputs  of 
the  D/A  converter. 


Chapter  4  Computational  Tasks 

Pour  DSP  applications  programs  have  been  implemented 
and  tested.  These  programs  were  selected  from  the  liter¬ 
ature  (Ref  7)  after  examining  many  such  programs  (Ref  3, 
5,  21,  31).  The  modifications  made  to  allow  them  to  run 
on  the  AFIT  HP  computer  were  minor  and  will  not  be 
discussed.  See  Ref  7  for  a  detailed  discussion  of  the 
algorithms  involved.  The  four  programs  that  were  imple¬ 
mented  were  chosen  for  two  reasons.  First,  they  are 
written  in  ANSI  standard  FORTRAN  IV,  so  that  few  changes 
were  expected.  The  final  versions  are  listed  in 
Appendix  D.  Secondly,  in  each  case,  the  author(s)  of  the 
program  included  a  test  problem.  By  running  the  test 
problem  on  the  HP  system  and  comparing  that  output  to 
the  output  provided  in  the  test  problem,  it  was  possible 
to  determine  if  the  program  was  not  functioning  properly. 
The  outputs  generated  for  all  four  programs  by  the  HP 
system  were  very  close  (only  a  few  differences  in  the  two 
least  significant  digits)  to  that  of  the  test  problem. 

The  differences  are  assumed  to  be  due  to  the  different 
accuracies  of  the  computers  on  which  the  programs  were 
run.  Based  on  this  assumption,  the  four  programs  are 
functioning  properly.  This  chapter  concludes  with  a  dis¬ 
cussion  of  a  method  for  implementing  real-time  processing 
( RTP ) .  RTP  involves  simultaneous  data  acquisition,  data 
processing,  and  display  of  results. 
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4.1  FFT/IFFT 

Many  variations  of  the  original  radix-2  Cooley- 
Tukey  algorithm  (Ref  19)  have  been  implemented.  This 
algorithm  is  commonly  called  the  fast  Fourier  transform 
( FFT ) .  The  algorithm  also  computes  the  inverse  of  the 
transform  and  is  called  the  inverse  fast  Fourier  trans¬ 
form  (IFFT).  The  motivation  for  using  the  FFT  is  that  it 
substantially  reduces  the  required  number  of  machine  com¬ 
putations  for  calculating  the  DFT.  An  algorithm  used 

2 

prior  to  the  FFT  involved  N  complex  multiplications  and 
2 

approximately  N  additions  (N  is  the  number  of  sample 
points).  On  the  other  hand,  the  basic  FFT  requires  at 
most  (l/2)x  Nlog2N  complex  multiplications,  where  N  is 
an  integral  multiple  of  two  (Ref  5:275). 

As  an  example  of  the  usefulness  of  the  FFT,  consider 
the  case  where  N=8192.  Approximately  2.5  minutes  of 
computation  time  is  required  to  compute  the  DFT  on 
an  IBM  370  using  an  algorithm  which  preceeded  the  FFT. 

For  the  same  data,  the  FFT  requires  approximately  1/4 
second  (Ref  5:264). 

The  FFT  is  an  essential  capability  for  a  general 
purpose  DSP  laboratory.  It  offers  a  computationally 
efficient  time-to-f requency  transformation.  It  is  also 
useful  for  other  DSP  computations  such  as  convolution, 
correlation,  power  spectral  density,  and  digital  filter 
design  (Ref  6).  The  program  FAFT4  (See  Appendix  D) 
calls  a  short  demonstration  version  of  the  FFT  (i.e., 
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FOURE).  FAFT4  generates  data  lor  a  complex  function, 
computes  the  theoretical  DFT  ( TDFT ) ,  computes  the  FFT 
(using  FOURE),  computes  the  IFFT  (based  on  values  obtained 
from  the  TDFT,  using  FOURE),  and  outputs  these  four  sets 
of  data.  The  program  has  been  modified  to  also  list 
these  sets  of  data  in  a  format  suitable  for  plotting 
with  Auto  Plot. 

FAFT4  is  not  intended  to  be  part  of  the  applications 
software  for  the  laboratory.  A  more  general  use  of  FOURE 
should  be  provided  which  allows  the  user  to  interactively 
specify  the  data,  the  number  of  data  points,  and  whether 
an  FFT  or  an  IFFT  is  to  be  computed. 

4.2  FIR  Filter  Design 

Digital  filtering  is  an  important  part  of  DSP  and 
is  in  common  use.  A  very  important  type  of  filter  is 
the  linear  phase  finite  impulse  response  digital  filter 
(Ref  5:185).  One  implementation  of  this  is  the  program 
FIRF2  (See  Appendix  D).  It  uses  the  Remex  exchange 
algorithm  with  minimum  weighted  Chebyshev  error  in 
approximating  a  desired  ideal  frequency  response.  The 
inputs  are:  1)  filter  length  (in  samples),  2)  type  of 
filter  (multiple  passband/stopband ,  differentiator,  or 
Hilbert  transformer),  3)  number  of  frequency  bands, 
specified  by  upper  and  lower  cutoff  frequencies, 

4)  desired  frequency  response  in  each  band,  5)  positive 
weight  function  in  each  band,  and  6)  grid  density.  The 
data  (provided  by  the  authors)  to  test  the  program  was 
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in  the  following  format: 

55,1,5,0 

0.0,0.05,0.1,0.15,0.18,0.25,0.3,0.36,0.41,0.5 

0.0, 1.0, 0.0, 1.0, 0.0 

10. 0,1.0,3.0,1.0,20.0 

Data  is  read  by  the  program  from  the  left  cassette  (LU4) 
of  the  graphics  terminal.  For  the  DSP  laboratory  the 
user  should  be  able  to  interactively  specify  the  data 
source  and  enter  the  data  from  the  keyboard,  if  that  is 
the  source. 

4.3  Autocorrelation  and  Covariance 

The  program  COVAU  implements  the  autocorrelation 
and  covariance  methods  of  linear  prediction  (LPC)  anal¬ 
ysis.  It  consists  of  three  parts:  subroutine  AUTO 
(autocorrelation),  subroutine  COVAF  (covariance),  and 
a  test  program  which  calls  the  two  subroutines  and 
generates  the  test  input  data.  AUTO  implements  a  form 
of  Robinson's  recursion.  COVAR  implements  a  form  of 
Cholesky  decomposition.  For  the  DSP  laboratory,  the 
user  should  be  able  to  interactively  specify  the  data 
source  and  method  of  LPC  analysis. 

4.4  HR  Filter  Design 

The  purpose  of  program  FWIIR  is  to  design  finite 
word-length  IIR  digital  filters.  It  consists  of 
eight  parts:  a  main  program,  which  acts  as  an  executive, 
and  seven  subroutines.  Subroutine  INDAT  reads  in  all 
the  problem  data.  FUNCT  computes  the  value  of  the  maxi¬ 
mum  weighted  error  function  on  the  coarse  grid,  and  at 


the  end,  on  a  finer  grid.  HANDJ  is  a  randomized  version 
of  the  Hooke  and  Jeeves  Pattern  Search.  EXPLR  is  a  local 
exploration  program  called  by  HANDJ .  SET  sets  one 
vector  equal  to  another.  SHUFF  randomly  orders  the  list 
of  coordinates  for  EXPLR.  UNI  generates  a  random  number. 

For  the  DSP  laboratory,  the  user  should  be  able  to 
interactively  enter  the  required  data,  or  specify  another 
data  source. 

4.5  Real-Time  Processing 

In  this  context,  real-time  processing  refers  to 
simultaneous  data  acquisition,  data  processing,  and 
display  of  results,  without  any  time  limit  or  interruptions 
of  these  three  functions.  One  basic  assumption  of  RTP 
is  that  the  data  can  be  processed  and  results  displayed 
at  least  as  fast  as  it  is  acquired.  Otherwise,  some  data 
would  be  overwritten  in  memory  by  new  data  before  the 
old  data  was  processed.  In  DSP  applications,  data  is 
typically  processed  in  time  frames  (Ref  32).  Each  time 
frame  is  represented  by  a  buffer  in  main  memory.  If  a 
time  frame  is  processed  as  a  group  of  data,  then  a  buffer 
cannot  be  processed  while  it  is  being  filled.  With  this 
constraint,  at  least  two  buffers  are  required. 

RTP  was  beyond  the  scope  of  this  project,  but  a 
possible  way  to  implement  it  on  the  HP  system  is  shown 
in  a  flowchart  (See  Fig.  5).  The  program  represented  by 
the  flowchart  is  terminated  by  pressing  the  "CLEAR  DISPLAY" 
button  on  the  front  panel  of  the  computer.  This  method 


of  program  termination  has  been  successfully  used  with 
another  program  (Appendix  1).  Each  number  in  the  flow¬ 
chart  corresponds  to  one  of  the  steps  listed  below.  SR 
is  the  switch  register  displayed  on  the  front  panel. 

1 .  Set  SR  to  1 

2.  Is  DCPC  channel  1  ready  (flag  bit  set)? 

3.  Is  DCPC  channel  2  ready  (flag  bit  set)? 

4.  Begin  DMA  data  acquisition  in  memory  buffer  #1 

5.  Begin  DMA  data  acquisition  in  memory  buffer  #2 

6.  Process  the  data  in  buffer  #1 

7.  Format  the  data  in  buffer  #1  for  display 

8.  Begin  DMA  output  from  buffer  01 

9.  Process  the  data  in  buffer  #2 

10.  Format  the  data  in  buffer  #2  for  display 

11.  Begin  DMA  output  from  buffer  #2 

12.  Is  SR=0? 

Some  computations  may  be  too  lengthy  to  permit  real¬ 
time  processing.  One  possible  solution  to  this  problem 
is  to  replace  portions  of  an  assembly  language  or  high- 
order  language  program  with  a  microprogram.  This  involves 
identifying  the  "bottle-neck"  portions  of  the  program 
and  replacing  them  with  microprogram  subroutines  which 
are  more  time  efficient.  The  system,  in  its  current  con¬ 
figuration,  supports  this  capability.  See  Ref  30  for  a 
complete  discussion  of  the  potential  for  improving 
program  execution  speed  through  microprogramming. 


41 


Chapter  5  User  Interaction 

Several  types  of  software  support  were  implemented 
to  make  the  DSP  laboratory  more  user-responsive,  or  to 
demonstrate  its  potential  in  this  regard.  Chapter  three 
discussed  the  generation  of  the  OS  and  the  software 
support  provided  for  the  A/D  and  D/A  converters.  This 
chapter  first  presents  a  program  which  demonstrates 
computer  control  of  the  graphics  terminal's  functions. 

This  is  followed  by  a  discussion  of  a  method  for  executing 
a  program  which  exceeds  the  logical  memory  size  of  the 
largest  available  user  partition  (17K  words).  In  this 
way,  the  user  is  not  constrained  to  command  the  execution 
of  several  smaller  programs  to  accomplish  the  same  task. 
Next,  two  methods  are  presented  whereby  one  program  can 
schedule  another  program.  A  scheduling  facility  is 
useful  for  implementing  a  DSP  executive  program  which 
can  schedule  the  DSP  applications  programs.  Finally,  a 
software  structure  is  recommended  for  a  large  system  of 
DSP  programs. 

5.1  Graphics  Terminal  Control 

Since  most  user  interaction  will  occur  with  the 
graphics  terminal,  an  effort  was  made  to  relieve  the 
user  of  many  of  the  manual  operations  needed  to  operate 
the  terminal.  The  program  APLOT  listed  in  Appendix  H 
demonstrates  a  method  by  which  any  graphics  terminal 
operation  that  can  be  performed  manually  by  the  user 
can  also  be  initiated  by  the  computer.  This  involves 
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the  use  of  escape  codes  (Ref  15).  The  decimal  equivalent 
of  the  ASCII  character  codes  are  sent  from  the  computer 
using  the  Read/Write  EXEC  call  (Ref  26:3-4).  APLOT  is 
written  in  FORTRAN  IV,  but  the  characters  can  also  be 
sent  from  an  assembly  language  program.  The  program  fills 
in  the  Auto  Plot  menu,  draws  the  axes  and  grid  and  plots 
the  program  generated  data  sent  from  the  computer. 

5.2  Segmented  Program 

A  segmented  program  consists  of  a  memory-resident 
main  program  and  one  or  more  disk-resident  subprograms 
(i.e.,  segments).  By  using  a  segmented  structure,  the 
total  program  size  is  only  limited  by  the  available 
disk  storage  space.  Data  is  passed  among  the  segments 
and  the  main  program  through  a  common  storage  area 
within  the  memory  partion  allocated  for  the  program. 

In  spite  of  the  size  advantage  of  segmented  programs, 
several  limitations  should  be  noted.  First,  the  program 
will  execute  slower  than  if  the  entire  program  resided 
in  main  memory,  because  a  segment  must  be  loaded  from 
the  disk  before  it  can  be  executed.  Therefore,  care  must 
be  exercised  to  insure  that  two  or  more  frequently 
referenced  sections  of  a  program  are  not  chosen  to  be 
disk-resident  segments.  Secondly,  the  only  way  to  return 
control  to  the  main  program  from  a  segment  is  for  the 
segment  to  call  the  main  program  (this  applies  to  FORTRAN 
IV  programs,  not  assembly  language).  However,  this  will 
cause  the  main  program  to  start  executing  from  the 
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beginning.  Therefore,  conditional  branching  to  the  seg¬ 
ments  should  occur  at  or  near  the  beginning  of  the 
main  program.  A  third  limitation  of  a  segmented  program 
is  that  if  one  segment  calls  another  segment,  the  calling 
segment  is  overlayed  by  the  segment  it  calls. 

Writing  a  segmented  program  involves  the  use  of  the 
program  segment  load  EXEC  call  (Ref  26:3-23).  An  example 
of  this  is  presented  in  Appendix  J.  The  following  FMGR 
commands  were  entered  to  load  the  program  for  execution: 

:MR, %MAIN1 
: MR , %SEGN1 
: MR , %SEGM2 
:MR, %SEGM3 

: RU , LOADR , 9  9 , 1 , ,00001 

To  save  the  loaded  program  for  future  execution,  the 

following  commands  were  entered: 

:SP,MAIN1 
: SP , SEGM1 
: SP , SEGM2 
: SP , SEGM3 

This  saves  the  program  as  a  type  6  file  (ready  for  execu¬ 
tion)  on  the  disk. 

5.3  Executive  Programs 

The  two  programs  discussed  in  this  section  are 
called  executive  programs  because  they  can  schedule  other 
programs,  each  of  which  resides  in  a  separate  disk  file. 
Both  are  listed  in  Appendix  K.  One  method  of  implementing 
a  DSP  executive  program  (program  SCHED)  involves  the 
program  schedule  EXEC  call  (Ref  26:3-24)  and  a  program 
suspend  EXEC  call  (Ref  26:3-22).  Data  is  passed  through 
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the  system  area  buffer  with  a  Class  I/O  -  Read/Write 
EXEC  call  (Ref  26:3-6)  and  a  Class  I/O  -  Get  EXEC  call 
(Ref  26:3-13).  Another  method  of  implementing  a  DSP 
executive  program  (program  MESGl)  involves  the  Message 
Processor  utility,  MESSS  (Ref  26:4-40).  Any  of  the 
system  commands  can  be  executed  from  a  FORTRAN  IV  program 
by  using  MESSS  (e.g.,  RU  -  run  command).  Data  is 
passed  in  the  same  manner  described  for  the  previous 
method . 

5.4  Recommended  Software  Structure 

The  software  structure  for  the  DSP  laboratory  (dem¬ 
onstrated  by  program  EXECT)  should  be  comprised  of  a 
user-oriented  executive  program  and  a  system  of  DSP 
applications  programs.  The  type  of  executive  demonstrated 
by  the  program  SCHED  is  considered  to  be  more  suitable 
for  the  DSP  application  than  that  demonstrated  by  MESGl. 
SCHED  gives  the  user  the  option  of.  immediate  scheduling 
so  that  the  output  of  the  scheduled  program  can  be 
seen  before  control  is  returned  to  the  executive.  This 
is  not  the  case  for  MESGl.  EXECT  uses  the  SCHED  method. 

In  addition,  several  of  the  programs  which  are  scheduled 
by  EXECT  are  segmented.  Therefore,  EXECT  demonstrates 
the  system's  potential  for  a  user-oriented  executive 
program  which  can  schedule  any  one  of  a  large  system  of 
DSP  applications  programs. 
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Chapter  6  Recommendations 

The  system,  in  its  present  state,  does  not  qualify 
as  a  user-responsive  general  purpose  DSP  laboratory, 
which  is  the  ultimate  goal.  This  is  because  not  all  of 
its  current  capabilities  have  been  exploited,  and  be¬ 
cause  additional  capabilities  need  to  be  added.  Both 
of  these  issues  will  be  addressed  in  this  chapter.  Im¬ 
plementation  of  the  recommendations  which  follow  should 
result  in  a  system  which  meets  the  design  goal  as 
stated  above. 

6.1  Additional  Hardware 

The  hardware  recommended  in  this  section  would  bring 
the  laboratory  much  closer  to  being  general  purpose. 

It  is  intended  to  support  educational  and  research 
activities  at  AFIT.  It  is  also  intended  to  support 
research  for  laboratories  at  WPAFB  which  might  solicit 
AFIT  students  for  thesis  projects  on  other  systems. 
Therefore,  a  diversity  of  equipment  is  desireable.  How¬ 
ever,  only  the  addition  of  a  graphics  printer  is  consider¬ 
ed  essential . 

6.1.1  Magnetic  Tape  Drive 

The  magnetic  tape  drive  can  be  useful  for  at  least 
two  reasons.  First,  information  stored  on  disk  can  be 
lost  in  the  event  of  a  disk  failure.  The  information 
could  be  later  retrieved  if  it  is  first  backed  up  on 
tape.  Secondly,  a  tape  can  be  used  to  transfer  large 
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amounts  of  information  from  one  computer  system  to 
another  (not  necessarily  HP).  For  certain  types  of 
research,  this  may  be  the  most  efficient  means  of  doing 
so . 

6.1.2  Faster  A/D  Converter 

The  A/D  converters  which  are  currently  configured 
have  a  measured  conversion  time  of  about  33  us.  Accord¬ 
ing  to  the  sampling  theorem,  a  minimum  of  two  samples/ 
cycle  of  the  input  signal  are  required  to  reconstruct 
the  signal.  This  corresponds  to  a  maximum  frequency  of 
about  15.15  kHz,  which  does  not  even  cover  all  of  the 
audio  spectrum  (20Hz  -  20kHz).  This  limits  the  types 
of  research  related  to  spectrum  analysis  that  can  be 
conducted  on  the  HP  system. 

Although  HP  makes  an  A/D  interface  (model  91000A) , 
the  sample  rate  is  only  20,000  samples/second.  By  the 
above  criterion,  the  maximum  useful  frequency  is  10kHz. 
Since  this  is  more  limited  than  the  existing  parts,  the 
HP  part  is  not  recommended.  A  better  alternative  is  the 
ACD1103  A/D  converter,  made  by  Analog  Devices.  It  per¬ 
mits  three  different  ccnversion  times:  1)  8  bits  in 

1  us,  2)  10  bits  in  1.5  us,  or  3)  12  bits  in  3.5  us 
(Ref  6).  These  are  maximum  values.  The  comparative 
maximum  useful  frequencies  are  500  kHz,  333  kHz,  and 
143  kHz  respectively. 

Only  the  conversion  time  specification  has  been 
considered  above  in  recommending  an  A/D  converter. 
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However,  since  some  applications  may  emphasize  other 
characteristics  such  as  accuracy  and  linearity,  they 
should  also  be  considered. 

6.1.3  HP  Fourier  Analyzer 

A  hardware  solution  to  real-time  processing  is  the 
HP  Fourier  Analyzer.  It  has  an  oscilloscope  display 
and  performs  operations  such  as  the  FFT/IFFT,  convolution, 
and  correlation.  Several  of  the  HP  users  at  WPAFB  have 
this  device.  Conversations  with  some  of  the  users  indi¬ 
cated  that  it  is  very  useful  for  real-time  DSP. 

6.1.4  Graphics  Printer 

A  graphics  printer  (GP)  is  an  essential  addition  to 
the  DSP  laboratory.  It  can  provide  a  hard  copy  for 
several  types  of  output.  For  example,  displays  from 
the  graphics  terminal's  memory  (i.e.,  graphics  memory) 
can  be  sent  directly  to  the  GP.  This  output  is  limited 
to  one  screen  size.  The  GP  can  also  serve  as  a  line 
printer.  Any  type  of  output  that  can  be  sent  to  a  line 
printer  can  also  be  sent  to  the  GP.  From  the  author's 
experience,  the  GP  should  operate  at  a  serial  data 
transfer  rate  cf  at  least  300  band  in  order  to  qualify 
as  user-responsive. 

6.2  DMA  for  D/A  Converter 

DMA  is  the  preferred  method  of  sending  data  to  the 
D/A  converter.  It  can  be  implemented  either  as  a  system 
driver  or  in  a  user  program.  For  RTP ,  DMA  is  essential 
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for  a  real-time  display.  No  software  has  yet  been 
written  to  provide  this  capability. 

6.3  Additional  Applications  Programs 

The  four  DSP  programs  which  have  been  implemented 
(See  chapter  four)  certainly  are  not  sufficient  for  the 
diverse  computational  applications  intended  for  the 
laboratory.  These  programs  were  published  by  IEEE 
(Ref  7),  and  are  available  on  standard  1/2  inch  tape 
from  IEEE,  along  with  the  other  DSP  programs  in  the 
text.  DSP  programs  can  also  be  found  in  many  other 
textbooks  or  periodicals.  All  of  the  computational 
tasks  listed  in  section  2.4  should  be  implemented. 

6.4  Real-Time  Processing 

Even  if  an  HP  Fourier  analyzer  is  configured,  it 
may  be  desired  to  perform  computations  which  the  an¬ 
alyzer  does  not  provide  for.  The  method  for  RTP 
described  in  chapter  five,  or  some  equivalent,  should 
be  implemented.  RTP  is  a  common  activity  in  DSP  and 
would  make  the  laboratory  more  general  purpose.  There¬ 
fore,  it  is  an  essential  capability. 

6.5  Operating  System 

RTE-III  has  made  a  significant  contribution  to  the 
establishment  of  the  HP  system  as  a  DSP  laboratory. 
However,  because  of  the  limitations  previously  discussed, 
RTE-IVB  is  recommended  as  an  upgrade  OS.  It  allows 
larger  memory  partitions  (e.g.,  28K)  for  user  programs 
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to  be  specified  (Ref  28:VI-13).  It  also  has  a  feature 
called  Extended  Meir.ory  Area  (EMA).  The  EMA  is  a 
large  area  of  memory  within  a  partition,  limited  only 
by  the  size  of  physical  memory  (Ref  28:VI-14).  These 
are  just  two  of  the  many  advantages  of  an  RTE-IVB  OS  up¬ 
grade.  To  do  this  will  require  the  HP  2108A  processor 
board  to  be  modified  or  replaced.  An  OS  upgrade  is  not 
considered  essential  because  RTE-III  has  been  demonstrated 
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Appendix  A:  RTE-III  Generation  Answer  File 

The  generation  answer  file  provides  responses  to 
input  requests  by  RT3GN.  In  this  way,  the  generation 
of  the  OS  can  proceed  without  interruption  for  inputs 
which  otherwise  would  be  entered  from  the  system  console. 
Listed  in  this  appendix  is  the  answer  file  used  in  the 
most  recent  OS  generation.  Changes  for  a  new  OS 
generation  can  be  made  by  simply  editing  the  file; 
the  filename  is  A81914. 
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♦  «  OF  IC  CLADCECC 

♦  "  LU  MA"C'  I ‘I'P 

♦  :s  PM 

♦  BUFFER  LIMIT C  CLOW  ?  H I GA'* 


13  *  DVR 33*  D 

1 3  *  DVR  05  * E • X* 1 3  *  T= 1  3  0 0  0 
10*  DVR 36 

c.'  O  *  D  V  P  0  O  *  t'  *  T  =  1  '3  0  0  0 
17*  DVP01 • E*  T =500 
4*  DVP43 

1  4  *  D  V  P  5  O  *  D  *  T  =  1  o  ij 


♦  EOT 

♦  EOT 

♦  EOT 

♦  EOT 

♦  EOT 

♦  EOT 

♦  EOT 
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15.  DVC'51 .  D.  T=1  00 
lS-DVESS.  I'.  1=1  00 


*  DEVICE  REFERENCE  TABLE 

♦ 

2'  0 
1  .  1 
0.  0 


0  «  o 
o,  0 
On  0 

o.  o 

4  >  0 
5- .  0 

0,  0 
0'  0 
ij<  0 

^ 

0«  0 
7.  0 
3*  0 
5>  0 
1  >  0 
1  <  2 
Ch  i'i 

o«  o 

"•  <  0 


♦  INTERRUPT  TAELE 

» 

4.  ENT«  lP0i,ie 
IS-EST-  1 


SO .  PPG.  pc  MPT 

0 

0 

i‘  — .. 

ve  : 

£0 

1.4. EG 
2 .  T .  E  G 
3.13.  EG 
4*14. EG 


*  EOT 

*  ZST  =:■? 


♦  LU1  SYS.  CONSOLE 

♦  LU2  SYS.  DISK  < LOWER 

♦  LUS  HUY.  DISC 

♦  LU4  L.  CTU  DRIVE 

♦  LUS  P.  CTU  DRIVE 

♦  LUG 

♦  LU7 

♦  LUS 

♦  LUS 

♦  LU10  TTY 

♦  LU11  RARE5,  TAPE  cEh] 

♦  LU12 

♦  LU 1  3 

♦  _U  1 4 

♦  LUIS  WCS  NODULE  12'rl 

♦  LUIS 

♦  LU17  MICROCIRCUIT  1. 

♦  LUIS  MICROCIRCUIT  1* 

♦  LUIS  D  A  CONVERTER 

♦  L'J£0  Did  SUBCHANNEL 

♦  LU2 1  DISC  SUBCHANNEL 

♦  LU22 

♦  LUS  3 

♦  LJ24 

♦  _U2S  ROWE”'  “AIL 


EDITP.  14 
LOfiDP • 14 
^SME.  1 

:  :pef.  i 

,-'E 
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Appendix  B:  Using  SWTCH  and  DSKUP 

SWTCH  is  a  program  which  is  run  to  switch  operating 
systems.  Below  are  listed  the  questions  that  the  pro¬ 
gram  will  ask  and  the  responses  that  were  entered  for 
the  last  generation.  SP  means  space;  CR  means  carriage 
return . 

FILE  NAME  OF  NEW  RTE  SYSTEM? 

S81914 

TARGET  CHANNEL  FOR  NEW  SYSTEM 

SP  CR 

TARGET  SUBCHANNEL  ( LOGICAL ) /UNIT  FOR  NEW  SYSTEM 

SP  CR 

NOW  IS  THE  TIME  TO  INSERT  CARTRIDGE  IN  TARGET 

SUBCHANNEL  UNIT. 

SP  CR 

SAVE  FILES  AT  TARGET? 

Y 

PURGE  TYPE  6  FILES 

N 

AUTO  BOOT-UP? 

Y 

READY  TO  TRANSFER.  OK  TO  PROCEED? 

Y 

DSKUP  is  an  off-line  utility  which  can  be  loaded 
from  a  cassette  tape.  To  load  it  from  the  left  cas¬ 
sette,  set  the  S  Reg.  to  041300.  Press  PRESET,  IBL, 
and  RUN.  When  it  is  finished  loading,  set  the  S  Reg. 
to  13  and  the  P  Reg.  to  2.  Press  PRESET  and  RUN.  Below 
are  listed  the  questions  that  the  program  will  ask  and 
the  responses  that  were  given  to  copy  the  operating 
system  to  a  removable  disk. 


DISC  BACKUP  UTILITY 
TASK? 

CO 

SOURCE  DISC  CHANNEL#  ? 

12 

SOURCE  DISC  TYPE? 

7906 

SOURCE  DISC  DRIVER#? 

O 

TYPE  OF  COPY? 

FR 

RTE  OR  DOS  DISC? 

RT 

FROM  CYLINDER#? 

O 

#  OF  TRACKS? 

256 

#  OF  SURFACES? 

2 

STARTING  HEAD#? 

2 

DEST  DISC  DRIVER#? 

O 

TO  CYLINDER#? 

O 

#  OF  SURFACES? 

2 

STARTING  HEAD#? 

O 

6144  WORD  BUFFER  DESIRED? 
YES 

VERIFY? 

YES 
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Appendix  C:  HP  users  at  WPAFB 


The  people  listed  below  were  visited  to  determine 
if  they  had  any  DSP  software  currently  in  use,  which 
might  be  useful  to  this  project.  Another  reason  for 
the  visits  was  to  observe  their  hardware  configurations 
with  respect  to  support  of  DSP  applications.  Unfortunately, 
no  useful  software  was  found.  Those  involved  with  DSP 
used  an  HP  Fourier  Analyzer. 


NAME 

PHONE 

BUILDING 

ORGANIZATION 

John  Bankovskis 

56361 

622 

AFWAL/AARI 

Bill  Griffin 

54016/52789 

20 

ASD/ENADA 

Jim  Leonard 

53050 

23 

AFWAL/AAFR/2 

L.  T.  Drzal 

52952 

32 

AFWAL/MLBM 

Bob  Ballard 

52493 

24C 

AFWAL/FIMN 

Irvin  F.  Luke 

52372 

18 

AFWAL/POOC/1 
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Appendix  D:  Computational  Tasks  Listings 

Four  source  code  DSP  program  listings  are  contained 
in  this  Appendix.  They  were  taken  from  the  literature 
(Ref  7)  and  modified  to  run  on  the  HP  system.  Also 
included  are  the  inputs  used  to  test  each  program  (if 
needed),  and  the  resultant  outputs.  FAFT4  computes 
the  FFT  and  IFFT  of  a  complex  function.  FIRF2  designs 
linear  phase  FIR  filters.  CCVAU  implements  autocorrela¬ 
tion  and  covariance  methods  of  linear  prediction  analysis 
FWIIR  designs  finite  word-length  HR  digital  filters. 

See  Ref  7  for  a  detailed  explanation  of  the  algorithms 
employed. 


n  o  n  >-  noo  onnoonn  onnoonooo  nnnonnonnnn  on 


PTN4,L 


PROGRAM  FAFT4 


MAIN  PROGRAM-  FOURSUBT 
AUTHOR:  C.  M.  RADER 

MIT  LINCOLN  LABORATORY,  LEXINGTON,  MA  02173 

INPUT:  FUNCTION  IS  GENERATED  BY  THE  PROGRAM  TO  TEST 

SUBROUTINE  FGURE 


COMPLEX  W,  C ( 32 ) ,  D,  E,  B<32),  QBC32),  A,  F(32) 

SEQUENCE  LENGTH  N=2#*MU  USES  MU=S  THUS  N=32 

METHOD  IS  TO  COMPUTE  DFT  OF  KNOWN  FUNCTION  A*#I ,  1=0 , 1 , . . . ,N-t 
OUTPUT  OF  PROGRAM  IS  LARGEST  DIFFERENCE  BETWEEN  DFT  COMPUTED  TWO 
WHEN  RUN  ON  AN  IBM  370,  THE  MAX  DIFFERENCE  WAS  0.3B9E-04. 

WHEN  RUN  ON  A  HONEYWELL  608QN,  THE  MAX  DIFFERENCE  WAS  Q.238E-06. 
A  MAX  DIFFERENCE  LARGE  COMPARED  TO  THE  ACCURACY  OF  THE  COMPUTER 
WOULD  PROBABLY  INDICATE  A  PROGRAM  ERROR. 

MU  =  S 

SET  UP  MACHINE  CONSTANT 

IOUTD  =  1 1MACH < 2 )  < - THIS  STATEMENT  IS  USELESS  TO  HP. 

SET  IOUTD  TO  THE  SYSTEM  CONSOLE  LU# . 


IOUTD  =  i 

NN  =  2**MU 

TP  I  =8 . #ATAN( 1 . ) 

TP  ION  =  TPI/FLOAT (NN) 

W  =  CMPLX<COS<TPION) ,-SIN<TPION) ) 

GENERATE  A¥*K  AS  TEST  FUNCTION 

A  =  (  9, .3) 

B  ( 1 )  =  (i.,0.) 

QB(  1 )  =  B  < 1 ) 

DO  10  K=2 , NN 

B(K)  =  A##CK-i> 

QB( K )  =  B(K) 

0  CONTINUE 

PRINT  COMPLEX  INPUT  SEQUENCE 
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WRITE  ( IOUTD , 9999 ) 

9999  FORMAT  ( IX , //////24H  COMPLEX  INPUT  SEQUENCE  ) 

WRITE  < IOUTD , 9998  >  ( I , QB ( I ) , 1= 1 , NN ) 

9998  FORMAT  (2(2X,  1H(  ,  13,  1H) ,  2E14.6)) 

C 

C  B(i)  CONTAINS  A**0  B(K)  CONTAINS  A*#K-lj  ETC 
C 

C  COMPUTE  DFT  OF  B  IN  CLOSED  FORM 
C 

D  =  <i  .  ,0  .  )  -  A*#NN 
DO  20  K= i , NN 

E  =  < 1  .  , 0  .  )  -  A*U#*<K-i> 

C  <  K )  =  D/E 
F(K )  =  C  <  K ) 

20  CONTINUE 
C 

C  DFT  OF  B  IS  (i-A#¥NN)/(l-AW*#K  ) 

C 

C  NOW  COMPUTE  DFT  OF  B  USING  FOURE 
C 

CALL  FOURE (B,  NN,  -1) 

C 

C  PRINT  FOURE  DFT  AND  THEORETICAL  DFT 
C 

WRITE  (IOUTD, 9997) 

9997  FORMAT  ( IX , //////l 1H  FOURE  DFT  ) 

WRITE  (IOUTD, 9998)  ( I , B ( I ) , I =1 , NN ) 

WRITE  (IOUTD, 9996) 

9996  FORMAT  ( IX ,//////17H  THEORETICAL  DFT  ) 

WRITE  (IOUTD, 9998)  ( I , C ( I ) , 1=1 , NN ) 

C 

C  FIND  MAX  DIFFERENCE  BETWEEN  B  AND  C 
C 

DD  =  0  . 

DO  30  1=1, NN 

DE  =  CABS ( C ( I ) -B ( I ) ) 

IF  (DD.GT.DE)  GO  TO  30 
DD  =  DE 

30  CONTINUE 
C 

C  COMPUTE  INVERSE  DFT  OF  C>  (F=C> 

C 

CALL  FOURE(F ,  NN ,  1) 

C 

C  PRINT  INVERSE  DFT 
C 

WRITE  (IOUTD, 999S) 

9995  FORMAT  ( IX , //////19H  FOURE  INVERSE  DFT  ) 

WRITE  (IOUTD, 9998)  ( I ,F( I ) , 1=1 ,NN) 

C 

C  COMPUTE  MAX.  DIFF  BETWEEN  INPUT  AND  INVERSE  OF  THEOR .  DFT 
C 


64 


GG  =  Q 

DO  40  1=1, NN 

GG1  =  CABS(QB( I >-F( I) ) 

IF  (GGi.LE.GG)  GO  TO  40 
GG  =GG 1 

40  CONTINUE 
C 

C  PRINT  MAXIMUM  DIFFERENCE 
C 

URITE  ( IOUTD , 9994 )  DD 
WRITE  (IOUTD, 9993)  GG 

9994  FORMAT  </42H  MAX  DIFF  BETWEEN  THEOR  AND  FOURE  DFT  IS  ,  E12.3) 
9993  FORMAT  (/SiH  MAX  DIFF  BETWEEN  ORIGINAL  DATA  AND  INVERSE  DFT  IS  , 
*  E12 . 3 , ////// ) 

URITE  (4,9000)  ( QB ( I ) , B ( I ) , C ( I > , F ( I > , 1=1 , NN > 

9000  FORMAT (8(2X,F7  S) > 

STOP 

END 

C 

C - 

C  SUBROUTINE  FOURE 

C  PERFORMS  COOLEY-TUKEY  FAST  FOURIER  TRANSFORM 

C - 

C 

SUBROUTINE  FOURE( DATA ,  N,  ISI > 

C 

C  THE  COOLEY-TUKEY  FAST  FOURIER  TRANSFORM  IN  ANSI  FORTRAN 
C 

C  DATA  IS  A  ONE-DIMENSIONAL  COMPLEX  ARRAY  WHOSE  LENGTH,  N,  IS  A 
C  POWER  OF  TWO.  ISI  IS  +1  FOR  AN  INVERSE  TRANSFORM  AND  -1  FOR  A 
C  FORWARD  TRANSFORM.  TRANSFORM  VALUES  ARE  RETURNED  IN  THE  INPUT 
C  ARRAY,  REPLACING  THE  INPUT. 

C  TRANSFORM(J)=SUM(DATA(I.)*W**((I-l)*(J-i))),  WHERE  I  AND  J  RUN 
C  FROM  1  TO  N  AND  U  =  EXP  (ISI*2*PI*SQRT(-1)/N> .  PROGRAM  ALSO 
C  COMPUTES  INVERSE  TRANSFORM,  FOR  WHICH  THE  DEFINING  EXPRESSION 
C  IS  INVTR(J)=(1/N#SIJM(DATA(I)#W#*((I-1)#(J-1  )))  . 

C  RUNNING  TIME  IS  PROPORTIONAL  TO  N*L002(N>,  RATHER  THAN  TO  THE 
0  CLASSICAL  N#*2. 

C  AFTER  PROGRAM  BY  BRENNER,  JUNE  1967  THIS  IS  A  VERY  SHORT  VERSION 
C  OF  THE  FFT  AND  T.S  INTENDED  MAINLY  FOR  DEMONSTRATION  PROGRAMS 
C  ARE  AVAILABLE  IN  THIS  COLLECTION  WHICH  RUN  FASTER  AND  ARE  NOT 
C  RESTRICTED  TO  POWFRS  OF  2  OR  TO  ONE-DIMENSIONAL  ARRAYS 
C  SEE  —  IEEE  TRANS  AUDIO  (JUNE  1967),  SPECIAL  ISSUE  ON  FFT 
C 

COMPLEX  DATA(N) 

COMPLEX  TEMP ,  U 

C  IOUTD  =  I 1MACH ( 2 )  ( THIS  STATEMENT  USELESS  TO  HP. 

C 

C  SET  IOUTD  TO  SYSTEM  CONSOLE  LU* . 

C 

IOUTD  =  1 
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c 

C  CHECK  FOR  POWER  OF  TWO  UP  TO  15 
C 

NN  =  1 

DO  10  1=1,15 
M=  I 

NN  =  NN*2 

IF  (NN.EQ.N)  GO  TO  20 
10  CONTINUE 

WRITE  ( IOUTD , 9999 ) 

9999  FORMAT  (30H  N  NOT  A  POWER  OF  2  FOR  FOURE  > 

STOP 

20  CONTINUE 

C 

PI  =  4.*ATAN<1. ) 

FN  =  FLOAT (N ) 

C 

C  THIS  SECTION  PUTS  DATA  IN  BIT-REVERSED  ORDER 
C 

J  =  i 

DO  80  1=1 ,N 
C 

C  AT  THIS  POINT,  I  AND  J  ARE  A  BIT  REVERSED  PAIR  (EXCEPT  FOR  THE 
C  DISPLACEMENT  OF  +1) 

C 

IF  (I-J)  30,  40,  40 
C 

C  EXCHANGE  DATA(I)  WITH  DATA(J)  IF  I.LT.J. 

30  TEMP  =  DAT  A  <  J ) 

DATA(J)  =  DATA  < I ) 

DATA(I)  =  TEMP 
C 

C  IMPLEMENT  J=J+1,  BIT-REVERSED  COUNTER 
C 

40  rt  =  N/2 

50  IF  (J-M)  70,  70,  60 

6  0  J  =  J  -  M 

M  =  (M+D/2 
GO  TO  50 
70  J  =  J  +  M 

80  CONTINUE 

C 

C  NOW  COMPUTE  THE  BUTTERFLIES 
C 

MMAX  =  1 

90  IF  (MMAX-N)  100,  130,  130 

100  ISTEP  =  2*MMAX 

DO  120  M=i,MMAX 

THETA  =  P IKcFLOAT  (ISI*(M-l))/FLOAT(  MMAX  > 

W  =  CMP  LX (COS (THETA) ,SIN(THETA>  > 


66 


DO  iiO  I=rt , N , ISTEP 
J  =  I  +  HMAX 
TEMP  a  U*DATA<J) 

DATA< J)  =  DATA< I )  -  TEMP 
DATA(I)  =  DAT  A  < I)  +  TEMP 
IiO  CONTINUE 
120  CONTINUE 

MMAX  =  ISTEP 
GO  TO  90 

130  IF  (ISI)  160 ,  140,  140 
C 

C  FOR  INV  TRANS  —  ISI=1  —  MULTIPLY  OUTPUT  BY  1/N 
C 

140  DO  150  1  =  1, N 

DATA(I)  =  DATA<I)/FN 
150  CONTINUE 
160  RETURN 
END 


FAFT4  Original  Output 


COMPLEX  INPUT  SEQUENCE 


( 

i) 

. IQOOOQE+Oi 

. QQQQQOE+QQ 

< 

3) 

. 720000E+00 

. 54Q  QQQE+Q  0 

< 

5) 

.22680QE+00 

. 77760  QE+Q  0 

( 

7) 

- . 256608E+QQ 

. 682344E+0  Q 

( 

9) 

- . 553223E+Q  0 

. 352719E+Q0 

< 

11) 

588789E+00 

-.447828E-Q1 

< 

13) 

- . 399745E+0  0 

- .350190E+00 

( 

15) 

- . 987 1 42E-0 1 

- . 467999E+Q  0 

( 

17) 

. 181645E+Q  0 

- .390265E+00 

( 

19) 

. 341S28E+0  0 

-.182902E+&0 

( 

21) 

. 344667E+0  0 

. 527352E-Q1 

( 

23) 

. 219683E+Q0 

.224090E+00 

< 

25) 

. 371635E-0 1 

. 279973E+0Q 

( 

27) 

124428E+00 

.  22.1649E+QQ 

( 

29) 

- . 209279E+QQ 

. 923963E-0 1 

< 

31 ) 

- . 200575E+00 

464851E-01 

FOURE  DFT 

< 

1) 

.  693974E+0  0 

.  349971E+Q 1 

< 

3) 

.  94Q297E+0 1 

- . 91350QE+01 

( 

5) 

.  1 13182E+01 

223416E+01 

< 

7) 

.  799557E+0 0 

1 13961E+0 1 

< 

9) 

.700862E+00 

- . 699565E+Q0 

( 

11) 

.  6531 09E+0  0 

- . 446246E+0  0 

< 

13) 

.  623788E+00 

- . 272086E+Q  0 

< 

15) 

.  602883E+0  0 

1357Q3E+Q0 

< 

17)  . 

.  586277E+0  0 

- . 179496E-0 1 

( 

19) 

.  57 1899E+0 0 

. 926Q65E-0 1 

( 

21) 

.  558492E+00 

. 204881 E  +  0  0 

< 

23) 

.  S45133E+0  0 

. 323364E  +  0  0 

( 

25) 

.  5310  Q 1E  +  00 

. 476678E  +  0  0 

< 

27) 

.  515362E+00 

.674849E+00 

< 

29) 

.  499467E+0  0 

980906E+00 

< 

31) 

.  49Q733E+QQ 

15770  8E  +  0 1 

( 

2) 

. 90  0  0  0  QE  +  Q  0 

. 3000Q0E+00 

( 

4) 

. 4860  OOE  +  Q  0 

70200 OE +00 

( 

6) 

29160QE-01 

. 76788QE  +  0  0 

( 

8) 

435650E+0Q 

537127E+00 

( 

10) 

- . 6  037 17E  +  0  0 

15 1480E+Q  0 

< 

12) 

- . 516475E+Q  0 

- . 216941E+  0  0 

< 

14) 

- . 2547 14E  +  0  0 

- . 435094E  +  0  0 

( 

16) 

5 15569E-0 1 

450813E+00 

( 

18) 

. 28Q56QE+00 

- .296745E+00 

( 

20  > 

. 362245E+Q  0 

- .621538E-01 

( 

22) 

. 29438QE+0  0 

. 150862E+00 

( 

24) 

. 13048SE+Q  0 

. 267586E+  0  0 

( 

26) 

- . 505448E-0 1 

. 263125E  +  Q  0 

( 

28) 

178400E+00 

.162156E+00 

( 

30  > 

216070E+00 

. 203731E-0  1 

( 

32) 

.1 66572E+Q 0 

- . 1020  0  9E  +  0  0 

( 

2) 

. 279227E+0 1 

. 805Q45E+0 1 

( 

4) 

. 186645E+0 1 

- . 3833S3E+0 1 

< 

6) 

904795E+0  0 

- . 153463E+  0 1 

( 

8) 

. 739606E+0  0 

-  .  882314E  +  0  0 

( 

10) 

. 673575E+00 

- . 558478E +0  0 

< 

12) 

. 63699 1E+0  0 

- . 352639E+0  0 

( 

14) 

. 612613E+0  0 

- . 20  0642E  +  0  0 

( 

16) 

. 59420 1E+0  0 

- . 753133E-0  1 

( 

18) 

. 5789  0  0E  +  0  0 

. 376520E- 0 1 

< 

20  ) 

. 565136E+0  0 

. 1 47983E+0  0 

( 

22) 

. 551859E+00 

.  264522E+00 

( 

24) 

53821 4E+00 

. 393257E  +0  0 

( 

26) 

. 523403E  +  0  0 

.  567132E  +  00 

( 

28) 

. 506926E  +  0  0 

.8031 0 1E+  0  0 

( 

30  ) 

. 49 1389 E+0  0 

. 121921E+01 

( 

32) 

. S17354E  +  0  0 

. 218383E+0 1 
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THEORETICAL  DFT 


< 

i  ) 

. 693974E  +  0  0 

34997 1E+0 1 

( 

2) 

279227E+Q  1 

3Q5Q4SE+0 1 

< 

2) 

940298E+0  1 

- . 91350 1E+Q 1 

< 

4) 

13664SE+01 

-  383383E+0 i 

< 

S) 

. U3182E  +  01 

- . 223416E+Q 1 

( 

6) 

904794E+Q0 

-  .  153463E  +  0 1 

( 

7) 

. 799557E  +  0  0 

- . 1 1396 1E+0 1 

( 

8) 

739606E+0 0 

- . 382315E  +  0  0 

( 

9) 

. 700862E+00 

- . 698S66E+0Q 

( 

10  ) 

. 673576E  +  0  0 

-  553479E+  0 0 

< 

11  ) 

653110E+00 

-  .  44624SE+Q0 

( 

12) 

. 63699 1E+Q  0 

-  .  352639E+  0  0 

< 

13) 

. 623788E+00 

- . 272086E+Q0 

( 

14) 

. 612613E+00 

-  .200642E  +  0Q 

< 

IS) 

. 602883E  +  0  0 

-  135703E+QQ 

( 

16) 

. 59420  OE  +  Q  0 

-  .  7531 40E-0 1 

< 

17) 

. 586277E+0  0 

1 79496E-Q 1 

( 

18) 

. 57890  0E  +  0  0 

. 37651 4E- 0 1 

< 

19) 

. 571899E  +  0  0 

. 926066E-Q 1 

( 

20  ) 

. 565136E+0C 

. 147983E+0  0 

< 

21) 

. 55B492E  +  0  0 

. 2Q  4  880E+Q  0 

( 

22  ) 

. 551 859E+0  0 

. 264522E  +  0  0 

< 

23) 

. 545134E+0  0 

. 328364E+Q  0 

< 

24) 

. S38215E+0  0 

. 398257E+0  0 

< 

25 ) 

. 531002E+00 

. 476677E  +  Q  0 

< 

26) 

. 523404E  +  0  0 

. 567132E+00 

< 

27) 

. 51 5363E+  0  0 

. 674849E  +  0  0 

( 

23) 

. 506927E  +  0  0 

. 808100 E+00 

< 

29) 

498468E  +  0  0 

. 980  9Q5E  +  0  0 

( 

30) 

. 491391E+0  0 

12 1920E  +  0 1 

< 

31) 

490734E+Q0 

. 157708E+01 

( 

32) 

. 5 1 7358E+  0 0 

. 218883E  +  0  1 

FOURE 

INVERSE  DFT 

< 

1 ) 

.  lOOOOOE+Oi 

-.774860E-Q6 

( 

2) 

90  0  0  0  0E+0  0 

. 299999E+0  0 

< 

3) 

.  720  0 0 QE+Q 0 

.5400  0  0E+  0  0 

( 

4) 

.  436000E  +  00 

. 7020  0  0E+0  0 

< 

5) 

.  22680 OE+Q 0 

. 77760  QE+Q  0 

( 

6) 

.  291604E- 0 1 

. 7673S0E+  0  0 

< 

7) 

- 

.  2566QRE  +  0  Q 

632344E+00 

( 

8) 

- 

.  435651 E  +  0  0 

. 537127E+0  0 

< 

9) 

- 

.  5S3224E+0  0 

352719E+0Q 

( 

10  ) 

.  6  03717E+0  0 

. 151480E+00 

< 

11) 

.  S88790E+aa 

-  .  447831E-01 

( 

12) 

.  516476E+00 

- .216942E+00 

< 

13) 

- 

.  399746E+0Q 

- . 350 19QE+00 

( 

14) 

- 

.  254714E+00 

435095E+00 

< 

IS) 

.  93714QE-0 1 

-  467999E+00 

( 

16) 

.  515572E-01 

450814E+00 

< 

17) 

.  181646E  +  QQ 

- . 39Q265E+00 

( 

13) 

.  280  56 1E+  0  0 

-  .  296745E  *-00 

< 

19) 

.  341528E  +  0Q 

- . 132902E+QQ 

( 

20) 

.  362246E+00 

-.621535E-01 

< 

21) 

.  344667E+  0 0 

. 52735 9E-01 

( 

22) 

.  294330E  +  00 

. 150862E+00 

< 

23) 

.  219683E  +  0  Q 

. 224090E  +  0  0 

( 

24) 

.  130  488E+0  0 

. 267S86E+0 0 

< 

25) 

.  371634E-01 

. 279974E  +  0  0 

( 

26) 

.  5 05450E- 0 1 

. 263 126E+  0  0 

< 

27) 

.  124428E  +  QQ 

. 221649E+0  0 

( 

28) 

.  178430E  +  00 

. 162156E+0  0 

( 

29) 

209278E+Q0 

. 923962E-01 

( 

30  ) 

- 

.  216069E  +  00 

20  3729E-0 1 

( 

31) 

- 

.  200574E  +  0  0 

- . 464854E-01 

( 

32) 

.  1 66571 E+0  0 

- . 1 020  0  9E  +  0  0 

MAX  DIFF  BETWEEN  THEOR  AND  FOURE  DFT  IS  ,  . iilE-04 

MAX  DIFF  BETWEEN  ORIGINAL  DATA  AND  INVERSE  DFT  IS  105E-0S 
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FAFT4  Output  in  Auto  Plot  Format 
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o  n  o  n  n  n  o  o  n  o  o  n  o  n  o  o  n  o  n  n  r,  o  r>  n  n  n  n  n  n  n  n  o  n  n  n  o  n  o  n  o  n  n  rj  o  n  n  o  on 


FTN4  ,  L 


PROGRAM  FIRF2 


MAIN  PROGRAM.-  FIR  LINEAR  PHASE  FILTER  DESIGN  PROGRAM 
AUTHORS--  JAMES  H.  MCCLELLAN 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING  AND  COMPUTER  SCIE 
MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
CAMBRIDGE,  MASS.  02139 

THOMAS  U.  PARKS 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
RICE  UNIVERSITY 
HOUSTON,  TEXAS  77001 

LAURENCE  R  RABINER 

BELL  LABORATORIES 

MURRAY  HILL,  NEU  JERSEY  07974 


INPUT 

NFILT —  FILTER  LENGTH 
JTYPE —  TYPE  OF  FILTER 

1  =  MULTIPLE  PASSBAND/STOPBAND  FILTER 

2  =  DIFFERENTIATOR 

3  =  HILBERT  TRANSFORM  FILTER 
NBANDS--  NUMBER  OF  BANDS 

LGRID —  GRID  DENSITY,  WILL  BE  SET  TO  16  UNLESS 

SPECIFIED  OTHERWISE  BY  A  POSITIVE  CONSTANT. 

EDGE ( 2*NBANDS ) —  BANDEDGE  ARRAY,  LOUER  AND  UPPER  EDGES  FOR  EACH 
WITH  A  MAXIMUM  OF  10  BANDS. 

F'X(NBANDS) —  DESIRED  FUNCTION  ARRAY  (OR  DESIRED  SLOPE  IF  A 
DIFFERENTIATOR)  FOR  EACH  BAND 

WTX ( NBANDS ) —  WEIGHT  FUNCTION  ARRAY  IN  EACH  BAND.  FOR  A 

DIFFERENTIATOR,  THE  WEIGHT  FUNCTION  IS  INVERSELY 
PROPORTIONAL  TO  F. 

SAMPLE  INPUT  DATA  SETUP: 

32,1,3,0 

0.0,0.1,0.2,0.35,0. 425 ,0.5 
0  .  0 , 1  .  0 , 0  .  0 
10  0,1.0,10.0 

THIS  DATA  SPECIFIES  A  LENGTH  32  BANDPASS  FILTER  WITH 
STOPBANDS  0  TO  0.1  AND  0  425  TO  0  5,  AND  PASSBAND  FROM 
0.2  TO  0.35  WITH  WEIGHTING  OF  10  IN  THE  STOPBANDS  AND  1 
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IN  THE  PASSBAND  THE  GRID  DENSITY  DEFAULTS  TO  16 
THIS  IS  THE  FILTER  IN  FIGURE  10. 

THE  FOLLOWING  INPUT  DATA  SPECIFIES  A  LENGTH  32  FULLBAND 
DIFFERENTIATOR  WITH  SLOPE  1  AND  WEIGHTING  OF  i/F. 

THE  GRID  DENSITY  WILL  BE  SET  TO  20. 

32,2,1,20 
0,0.5 
1  .  0 
1  .  0 


COMMON  P 12 , AD , DEU , X , Y , GR ID , DES , WT , ALPHA , IEXT , NFCNS , NGR ID 
COMMON  NITER, TOUT 

DIMENSION  IEXT < 34 > ,AD<34) , ALPHA (34) ,X(34) ,Y(34) 

DIMENSION  H ( 34 ) 

DIMENSION  DES < 544 )  ,GRID(544> (UTC544> 

DIMENSION  EDGE <20 ) ,FX(10) , WTX < 1 0 ) , DEV  I  AT ( 1 0  ) 

DOUBLE  PRECISION  PI2.PI 
DOUBLE  PRECISION  AD,D£V,X,Y 
DOUBLE  PRECISION  GEE.D 
INTEGER  BD1 , BD2 , BD3 , BD4 
DATA  BD1,BD2,BD3 .BD4/1HB, 1HA.1HN, 1HD/ 

INPUT=I1MACH< 1 ) 

I0UT=I1MACH<2) 

INPUT=4 
IOUT=i 

WR ITE ( I OUT ,9000) 

'? 0  0  0  FORhAT  <  22  (  /  )  ) 

P 1  =  4 . 0#DATAN  < 1  ODO) 

PI2=2 . 0D  0  0  *P  I 

C 

S  THE  PROGRAM  IS  SET  UP  FOR  A  MAXIMUM  LENGTH  OF  64,  BUT 
C  THIS  UPPER  LIMIT  CAN  BE  CHANGED  BY  REDIMENSIONINC,  THE 
C  ARRAYS  IEXT,  AD,  ALPHA,  X,  Y,  H  TO  BE  NFMAX/2  +  2. 

C  THE  ARRAYS  DES,  GRID,  AND  WT  MUST  DIMENSIONED 
C  lo<  nFmAX/2  *■  2). 

C 

NFMAX=64 
ISO  CONTINUE 
JTYPE=0 
C 

C  PROGRAM  INPUT  SECTION 
C 

READ (INPUT,*)  NFILT , JTYPE , NBANDS , LGR ID 
IFlNFILT.EQ  0>  STOP 
110  FORMAT (415) 

IF (NFILT , LE . NFMAX  OR . NFILT . GE , 3 )  GO  TO  115 

CALL  ERROR 

STOP 
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MS  IF  <  NBANDS  L£  .  0  )  NBANDS=i 

C 

C  GRID  DENSITY  IS  ASSUMED  TO  BE  16  UNLESS  SPECIFIED 
C  OTHERWISE 
C 

IF (LGRID . L£ . 0 )  LGRID=16 
JB=2*NBANDS 

READ( INPUT ,*)  ( EDGE ( J ) , J=1 , JB ) 

120  FORMAT  ( 4F1S  9) 

READC INPUT ,*)  ( FX < J) , J=i , NBANDS ) 

READ (INPUT,*)  < WTX<  J ) ,  J=i , NBANDS) 

Ir(JTYP£  .GT. 0  AND . JTYPE . L£ .  3)  GO  TG  125 

CALL  ERROR 

STOP 

125  NEG=i 

IF( JTYPE . EQ . 1 )  NEG=0 
NODD=NFILT/2 
NGDD=NFILT-2#N0DD 
NFCNS=NFILT/2 

IF<.NQDD.£Q.  1  .AND  NEG.EQ.O)  NFCNS=NFCNS+ 1 

0 

0  SET  UP  THE  DENSE  GRID.  THE  NUMBER  OF  POINTS  IN  THE  GRID 
C  IS  (FILTER  LENGTH  +  1 ) *GR ID  DENSITY/2 

GRID(1)=EDGE(1) 

DELF=LGRID*NFCNS 
DELF=0 ,5/DELF 
IF  <  NEG . EQ . 0 )  GO  TO  135 
IF(EDGE(1)  .LT.DF.LF)  GRID(i)=DELF 
135  CONTINUE 
J  =  1 
L=i 

LBAND* 1 

140  FUP=EDGE(L+i) 

145  TEnP=GR ID ( J ) 

C 

C  CALCULATE  THE  DESIRED  MAGNITUDE  RESPONSE  AND  THE  WEIGHT 
C  FUNCTION  ON  THE  GRID 
C 

DES< J)=EFF (TEMP ,FX ,WTX,LBAND , JTYPE) 

WT( J)=WATE( TEMP ,FX,WTX ,LBAND, JTYPE) 

J=J+1 

GR ID( J ) =TEMP  +  DELF 
IF(GRID(J) .GT.FUP)  GO  TO  ISO 
GO  TO  14S 

ISO  GRID ( J-l ) =FUP 

DES ( J-l ) =£FF ( FUP , FX , UTX , LBAND , JTYPE ) 

WT( J-l )=WATE(FUP , FX,WTX, LBAND , JTYPE) 

LBAND=LBAND+1 

L=L+2 

IF(LBAND .GT . NBANDS)  GO  TO  160 


73 


r 


GRID< J)=EDGE<L) 

GO  TO  140 
160  NGRID=J-1 

IF<NEG .NE .NODD)  GO  TO  165 

IF (GRID(NGRID) .GT . <0 ,5-DELF) >  NGR ID=NGR ID-i 
165  CONTINUE 
C 

C  SET  UP  A  NEW  APPROXIMATION  PROBLEM  WHICH  IS  EQUIVALENT 
C  TO  THE  ORIGINAL  PROBLEM 
C 

IF(NEG)  170,170,180 
170  IF (NODD .  EQ  1 )  GO  TO  200 
DO  175  J=1 , NGR I D 
CHANGE=DCOS<PI*GRID< J)  ) 

DES(J)=DES(J) /CHANGE 
175  WT<J)=WT(  J)*CHANGE 
GO  TO  200 

100  IF ( NODD . EQ . 1 )  GO  TO  190 

DO  185  J=1 , NGR ID 
CHANGE=DSIN<PI*GRID<J> ) 

DES< J)=DES< J) /CHANGE 
185  WT ( J)=UT  < J)*CHANGE 

GO  TO  200 

190  DO  195  J-l  ,NGRID 

CHANGE- DS I N< P I2*GR ID ( J ) ) 

DES(J)=DES(J) /CHANGE 
195  WT  (  J ) =UT ( J  >  ^CHANGE 

C 

C,  INITIAL  GUESS  FOR  THE  EXTREMAL  FREQUENCIES— EQUALLY 
C  SPACED  ALONG  THE  GRID 

C 

200  TEMP=FLOAT(NGRID-l ) /FLOAT (NFCNS) 

DO  210  J=1 , NFCNS 
XT= J-l 

2l0  IEXT< J)=XT#TEMP+1 . 0 
I EXT ( NFCNS+i ) =NGR ID 
NMi=NFCNS-l 
NZ=NFCNS+ 1 
C 

C  CALL  THE  REMEZ  EXCHANGE  ALGORITHM  TO  DO  THE  APPROXIMATION 
C  PROBLEM 
C 

CALL  REMEZ 

C 

C  CALCULATE  THE  IMPULSE  RESPONSE. 

C 

IF(NEG)  300,300,320 
300  IF ( NODD . EQ  .  0 )  GO  TO  310 

DO  305  J=i,NMi 
NZM J=NZ- J 

305  H  (  J )  =0  .  5* ALP HA  <NZMJ ) 


74 


H(NFCNS)=ALPHA(i ) 

GO  TO  350 

310  H(  1  )  =  0 . 25j|tALP  HA  ( NFCNS  ) 

DO  315  J=2 , NMi 
NZM J=NZ-J 
NF2J=NFCNS+2-J 

315  H<J)  =  0 .  25#  (  ALP HA <  NZM J ) + ALP HA  <  NF2J ) ) 

H<  NFCNS ) =0 ,5#ALPHA(l)+0 ,25#ALPHA(2) 

GO  TO  350 

320  T.F(NODD.EQ.Q)  GO  TO  330 
H ( 1 )  =  0 . 25* ALPHA  (NFCNS ) 

H  <  2)  =  0 . 25#ALPHA  <  NMI ) 

DO  325  J=3 , NMi 
NZH J =NZ-J 
NF3J=NFCNS+3-J 

325  H< J)  =  0 . 25#  <  ALPHA  (NZMJ)-ALPHA(NF3J> ) 

H<  NFCNS ) =0 . S#ALPHA ( 1 )  -0 ,25#ALPHA<3) 

H  <  NZ ) =0 . 0 
GO  TO  350 

330  H < 1 ) =0 . 25# ALP HA  <  NFCNS ) 

DO  335  J=2,NM1 
NZMJ=NZ-J 

NF2 J=NFCNS+2- J 

335  H ( J ) =0 . 2S#  (  ALPHA (  NZM J ) -ALPHA ( NF2 J ) ) 

H ( NFCNS )  =0 .5*ALPHA<i)-Q . 25* ALP HA ( 2 ) 

C 

C  PROGRAM  OUTPUT  SECTION. 

C 

350  WRITE (IOUT, 360) 

360  FOPMAT(iQ(/)  ,  7  0  ( 1H#  ) // 1SX  ,  29FINI  TF_  IMPULSE  RESPONSE  (FIR)/ 

113X,34HLINEAR  PHASE  DIGITAL  FILTER  DESIGN/ 

21.7X  ,  24HREMEZ  EXCHANGE  ALGORITHM/) 

IF< JTYPE.EQ. 1 )  URITE(IOUT,365> 

565  FORMAT ( 22X , 15HBANDP ASS  FILTER/) 

IF (JTYPE.EQ. 2)  UR ITE ( IOUT , 370 ) 

370  FORMAT (22X , 14HDIFFERENTI ATOR/ ) 

IF ( JTYPE . EQ . 3 )  WRITE( IOUT, 375) 

375  FORMAT ( 20 X, 19HHILBERT  TRANSFORMER/) 

WRITE(I0UT,378)  NFILT 
373  FORMAT(20X,16HFILTER  LENGTH  =  ,13/) 

WRITE (IOUT, 330) 

300  FORMAT ( 1SX ,28H#*##*  IMPULSE  RESPONSE  **#**) 

DO  301  J=l, NFCNS 
K=NFILT+i-J 

IF(NEG.EQ.O)  WRITE( IOUT, 332)  J,H(J),K 
IF(NEG.EQ.l)  UR ITE ( IOUT , 383  >  J,H(J),K 

331  CONTINUE 

332  FORMAT ( 13X , 2HH ( ,12, AH)  =  ,E15.8,5H  =  H(,I3,1H)) 

383  F0RMAT(13X, 2HH ( , 12 , 4H )  =  ,E15  8,6H  =  -H(,I3,1H)) 

IF(NEG. EQ . 1 . AND .NODD. EQ . 1 )  UR ITE ( IOUT , 384 )  NZ 
334  FORMAT < 13X,2HH( , I2,8H)  *  0.0) 


DO  450  K=i,NBANDS,4 
l(UP=K+3 

IF (KUP  GT.NBANDS)  KUP=NBANDS 

WRITE( IOUT ,385)  ( BDi. , BD2 , BD3 , BD4 , J , J=K , KUP ) 

385  FORMAT  (/24X,4(4Ai , 13 , 7X ) ) 

WRITE  (IOUT,  390  )  ( EDGE  ( 2*  J-i  )  ,  J=l<  ,KUP) 

390  FORMAT (2X,iSHL0UER  BAND  EDGE , 5F14 . 7 ) 

WRITE (IOUT,  395)  ( EDGE ( 2* J ) , J=K , KUP ) 

395  FORMAT (2X, 15HUPPER  BAND  EDGE, 5F 14.7) 

IF ( JTYPE . NE . 2)  WR ITE( IOUT , 400 )  ( FX ( J ) , J=K , KUP ) 

400  FORMAT ( 2X , 13HDESIRED  VALUE , 2X , 5Fi 4 . 7 ) 

IF ( JTYPE . EQ . 2)  WR ITE ( IOUT , 405 )  ( FX ( J ) , J=K , KUP ) 

40S  FORMAT (2X, 13HDESIRED  SLOPE , 2X, 5F i 4 . 7 ) 

WRITE(IOUT,4iO)  ( WTX ( J ) , J=K , KUP ) 

4  i  0  FORMAT ( 2X , 9HWEIGHTING , 6X , SFi 4 .7) 

DO  420  J=K  ,KUP 
420  DEVIAT(J)=OEV/WTX(J) 

WRITE (I OUT ,42S)  (DEVI AT( J ) , J=K , KUP ) 

425  F0RMAT(2X, 9HDEVI ATION ,6X,5Fi4.7) 

IF ( JTYPE , NE. 1)  GO  TO  450 
DO  430  J=K  , KUP 

430  DEV I AT ( J ) -20 . 0*ALOGT ( DEVIAT ( J ) +FX ( J ) ) 

UR ITE (IOUT, 435)  (DEVIAT(J) ,J=K  ,KUP) 

435  FORMAT (2X, 15HDEVI ATION  IN  DB,5Fi4.7) 

450  CONTINUE 

DO  452  J=i , NZ 
IX=IEXT ( J ) 

452  GRID ( J) =GR ID ( IX ) 

UR ITE (IOUT, 455)  ( GR ID ( J ) , J=i , NZ ) 

455  FORMAT ( /2X , 47HEXTREMAL  FREQUENCIES— MAXIMA  OF  THE  ERROR  CURVE/ 

1  <  2X , 5F12 . 7 ) ) 

WRITE( IOUT, 460 ) 

460  FCRMAT(/iX,70(iH*),10(/)  ) 

C  GO  TO  100 

STOP 
END 
C 

C - 

C  FUNCTION-.  EFF 

C  FUNCTION  TO  CALCULATE  THE  DESIRED  MAGNITUDE  RESPONSE 
C  AS  A  FUNCTION  OF  FREQUENCY. 

C  AN  ARBITRARY  FUNCTION  OF  FREQUENCY  CAN  BE 

C  APPROXIMATED  IF  THE  USER  REPLACES  THIS  FUNCTION 

C  WITH  THE  APPROPRIATE  CODE  TO  EVALUATE  THE  IDEAL 

C  MAGNITUDE.  NOTE  THAT  THE  PARAMETER  FREQ  IS  THE 

C  VALUE  OF  NORMALIZED  FREQUENCY  NEEDED  FOR  EVALUATION 

C - 

C 

FUNCTION  EFF (FREQ ,FX , WTX ,LBAND , JTYPE ) 

DIMENSION  FX(5) ,UTX(5) 

IF (JTYPE . EQ . 2)  GO  TO  1 
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EFF=FX  <  LBAND ) 
RETURN 

1  EFE=FX ( LBAND ) *FREQ 

RETURN 
END 


FUNCTION  =  WATS 

FUNCTION  TO  CALCULATE  THE  WEIGHT  FUNCTION  AS  A  FUNCTION 
OF  FREQUENCY.  SIMILAR  TO  THE  FUNCTION  EFF,  THIS  FUNCTION  CAN 
BE  REPLACED  BY  A  USER-WRITTEN  ROUTINE  TO  CALCULATE  ANY 
DESIRED  WEIGHTING  FUNCTION. 


FUNCTION  WATE<  FREQ , FX , WTX ,  LB AND , JTYPE) 
DIMENSION  FX(5),WTX(S) 

IF (JTYPE . EQ . 2 )  GO  TO  1 
WATE=WTX< LBAND) 

RETURN 

IF(FX(LBAND) .LT. 0 . 0001)  GO  TO  2 
WATE=WTX  < LBAND )  /FR  EG) 

RETURN 

WATE=UTX< LBAND) 

RETURN 

ENI) 


SUBROUTINE  ERROR 

THIS  ROUTINE  WRITES  AN  ERROR  MESSAGE  IF  AN 
ERROR  HAS  BEEN  DETECTED  IN  THE  DATA. 


SUBROUTINE  ERROR 
COMMON  IOUT 
WRITE< IOUT  ,  i ) 

i  FORMAT ( 44H  ************  ERROR  IN  INPUT  DATA  **********> 

RETURN 
END 
C 

C - 

C  SUBROUTINE:  REMEZ 

C  THIS  SUBROUTINE  IMPLEMENTS  THE  REMEZ  EXCHANGE  ALGORITHM 
C  FOR  THE  WEIGHTED  CHEBYSHEV  APPROXIMATION  OF  A  CONTINUOUS 
C  FUNCTION  WITH  A  SUM  OF  COSINES  INPUTS  TO  THE  SUBROUTINE 
C  ARE  A  DENSE  GRID  WHICH  REPLACES  THE  FREQUENCY  AXIS,  THE 
C  DESIRED  FUNCTION  ON  THIS  GRID,  THE  WEIGHT  FUNCTION  ON  THE 
C  GRID,  THE  NUMBER  OF  COSINES,  AND  AN  INITIAL  GUESS  OF  THE 
C  EXTREMAL  FREQUENCIES.  THE  PROGRAM  MINIMIZES  THE  CHEBYSHEV 
C  ERROR  BY  DETERMINING  THE  BEST  LOCATION  OF  THE  EXTREMAL 
C  FREQUENCIES  (POINTS  OF  MAXIMUM  ERROR)  AND  THEN  CALCULATES 
C  THE  COEFFICIENTS  OF  THE  BEST  APPROXIMATION 

C - 
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non 


COMMON  NITER, IOUT 

DIMENSION  IEXT ( 34 ) ,AD<34) ,ALPHA<34) ,X<34> ,Y(34) 
DIMENSION  DES < 544 ) ,GRID<S44> ,WT<544> 

DIMENSION  A<34) ,P<65) ,Q<6S> 

DOUBLE  PRECISION  PI2,DNUM, DDEN , DTEMP , A,P , Q 
DOUBLE  PRECISION  DK , DAK 
DOUBLE  PRECISION  AD,DEV,X,Y 
DOUBLE  PRECISION  GEE  ,  D 

THE  PROGRAM  ALLOWS  A  MAXIMUM  NUMBER  OF  ITERATIONS  OF  25 

ITRMAX=25 
DEVL=-i . 0 
NZ=NFCNS+1 
NZZ=NFCNS+2 
NITER=0 
100  CONTINUE 

IEXT<NZZ)=NGRID+i 

NITER=NITER+1 

IF<NITER .GT . ITRMAX)  GO  TO  400 
DO  110  J=t,NZ 
JXT  =  IEXT  ( J) 

DTEMP  =  GR ID  <  JXT ) 

DTEMP=DC0S<DTEMP*PI2> 

110  X(J)=DTEMP 

JET=<NFCNS-t)/15+l 
DO  120  J=i ,NZ 
120  AD< J)=D< J,NZ, JET) 

DNUM=0  0 
DDEN=0 . 0 
K  =  1 

I> '  130  J  =  1  , NZ 
L=IEXT(J) 

DTEMP=AD  <  J ) #DE3  <L ) 

DNUM=DNUM+DTEMP 

DTEMP =FLQ AT  <K ) #AD ( J ) /WT ( L) 

DDEN=DDEN+DTEMP 

130  K=-K 
DEV=DNUM/DDEN 
WRITE < IOUT , 131 )  DEV 

131  FORMAT  < IX , 12HDEVIATI0N  =  ,F12.?> 

NIJ=i 

IFCDEV  GT  .  0  .0)  NU=-1 
DEV  =  -FLQAT  <NU)*DEV 
K=NU 

DO  140  J=i , NZ 
L=IEXT< J) 

DTEMP=FLOAT( K ) *DEV/UT  <L> 
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Y(J)=DES(L>+DTEMP 
140  K  =  -K 

IF ( DEV . GT . DEUL )  GO  TO  150 
CALL  OUCH 
GO  TO  400 
ISO  DEUL=DEV 
JCHNGE=0 
K 1= IEXT ( 1 ) 

KN2= IEXT  <  N2 ) 

KLOU=0 
NUT=-NU 
•  J  =  1 
C 

C  SEARCH  FOR  THE  EXTREMAL  FREQUENCIES  OF  THE  BEST 
C  APPROXIMATION 
C 

200  IF ( J . EQ . N22  >  YN2=C0MP 
IF(J.GE.NZZ)  GO  TO  300 
KUP=IEXT( J+i> 

L=IEXT ( J ) +1 
NUT  =  -NUT 

IF ( J  .  EQ . 2 )  Yi=COMP 
COMP=DEU 

IF(L.GE.KUP)  GO  TO  220 
ERR=GEE  <  L , N2  > 

ERR  =  <  ERR-DES ( L ) ) #UT<  L  > 

DTEMP=FLOAT< NUT >*ERR -COMP 
IF ( DTEMP . L£ . 0 . 0 )  GO  TO  220 
COMP=FLOAT <  NUT) ¥ERR 
210  L=L+1 

IF < L .GE.KUP)  GO  TO  215 
ERR=GEE(L,NZ) 

ERR=(ERR-DES(L) )*WT(L> 

DTEMP=FI.OAT  (  NUT  >  #ER  R  -COMP 
IF < DTEMP . LE . 0 . 0 >  GO  TO  215 
CQMP  =  FLQAT <  NUT ) KERR 
GO  TO  210 
215  IEXT  < J)=L-i 

J  =  J+1 
KLOU=L-i 
JCHNGE=JCHNGE+1 
GO  TO  200 
220  L=L- 1 
225  L=L-i 

IF(L.LE.KLOU)  GO  TO  250 
ERR=GEE(L,N2) 

ERR=(ERR-DES(L) >KUiT<L> 

DTEMP=FLOAT< NUT) KERR -COMP 
IF (DTEMP  GT.0.0)  GO  TO  230 
IF ( JCHNGE . LE  0)  GO  TO  225 
GO  TO  260 
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GO  TO  200 
250  L=IEXT ( J ) +1 

IF ( JCHNGE . GT . 0 )  GO  TO  215 
255  L=L+t 

IFCL.GE .KUP)  GO  TO  260 
ERR=GEE<L,NZ) 

ERR  =  ( ERR-DES  <L  ) )  #WT<L. ) 
DTEMP=FLOAT ( NUT ) #ERR-COMP 
IF (DTEMP . LE . 0 . 0 )  GO  TO  255 
COMP =FLOAT( NUT )*ERR 
GO  TO  210 

26 1)  K  LOU  =  IEXT  <  J ) 

J  =  J  +  1 
GO  TO  200 

300  IF(J.GT.NZZ)  GO  TO  320 

IF (K 1 . GT . IEXT ( 1) )  K i  =  IEXT ( 1 ) 
IF(KNZ.LT . IEXT(NZ) )  KNZ=IEXTCNZ) 
NUTi=NUT 
NUT=-NU 
L  —  0 

KUP  =K 1 

COMP=YNZ*( 1.00001) 

LUC!<  =  1 
310  L=L+1 

IFCL.GE.KUP)  GO  TO  315 
ERR=GEE ( L , NZ ) 

ERR= ( ERR-DES ( L) ) #UT(L  ) 

DTEMP  =FLOAT ( NUT ) #ER  R -COMP 
IF ( DTEMP . LE  .  0 . 0 )  GO  TO  310 
COMP=FLOAT(NUT)#ERR 
J=NZZ 
GO  TO  210 
315  LUCK=6 

GO  TO  325 

320  IFCLUCK .GT.?)  GO  TO  350 
IF ( COMP  .  GT . Y1 )  Y l=COMP 
K1 = IEXT ( NZZ ) 

325  L=NGR ID+1 
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KLOW=KNZ 
NUT =-NUT 1 
COMP=Yi#(i . 00001) 

330  L=L-i 

IF  ( L  . LE . KLOW )  GO  TO  340 
ERR=GE£(L,NZ> 

£RR=(ERR-DES(L) )#UT (L) 

DTEMP=FLOAT  <NUT)*ERR-COMP 
IF  (DTEMP , LE . 0 . 0 )  GO  TO  330 
J=NZZ 

COMP=FLOAT (NUT ) *ERR 
LUCK=LUC!<  +  i0 
GO  TO  235 

340  IF (LUCK .EQ .6)  GO  TO  370 
DO  345  J=1 ,NFCNS 
NZZMJ =nzz-j 
NZMJ=NZ-J 

34S  IEXT (NZZMJ)=IEXT (NZHJ) 

IEXT ( 1 ) =K  1 
GO  TO  100 
350  KN=IEXT ( NZZ) 

DO  360  J=i ,NFCNS 
360  IEXT( J ) = IEXT ( J+i ) 

IEXT ( NZ ) =KN 
GO  TO  100 

70  IF(JCHNGE.GT.O)  GO  TO  100 

C  CALCULATION  OF  THE  COEFFICIENTS  OF  THE  BEST  APPROXIMATION 
C  USING  THE  INVERSE  DISCRETE  FOURIER  TRANSFORM 
C 

400  CONTINUE 

NMi=NFCNS-i 
FSH=1 . 0E-06 
GTEMP=GR ID ( 1 ) 

X(NZZ)=-2. 0 
CN=2*NFCNS-i 
DELF=1 . 0/CN 
L=i 
KKK  =  0 

JF(GRI.D(i  )  .LT  .  0  .  01  .  AND  .GRID(NGRID)  .  GT  .  0  49)  KKK  =  1 
IF(NFCNS .LE . 3)  KKK=i 
IF(KKK  .EQ. 1)  GO  TO  405 
DTEMP=DC0S(PI2*GRID(  1. )  ) 

DNUM=DCOS(PI2*GRID(NGRID) ) 

AA=2  0/ (DTEMP-DNUM) 

BB=- (  DTEMP-t-DNUM )  /  (  DTEMP -DNUrt  ) 

405  CONTINUE 

DO  430  J=i ,NFCNS 
FT=FLOAT( J-i) 

FT=FT*DELF 

XT=DC0S(PI2*FT) 
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IF  (KKK .EQ. i)  GO  TO  410 
XT=<  X  r-BB ) / AA 
XT1=SQRT<1 .0-XT*XT) 
FT=ATAN2<XTi,XT)/PI2 
410  XE=X<L) 

IF(XT.GT.XE)  GO  TO  420 
IF((XE-XT) .LT.FSH)  GO  TO  415 
L=L+i 
GO  TO  410 
415  A(J)=Y(L> 

GO  TO  42S 

420  IF< (XT-XE) .LT .FSH)  GO  TO  41S 
GR ID  < 1 ) =FT 
A(J)=GEE(i ,NZ) 

42S  CONTINUE 

IF(L.GT.i)  L=L-1 
430  CONTINUE 

GR ID (  1 )=GTEMP 

DDEN=P I2/CN 

DO  S10  J= 1 , NFCNS 

DTEMP=0 . 0 

DNUM=J-1 

DNUH=DNUM*DDEN 

IF (NM1 .LT. 1)  GO  TO  505 

DO  500  K  =  1 , NM1 

DAK=A(K+i) 

DK  —  K 

500  DTEHP=DTEMP+DAK*DCOS<DNUM*DK> 

SOS  DTENP=2.0*DT£MP+A<1) 

510  ALPHA  <J)=.DTEffP 

DO  550  J=2, NFCNS 
SS0  ALPHA< J)=2.0*ALPHA(J)/CN 
ALPHA < 1)=ALPHA< 1 )/CN 
I F  <  K  K  K  EQ  .  1 )  GO  TO  545 
P  < 1 ) =2 . 0* ALP HA  <  NFCNS ) ALPHA ( NHi ) 
P (2) =2. 0*AA*ALPHA< NFCNS) 

Q< i)=ALPHA< NFCN5-2)- ALP  HA (NFCNS) 

DO  540  J=2  ,  Ntl  1 
IF(J.LT.NMi)  GO  TO  515 
AA=0 ,S*AA 
BB=0 .S*BB 
515  CONTINUE 

P(J+1)=0 . 0 
DO  520  K=1 , J 
A  ( K ) =P ( K  ) 

S20  PCO=2.0*BB*A(K) 

P(2)=*P(2)+A(1)#2.0*AA 

JH1=J-1 

DO  S25  K=i , JM1 

525  P(K)=P(K)+Q(K)+AA*A(K+1> 

JP 1* J+l 
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DO  530  K=3,JP1 
530  P(K)=P(K)+AA*A(K-i) 

IF ( J . EQ . NM1 )  GO  TO  540 
DO  535  K=i ,  J 
535  QUO  =  -A ( K  ) 

NFiJ=NFCNS-l-J 

Q(i >=Q(1)+ALPHA(NF1J) 

540  CONTINUE 

DO  543  J=l, NFCNS 
543  ALPHA(J)=P ( J) 

545  CONTINUE 

IF (NFCNS  . GT . 3)  RETURN 
ALPHA(NFCNS+1 >  =  0  .  0 
ALPHA (NFCNS+2)=Q . 0 
.RETURN 
END 

n 

w 

c - 

C  FUNCTION i  D 

C  FUNCTION  TO  CALCULATE  THE  LAGRANGE  INTERPOLATION 
C  COEFFICIENTS  FOR  USE  IN  THE  FUNCTION  GEE. 

C - 

C 

DOUBLE  PRECISION  FUNCTION  D(K,N,M> 

COMMON  P 1 2 , AD , DEV , X , Y , GR I D , DES , UT , ALP  HA , I EXT , NFCNS , NCR  I D 
DIMENSION  IEXT (34 > ,AD(34) ,ALPHA(34> ,X(34> ,Y<34> 

DIMENSION  DES < 544 J , GRID (544) ,WT(S44> 

DOUBLE  PRECISION  AD,D£V,X,Y 

DOUBLE  PRECISION  Q 

DOUBLE  PRECISION  PI2 

D=1  .  0 

a=x<K  > 

DO  3  L=i,M 
DO  2  J=L,N,M 
IF  ( J-K) 1,2,1 

1  D=2. 0*D*(Q-X(J) ) 

2  CONTINUE 

3  CONTINUE 
D=1 . 0/D 
RETURN 
END 

C 

C - 

C  FUNCTION > GEE 

C  FUNCTION  TO  EVALUATE  THE  FREQUENCY  RESPONSE  USING  THE 
C  LAGRANGE  INTERPOLATION  FORMULA  IN  THE  BARYCENTRIC  FORM 


C 

DOUBLE  PRECISION  FUNCTION  GEE (K  ,N> 

COMMON  P 12, AD, DEV, X,Y, GRID, DES, UT, ALPHA, TEXT, NFCNS, NGR ID 
DIMENSION  IEXT (34) ,AD(34> , ALPHA ( 34 > , X ( 34 > , Y ( 34  ) 


o  n  n  n  n  o  n  r 


DIMENSION  DES<  544 ) ,GRID<544> ,UTC544) 
DOUBLE  PRECISION  P,C,D,XF 
DOUBLE  PRECISION  PI2 
DOUBLE  PRECISION  AD,DEV,X,Y 
P=0  0 

XF=GRID<K  ) 

XF=DC0S<PI2*XF) 

D=0  .  0 

DO  1  J=*i,N 
C=XF-X( J) 

C=AD< J)/C 
D=D+C 

1  P=P+C*Y<J) 

GEE=P/D 

RETURN 

END 

C 


SUBROUTINE.  QIJCH 

WRITES  AN  ERROR  MESSAGE  WHEN  THE  ALGORITHM  FAILS  TO 
CONVERGE .  THERE  SEEM  TO  BE  TWO  CONDITIONS  UNDER  WHICH 
THE  ALGORITHM  FAILS  TO  CONVERGE.  (1)  THE  INITIAL 
GUESS  FOR  THE  EXTREMAL  FREQUENCIES  IS  SO  POOR  THAT 
THE  EXCHANGE  ITERATION  CANNOT  GET  STARTED,  OR 
(2)  NEAR  THE  TERMINATION  OF  A  CORRECT  DESIGN, 

C  THE  DEVIATION  DECREASES  DUE  TO  ROUNDING  ERRORS 

C  AND  THE  PROGRAM  STOPS.  IN  THIS  LATTER  CASE  THE 

C  FILTER  DESIGN  IS  PROBABLE  ACCEPTABLE,  BUT  SHOULD 

C  BE  CHECKED  BY  COMPUTING  A  FREQUENCY  RESPONSE. 

C - 

C 

SUBROUTINE  OUCH 
COMMON  NITER, IOUT 
WRITE (IUUT,  UNITER 

l  FORMAT  ( 44H  *##*#*##*#*#  FAILURE  TO  CONVERGE  * *»*:** ****  ,  3 C  /  )  , 

14QHPR0BABLE  CAUSE  IS  MACHINE  ROUNDING  ERROR, 3</>, 

222HNUMBER  OF  ITERATIONS  «,I4,3</>, 

330HIF  THE  NUMBER  OF  ITERATIONS  EXCEEDS  3,3(/>, 

4t.lH  THE  DESIGN  MAY  BE  CORRECT,  BUT  SHOULD  BE  VERIFIED  WITH  AN  FFT ) 
RETURN 
END 
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FIRF2  Output 


DEVI AT  ION  = 
DEVIATION  = 
DEVIATION  = 
DEVIATION  = 
DEVIATION  = 
DEVIATION  = 
DEVIATION  = 
DEVIATION  = 


.000734616 
.006314959 
. 021566622 
.026202958 
- . 032680455 
- .034435450 
-  034448376 
- . 034448590 


BANDPASS  FILTER 
FILTER  LENGTH  =  55 


*****  IMPULSE  RESPONSE  ***** 


H<  1) 

= 

. 1 0662768E-02 

= 

H  ( 

55) 

H  (  2) 

= 

. 63777557E-02 

= 

H  ( 

54) 

H  (  3) 

= 

. 35755956E-02 

= 

H  ( 

53) 

H  (  4) 

= 

- . 90673427E-02 

H( 

52) 

H  (  S) 

= 

-  .  90907477E-Q2 

= 

H  ( 

51) 

H  (  6) 

= 

.29156702E-02 

= 

H  ( 

50) 

H<  7) 

= 

. 39638160E-02 

= 

H  ( 

49) 

H  (  8) 

= 

. 1 1 171967E-0 1 

= 

H  ( 

48) 

H (  9) 

= 

. 11646770E-01 

= 

H  ( 

47) 

H<  10) 

= 

- . 99630672E-02 

= 

H( 

46) 

H  <  1 1 ) 

= 

- . 92384685E-02 

= 

H< 

45) 

H(12> 

= 

-.20406403E-01 

= 

H( 

44  ) 

H(  13) 

= 

-  19460432E-0 1 

= 

H( 

43) 

H<  14) 

= 

. 31243  082E- 0 1 

= 

H( 

42) 

H  ( 15 ) 

= 

.63045118E-02 

= 

H  ( 

41) 

H  ( 16 ) 

= 

- .20482879E-01 

= 

H  ( 

40) 

Hv  17) 

= 

.  65741 474E- 02 

= 

H  ( 

39) 

H  <  13  > 

= 

-  11201 589E- 02 

= 

H< 

38) 

H<  19) 

= 

4 1956394E-0 1 

= 

H  ( 

37) 

H  <  20  ) 

= 

. 35784297E- 0 1 

— 

H  < 

36) 

H(2i  > 

= 

. 34744851E-0 1 

= 

H  ( 

35) 

H(22) 

= 

. 7 1 496233E-0 1 

= 

H  ( 

34) 

H  <  23 ) 

= 

-. 17 138839E+00 

= 

H  ( 

33) 

H  ( 24 ) 

= 

- . 18255037E+00 

H< 

32) 

H  ( 25  > 

= 

. 74059084E-0 1 

= 

H  ( 

31  ) 

H  ( 26 ) 

= 

- . 103174 25 E+00 

= 

H  ( 

30) 

H  ( 27 ) 

= 

. 25716718E-0 1 

= 

H< 

29) 

H  ( 28 ) 

= 

378 13556E  +  0  0 

S 

H  ( 

28) 
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BAND  1 

BAND  2 

BAND  3 

BAND  4 

LOWER  BAND  EDGE 

.0000000 

.1000000 

1800000 

3000000 

UPPER  BAND  EDGE 

. 0500000 

. 1500000 

2500000 

.3600000 

DESIRED  VALUE 

0000000 

1.0000000 

0000000 

1.0000000 

WEIGHTING 

10  .  00000Q0 

1.0000000  3 

0000000 

1.0000000 

DEVIATION 

. 0034449 

. 0344486 

0114829 

0344486 

DEVIATION  IN  DB 

-49 . 2565765 

.2941777  -38 

7990036 

2941777 

BAND  5 

LOWER  BAND  EDGE 

. 4100000 

UPPER  BAND  EDGE 

.5000000 

DESIRED  VALUE 

0000000 

WEIGHTING 

20.0000000 

DEVIATION 

. 0017224 

DEVIATION  IN  DB 

-55 . 2771693 

EXTREMAL  FREQUENC 

IES — MAXIMA 

OF  THE  ERROR  CURVE 

.0000000 

0167411 

0323661  0446429 

0500000 

.1000000 

1089285 

1267856  1424105 

. 1500000 

.1900000 

1855803 

1978571  .2134820 

. 23C2230 

.2436153 

2500000 

3000000  .3122771 

. 3323668 

.3502244 

3600000 

4100000  .4155305 

4289737 

.4457151 

4635727 

4814302  .5000000 

********************************************************************** 


AD-A118  066  AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON  AFB  OH  SCHOO— ETC  F/G  17/2 

A  GENERAL  PURPOSE  MINI-COMPUTER  BASED  DIGITAL  signal  PROCESSING— ETC (U) 
JUN  82  M  TODD 

UNCLASSIFIED  AFIT/GE/EE/82J-14  NL 


CTM4«  L 


PROGRAM  covau 


c  MAIN  PROGRAM:  TEST  PROGRAM  POP  COVAP  AMD  AUTO 
C  AUTHORS:  A.  H.  GRAY?  JR.  AMD  J.  D.  MARKEL 

S  GRAY  -  UN I V .  OF  CALIF..  SANTA  BARBARA.  CA  93109 

S  EOTH  -  SIGNAL  TECHNOLOGY.  INC..  15  M.  DE  LA  GUERRA 

C  STREET.  SANTA  ERABRPA.  CA  93101 


D I MENS  I ON  X < 100>  *•  A<S1 > .  RC <2 1 > 

IOUTD  =  1 
X<1>  =  1 . 

X  <2>  =  1.44 

X  (  3  1  =  1 . 44* (v  )  —  1 .  it. 

DO  10  J=4 .40 

X  t  J>  *  1 . 44*X  ( J— 1>  -  1 . 2$*X  <J-2>  +  ..51*XCJ-'S> 
10  CONTINUE 
9999  CCRMAT  I  S .  i'E  1 5 .  €  > 

339 3  FORMAT '.E12. 

N  =  40 
DO  20  M=2.  3 

CALL  COVAR  <N .  X .  M  >  A .  ALPHA .  R:C  > 

MR  =  M+l 

i.i i R  ITS  <1 OUTD .  9939''  < I .  A  <I>  .  PC  < I >  .  1=1.  ‘IP ) 

MR  I  TE  <  I OUTD . 9992 "■  ALPHA 
CALL  AUTO  CM.  X.  M.  A.  ALPHA.  RC"' 

MR  I  TE  < I  OUTD.  9999  ■  <1  •  A  <  I  >  .  RC  >  I  •  .  I  —  1 .  ■ 

MR  I TE  < I OUTD . 9992  '  ALPHA 
20  CONTINUE 

DO  30  JE=1 . 40 
J  =  44— JB 
j>  =  x<  J-3> 

20  CONTINUE 

DO  4  0  J=1 >  3 
JE  =  43+ J 
X<J.J  =  0. 

X <  JB>  =  0. 

40  CONTINUE 

M  =  3 
N=4* 

CALL  COVAR 1 N . X . M . A . ALPHA . Rc ' 

MR  =  M+l 

MR ITE-'I OU TD .  9999>  < I .  A  (I  1  .  PC  \ I  ">  .  I  =  1 .  MR'* 

MR  I  TE  (  I  OUTD .  9992  )  ALPHA 
CALL  AUTO  <N.  X.  M»  A.  ALPHA .  RC.* 

MR  I  TE  -T  IOUTD.  9999)  <  I  <  A  ( I  >  .  PC  <1 )  .  1  =  1 .  f,?'* 

MR I TE (IOUTD. 9993)  ALPHA 

:top 

END 
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INPUTS:  N  -  NO.  DP  DATA  POINTS 

X  CN.>  -  INPUT  DATA  SEQUENCE 
M  -  ORDER  OF  FILTER  CM <21. SEE  NOTEV> 
OUTPUT::  A  -  FILTER  COEFFICIENTS 

AL°HA  -*  RESIDUAL  "ENERGY". 

GRC  -  "GENERALISED  REFLECTION  COEFFICIENT' 


♦PROGRAM  LIMITED  TO  M*20  BECAUSE  OF  THE  DIMENSIONS 
B  CM*  CM- 1  '•  •'  2>  »  EETA  CM )  >  AND  CC  CM*  1  > 


0 1  MENS  I  ON  X  <  1  .  A  C 1  ,  GRC  <  1  > 

D I  MENS  I  ON  E  •:  1 9  ®>  «  EE  TA  < 2  0)  » CC  C2  1  '• 

MP  =  M+l 
MT  =  CMP*M>  1 
DO  10  J=1 • MT 
B  CJ>  =  0. 

coral  hue 

ALc'HA  -  0. 

:  c  •  i  >  =  o . 

:c 'S'  =  o. 

DO  20  NPaMP.N 
•Nc  l  *  NP-1 

ALPHA  «  ALPHA+X  •:'NP'>  ♦XCNP) 

C C  <1»  =  CC  <  1  >  +  .  <NP.»  ♦X  CMP  1  •• 

C C  C£>  =  CC  <2>  *X  CMP  1  >  *X 1  NP 1  > 
corai  HUE 

bci>  *  i. 

BETA  Cl)  =  CCC2) 

GcC''l  '1  =  —  CC  1  1  >  CC  >'2) 

A  C 1  >  =  1 . 

A  (2)  =  GRC Cl) 

ALPHA  =  ALPHA+GRC  1  1  '•  ♦CC  C 1  > 

MF  =  M 

DO  ISO  M INC-2* MF 
DO  SO  J=1 «MINC 
JP  =  MINC+2-J 
ra  =  mp+i-jp 
N2  =  N+l-MINC 
NS  =  N+2-JP 
N4  =  MP-MINC 

CC  C  JP)  =  CC  <  J? - 1  ■  +X  CN4  '  ♦X  CM  1  "•  -X  <N2>  *X  CN3> 
coraiNUE 
CC  <■  1 )  =  0. 

DO  40  NP*MP» N 
HI  =  NF-MINC 


40  CONTINUE 

MS'JB  =  CM  I NC ♦M I NC  — M  I NC  >  £ 

NN1  =  MINC-1 
MM1  =  MINC-1 
N1  »  MSUB+MINC 
B  cm  i  >  =  i . 

DO  30  IP=1.MM1 

I CUB  *  CIP-*IP-IP'’ '£ 

IF  <BETA  <IP>  )  150,150,50 
50  GAM  =  0. 

DO  60  J=1 j IP 
N1  =  ISUB+J 

CAM  =  G AM+CC  < J+ 1 > ♦ B  < N 1 • 

40  CONTINUE  v 

SAM  =  SAM/' BETA  CIP) 

DO  70  JP=1. IP 
N 1  =  MSUB+JP 
M2  =  ISUB+JP 
p  <  N 1  ;*  «  B  CN1)  -GAM»B  -M2) 

70  CONTINUE 

30  CONTINUE 

BETA '-MI NO  =  0. 

DO  30  J=1  ?  M INC 
N1  -  MSUB+J 

BETA  CM  I  NO  *  BETA  ''MINC)  +CC  1  J+ 1  >  *B  CM1  > 

90  CONTINUE 

IF  i.;BETA  'NINO  15  0,  150.  100 
100  3=0. 

DO  110  IP=1 . MINC 

:  =  c+ccap>»A' ip> 

110  CONTINUE 

3P  C  -M I  NC )  =  -S  BETA  •'MINC'1 
00  120  !P=2.MINC 
M2  =  MSUB+IP-i 

9  <Tjp>  =  A  1  IP:-  +-3PC  ’  MINC  >  ♦  E  CM2' 

120  CONTINUE 

A  CNINC+l >  =  SRC  CM  INC) 

C  =  SRC  CM  INC)  ♦SRC  CHINO  ♦BETA  -..MINC) 

ALPHA  =  ALPHA- S 
IF  '.ALPHA)  150.  150.  130 
130  CONTINUE 
140  RETURN 
150  CONTINUE 

•I  WARNING  -  SINGULAR  MATRIX 

I0UTD  =  1 
WR I TE<  I OLITD. 

9999  FORMAT C34H  WARNING  -  SINGULAR  MATRIX  -  COVAR) 
GO  TO  140 
END 


S  SUBROUTINE:  AUTO 
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in  if  i  r  - 


•:  A  SUBROUTINE  FOP  IMPLEMENTING  THE  AUTOCORRELATION 
:  METHOD  OF  LI  HEAP  PREDICTION  ANALYSIS 


SUBROUTINE  AUTO  '  H « X,  M.  A?  ALPHA.  PO 

:  INPUTS:  N  -  NO.  OF  DATA  POINTS 

C  XCN>  -  INPUT  DATA  SEQUENCE 

C  M  -  ORDER  OF  FILTER  <M<£1.  SEE  NOTE*-) 

C  OUTPUTS:  A  -  CILTEP  COEFFICIENTS 

C  ALPHA  -  RESIDUAL  "ENERGY 

C  PC  -  REFLECTION  COEFFICIENTS 

C 

■S  ♦PROGRAM  LIMITED  TO  M<£1  BECAUSE  OF  DIMENSIONS  OF  PC.) 

D I MENS  I ON  X  < 1 >  .  ft  C  i  > .  PC < 1 > 

DIMENSION  PC£1.) 

C 

rip  =  m+i 
DO  S  0  )■>  1  *  r-,c 

p  i-;iO  =  0. 

NK  =  N-K+l 

DO  10  np=i  .m 
N 1  =  NP+^-1 

R  NO  =  P  CK>  +X  CNP*  *X  <  N  1  > 

10  CONTINUE 

SO  CONTINUE 

S!-C<1>  =  -P  02>  •  R  <  1 ') 

ftci>  =  l. 

ACS)  =  PC'-'1> 

ALPHA  a  R  1  :•  +P  ■  £)  ♦PC  <  1 .) 

DO  50  M INC= I ? M 
S=0. 

DO  50  IP-=l,MINC 
N 1  =  MINC-IP+i 
S  -  S+R  ''Ml '  ♦A  '.IP"1 
50  CONTINUE 

P'S  CM  I  NS.)  =  -S-  'ALPHA 
MH  a  M INC '2+1 
DO  4  0  I  ='=■£■  »MH 
IB  =  MINC-IP+E 
AT  =  A  < I p) +RC CM I NC > »A C I E  > 
ft  ■'  I  £  •  =  ft  C  IB>  +RC  CM  I  NO  +ft  <  I P'1 
A<  IP)  =  AT 
40  CONTINUE 

AcriINC+1  >  =  RC  CM  I  NO 
ALPHA  =  ALPHA+PC CM INC) ♦S 
I F  C A L P H A )  70*70.50 
0  CONTINUE 

0  RETURN 

0  CONTINUE 

WARNING  -  SINGULAR  MATRIX 

IOUTD  =  1 
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l.F  I TE  a  aiJTD 
POPART  -::33H  WAR  IN  IMG 
60  TO  60 


1 

4. 

. 1 00000E+01 

-. 391 1 17E+00 

D 

w 

1 0S020E+01 

. 1 977 1 4E+0 0 

“S 

_« 

. 197714E+00 

. GOOOOOE+OO 

.  172 

093E+01 

1 

.  1 OOOOOE+Ol 

-.  905549E+00 

■“> 

C 

-. 1 18200E+01 

. 233193E+00 

•“V 

It 

. 233198E+00 

. 00000 OE+OO 

.  303 

339E+01 

1 

.  1 00000E+01 

-. 9354S7E+00 

L_ 

144000E+01 

. 313708E+00 

T; 

. 125999E+01 

-. 309995E+00 

4 

309995E+00 

. 00 000 OE+OO 

.  132 

125E-04 

1 

. 1 OOOOOE+Ol 

-. 905549E+00 

£. 

133303E+G1 

. 233193E+00 

. 1 1 0323E+01 

-.71 0033E+00 

4 

—  .710  Q38E+ 0 0 

. OOOOOOE+OO 

-  152 

9 1 3E+  0 1 

1 

.  1 OOOOOE+Ol 

-. 905549E+00 

— . 1 333  03E+0 1 

. 2S3193E+00 

■-» 

. 1 1 0S23E+01 

-. 710039E+00 

4 

-.71 0039E+00 

. OOOOOOE+OO 

.  152 

913E+01 

1 

4 

.  1 OOOOOE+Ol 

-. 905549E+00 

•2> 

-. 133308E+01 

. 23  3 1 93E+  0 0 

.  110SS3E+01 

-. 71 0023E+00 

4 

-.71 0038E+ 00 

. OOOOOOE+OO 

.  152913E+01 
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•-L 


FTN4. L 


PROGRAM  ?  FWI  IR  NAMED  COMMON 


C  MAIN  PROGRAM:  FIJI  IP 

C  AUTHORS:  KENNETH  STEIGLITZ  AND  BRUCE  D.  LADENDOPF 

C  PRINCETON  UNIVERSITY.  PRINCETON.  NJ  08540 

C  VERSION  OCTOBER  15.  1978 

C  INPUT:  A  GRID  OF  FREQUENCY  POINTS  AND  A  SET  OF 

C  '"ESSENTIALLY'1  INFINITE  PRECISION  COEc- 

C  “tCIENTS.  SEE  SUBROUTINE  I  HD AT  FOR 

S  DETAILED  INPUT  INFORMATION. 

C 

C  THIS  PROGRAM  DESIGNS  FINITE  WORD-LENGTH  IIP  DIGITAL 
•:  FILTERS.  the  method  uses  randomized  pattern  search 
C  OF  HOOKE  AND  JEEVES.  AND  IS  DESCRIBED  IN  "DESIGNING 
C  SHORT-WORD  RECURSIVE  DIGITAL  FILTERS?"  BY  KENNETH 
C  STEIGLITZ.  IN  c'ROC.  9TH  ANNUAL  ALLEPTON  CONF.  ON  CIRCUIT 
S  AND  SYSTEM  TH. .  PP. 778-788  OCT.  1971.  REPRINTED  IN 
C  DIGITAL  SIGNAL  PROCESSING  II?  EDITED  EY  THE  DIGITAL 
S  SIGNAL  PROCESSING  COMMITTEE  OF  THE  GROUP  ON  ASSR.  IEEE 
C  PRESS?  M.  V.?  1978. 


DOUBLE  Pc  EC  I S I  ON  ljf.90 » ?  Y-'SO'1  « WEIGHT  'SO-1  ?  X  >'36 1 
DOUBLE  PRECI  SION  SDELTA? DELTAS . DUMMY «  RHC . EST. F 
LOGICAL  P= INT?  TUOPT? FIXCuF 

COMMON ''ABC-- X?  NSET.  NT  WORT.  NEIT1 .  MBITS.  NSECT «  N 
COMMON  'RAW-'U . V >  WE I GHT . M ?  KOUNT 
COMMON--' STEER  PRINT.  TUOPT.  FIXCOF 
COMMON  -' MACH  -  IND*  IOUTD 

S  SET  UP  MACHINE  CONSTANTS 

WRITE  a  ?  1  00.) 

100  FORMAT  <••'■’ WHAT  IS  THE  INPUT  LU«*M'- 
READ  a  <  ♦  >  IND 
W  *  I  T  E  1 . 1 .  c  U  0 

£00  FORMAT  <V" WHAT  IS  THE  OUTPUT  LU«?"> 

READ <1 ? ♦>  IOUTD 

l“ 

S  READ  IN  F- POE  LEM  DATA 
CALL  INDAT 

C  SET  SMALL  AND  LARGE  DELTA 

SDELTA  =  5D0'1  ♦♦MBIT  1 

DELTAS  =  i:.5B0>  ♦♦MBITS 

C  SET  RANDOM  NUMBER  GENEc'ATOc'  BY  CALLING  IT  NSET  TIMES 
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DO  10  J*1?NSET 

DUMMY  =  DBLEMJNI  <0>  > 

10  CONTINUE 

;  TURN  OPP  THE  PRINTING  OPTION  FOP  FIJNCT 

C 

PRINT  =  . FALSE. 

S  POUND  OFF  THE  COEFFICIENTS  TO  NEIT1  BITS 
DC  £0  J= 1 ? N 

X  <  J>  *  DS 1 GN  <.SDELTA»FLOAT  •'  I D I  NT  '..DABS  <X  <  J  •  >  ✓  SDELTA 
+  +.  5D0)  )  ?  X  l'J>  ) 

£0  CONTINUE  v 

UP  I TE  <  I  OUTD?  •?•?•?•?■>  HB I T 1 

9999  FORMAT  O' '37HTHE  INITIAL  COEFFICIENTS?  POUNDED  TO  ?  IS? 

+  1  OH  BITS?  APE) 

DO  30  J=1 ?  NSECT 

WR  I  TE  •'  I OUTD  ?  9998)  J  ?  X  <4*.j-  3;.  ,  x  ?  X  <4»J- 1 :>  .  X  '..4*  J) 

30  CONTINUE 

9993  FORMAT  <?HSECTI ON ?  1 3  - 1 H  ?  £D£ 0 .  1  0  '  1 H  ?  £D£ 0 .  1  0 ) 

Q 

C  SET  PARAMETER  VALUES  FOR  H-AND- J?  THESE  MAY  BE  CHANGED  IF 
C  DESIRED.  H-AND- J  MULTIPLIES  THE  STEP  SIZE  BY  PHO!  TAKE  FHO- 

PHO  =  .  SB0 

C  H-AND- J  STOPS  IF  THE  FUNCTION  VALUE  FALLS  BELCH  EST5 
S  TAKE  EST  =  0. 

EST  =  0. 

S  H-AND- J  STOPS  IF  HE  EXCEED  "LIMIT"  FIJNCT  CALLS! 

C  TAKE  LIMIT  =  10-000 

LIMIT  =  10000 
MR  ITE  •'  IOUTD?  9997') 

999 7  FORMAT  A7HNEXT  FOLLOWS  A  REPORT  FROM  THE  SEARCH  ALGOR  I 
CALL  HAND  J <N -  DELTAS?  SDELTA? PHO?  E  ST  ?  L I MI T?  X  > 

WP  I  TE  <  I  OUTD  ?  9993)  NB I T 1 

9993  FORMAT  K-3  0HTHE  FINAL  COEFFICIENTS?  HAVING?  13? 

♦  10H  BITS?  ARE) 

DO  40  J=l? NSECT 

WRITE1  I  OUTD  ?  9993)  J?  X 'S4*  J— 3)  ?X'-4*J— £>  ?  ’ !  \4*  J—  1 ''  -  sv4*J> 
40  CONTINUE 

•S  ®PINT  OUT  THE  FINAL  FINE  GRID  AND  QUIT 

PRINT  =  .TRUE. 

CALL  FIJNCT  <N?  X?  F> 

STOP 

END 

C - 

C  SUBROUTINE:  INDAT 
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S  THIS  SUBROUTINE  READS  IN  THE  PROBLEM  DATA: 

I  CARD  1  HAS  THE  NUMBER  OF  COEFFICIENTS  PER  STAGE  •S  0=  4... 

C  THE  INITIAL  SETTING  OF  THE  RANDOM  NUMBER  GENERATOR? 

■I  WHETHER  TWO-OPT  SHOULD  EE  USED  -0  IS  NO?  1  IS  YES;'-  THE 
C  DESIRED  BIT  LENGTH  OF  THE  COEFFICIENTS?  AND  THE  INITIAL 
S  SEARCH  DELTA  <IN  EITS>  FOR  THE  HOOKE  AND  JEEVES  ALGORITHM. 

C  THE  NEXT  SET  OF  CARDS  SPECIFIES  THE  GRID?  ONE  GRID  POINT 
S  PER  CARD:  EACH  CARD  HAS  FIRST  THE  FREQUENCY  IN  FRACTIONS 
S  OF  THE  NYQIJIST  FREQUENCY?  NEXT  THE  DESIRED  TRANSFER  FiJNC- 
C  T I ON ?  AND  THIRD  THE  TOLERANCE.  TO  INDICATE  THE  END  OF  THE 
C  GRID  POINTS.  SPECIFY  A  FREQUENCY  OF  1.  THE  M  CARDS  ARE 
S  COUNTED?  M.LE.100.  THE  PRESENT  SUBROUTINE  FIJNCT  ASSUMES  POP 
C  THE  A  CALCULATION  ''THAT  THE  DESIRED  TRANSFER  FUNCTION  TAKES 
C  ON  ONLY  THE  VALUES  0  OP  1.  THAT  IS?  THE  FILTER  HAS  ONLY  PASS 
S  AND  STOP  BANDS?  AND  EVERY  FASS  BAND  HAS  DESIRED  VALUE  OF  1? 

S  AND  the  same  tolerance. 

:  THE  next  card  has  the  number  of  second-order,  sections?  nsect. 

:  THE  NEXT  CARDS  HAVE  THE  "INFINITE  PRECISION"  COEFFICIENTS? 

S  ONE  PER  CARD.  THE  COEFFICIENTS  POP  SECTION  1  ARE  GIVEN 
S  FIRST?  NEXT  SECTION  E?  ETC.  WITH  IN  EACH  SECTION  THEY  ARE 
C  IN  THE  ORDER:  NUMERATOR  COEFFICIENT  OF  Z**-.-l"'?  NUMERATOR 
C  COEFFICIENT  OF  Z*-*-OE>  <IF  4  PER  SECTION'*  «  DENOMINATOR 
S  COEFFICIENT  OF  I**?::-l>?  DENOMINATOR  COEFFICIENT  OF  ;>♦  •'-£>  . 

C  THERE  ARE  3 ♦NSECT  OR  4*MSECT  COEFFICIENT:  READ  IN?  DEREN- 
S  DING  ON  WHETHER  THERE  APE  3  OR  4  COEFFICIENT:  =  ER  SECTION? 

C  NSECT. LE.?.  IN  THE  CASE  OF  S  COEFFICIENTS  C'ER  SECTION?  THE 
C  SECOND  NIJMEERATOR  COEFFICIENT  IS  SET  TO  1  AND  FROZEN. 


SUBROUTINE  IMDAT 

DOUBLE  PRECISION  i.lc'30')  ?Y‘SO«  ?  HEIGHT  1  SO  *  ? 

LOGICAL  R= INT. TMCPT?  PIXCOP 

COMMON  ■AB.C^X?.n:ET?  NTMOP'T?  NBIT1  ?  MBITS?  NSECT*  N 
COMMON  ■'  R'A'.l  '..I  ?  V  *  WE  I GHT  ?  M  *  ‘  .CUNT 
COMMON  STEER  'PRINT?  T'.lOPT?  FIX'COF 
COMMON  MACH  IND?  IO'JTD 

C  INITIALIZE  ‘■'O'JNT?  THE  NUMBER  OF  FUNCTION  EVALUATIONS 
OUNT  =  0 

C  READ  AND  JRITE  CARD  I  PARAMETERS 


READ'  IND?  ♦'  NCCEF  ?  NSET  ?  NTWCF'T  •  NB  IT  1  ?  MBITS: 

MR  I  TE  CIOUTD-  G9'?S"‘  NCOEF -  NSET?  NTWOFT ?  MBIT  1  ?  MBITS 
?'?■?•<  FORMAT'  ZSH  ♦♦♦♦♦  INPUT  DATA  ♦♦♦♦♦ '9HNUMBER  Oc 

*  1  EH  COEFFICIENTS  PE-10HR  STAGE  IS?  13- 

+  Z5HCAND0M  INITIALIZATION  IS  -IS  SHTWOPT  IS? 13’ 

+  4  SHCINAL  AND  INITIAL  ROUNDING  "•  PRECISION  IS? 

♦  I S  ?  TH  BITS  23H INITIAL  SEARCH  DELTA  IS?I3?5H  BIT 
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:  :et  logical  variable  tmgpt  ir  tmo-cpt  i:.  called  for 

T.JCPT  =  .  "ALSE. 

I- fNTWDPT.EQ. 1>  TWOPT  =  . TCUE. 

C  READ  AND  WRITE  THE  GRID 
MR  I TE  C I OUTD »  9997') 

•?•???  FORMAT  C- 19HGP ID  SPECIFICATIONS-'' 

+  £9HP0INT  NO.  FREQUENCY? 

+  4 OH  DESIRED  Y  TOLERANCE.'1 

M  =  0 

10  M  =  M+l 

READ  CIND?  9996)  W  CM:-  ?  Y  -M)  ?  WEIGHT  CM) 

999G  FORMAT  C3F 1  0.  6) 

MR  I TE  <1 OUTD  ?  9995)  M ?  W CM > ?  Y CM> ? WE  I GHT  CM) 

9995  FORMAT C 1 9 - 3D£ 0 . 1 0) 

WEIGHT (N)  =  1. DO -ME  I GHT CM) 

I F  CM  CM)  .  LT .  1 .  D 0)  GO  TO  1  0 

C  READ  AND  WRITE  THE  COEFFICIENTS 

READ  CIND? ♦  )  NSECT 
WRITE CIOUTD?  9994  -  NSECT 

9994  FORMAT  C-'S9H INITIAL  HIGH  PRECISION  COEFFICIENTS  rOR. 

*'  IS-  1SH  SECTIONS  ARE) 

N  =  4  ♦NSECT 

I F  NCOEF  .  ECi  •  3)  GO  TO  £0 
FIXC OF  =  .FALSE. 

READ  C  IND?  999;.)  <X  CJ)  ?  .J=  1  ?  M) 

FORMAT CF£0. 16) 

GO  TO  4  0 
FIXCOe  =  .TRUE. 

READ  CIND?  9993)  C.<C4*J-3)  *  X  C4^  J-l)  ?XC4*J)  ?  J  =  1  ?  NSECT) 

DO  30  J=l« NSECT 
X C4*  J— £ )  =  1.D0 
CONTINUE 
CONTINUE 
DO  50  J=l? NSECT 

MR  I  TE  C I  OUTD  ?  99 9£)  J  ?  X  <«*♦  J-3)  >  X  C4*  J-£ '  ?  I  <4*  J-l)  ?  X  C4  *  J) 
CONTINUE 

FORMAT  C7HSECT I  ON  ?  1 3  1 H  ?  £B£  0 .  1  0-'  1 H  ?£B£0.  10- 
RETURN 
END 


r*  9  3 
£0 

SO 
4  0 

50 

999  £ 


C  SUBROUTINE:  FUNCT 

3  THIS  COMPUTES  THE  VALUE  OF  THE  ERROR  FUNCTION?  USING  THE 
C  GRID.  IF  PRINT. EQ.  .TRUE.?  IT  COMPUTES  THE  ERROR  FUNCTION 
C  AT  A  FINER  GRID  C  10  POINTS  FOR  EACH  GRID  POINT'?  AND 
S  PRINTS  THE  RESULTS?  THIS  IS  USED  TO  PRINT  OUT  THE  FINAL 
C  TRANSFER  FUNCTION. 

C  NOTE  THAT  THE  CALCULATION  OF  THE  CONSTANT  MULTIPLIER  A 
C  ASSUMES  THAT  ALL  THE  PASSBAND I  HAVE  THE  SAME  DESIRED  VALUE 


•:  OF  1.  THIS  MUST  EE  CHANGED  TO  A  BAND-BY— BAND  CALCULATION 
c  IF  THERE  IS  MOPE  THAN  ONE  FASCBAND  DESIRED  VALUE. 


SUBROUTINE  F1JNCT  <N«  X*  F> 

DOUBLE  PRECISION  W  030  *  V  030  .  WEIGHT  •SO'1 
DOUBLE  PREC I  S I  CM  COSTS  03  0>  .  CQSST4  •:  3  0> 

DOUBLE  PREC I S I ON  SUE 1  03 O -  SUES  03 0 > . SUES  o=  o> 

DOUBLE  PRECISION  SUE6O30>  ,SUB?<30>  .  SUB3<30> 

DOUBLE  PRECISION  X<33>  *  YHT<30>  »  E<30> 

DOUELE  PREC  I S  I  ON  YHTHT  <  1D0.11  :• 

DOUBLE  PREC  I  S I  ON  ,A.  P » P I  .•  NIJM «  DEN .«  Y 1 .«  YE .«  FCEQ «  CCC  .«S3S>fil.>ER 

LOGICAL  PRINT. TWOPT. FIYCOF 

COMMON P fl W / W » '  V «  WEIGH*.  M.  KCUNT 

COMMON ''STEER -"PRINT.  TWOPT.  FIXCGF 

COMMON-'MRCH I ND .  I OU TD 

p 

C  CHECK  THAT  THE  POLES  OF  THE  FILTER  ARE  INSIDE  OR  ON  THE 
C  UNIT  CIRCLE;  no  CHECKS  APE  MAKE  ON  NUMEEATOP  CCECFICIENTS 

K  =  N  4 
DO  10  J=1 >  K 
J4  =  J-*4 

IF  <X  OJ4.)  .  LE.  1 .  DO.  AND.  1 .  DO+X  <  J4  ■  .  GE.  DABS  <X  ■  J4- 1  >  •  > 

+  GO  TO  10 
-  =  1.D1S 
RETURN 
10  CONTINUE 

=  1  =  4 . D  0*  DA T AN  < 1 . D 0> 

I F  ‘'KQUNT .  ME .  0>  GO  TO  30 

C  COMPUTE  AND  CAVE  THE  TRIG  FUNCTIONS  AT  THE  G- ID  POINTS 
C  THE  FIRST  TIME  THROUGH 

DO  30  1=1. M 
COSTS < I • 

C0SST4  ■:  I) 

30  CONTINUE 

*  MOVE  THE  CALCULATION  OF  CONSTANTS  OUTSIDE  THE  INNER  LOOP 

30  DO  40  J=1.K 

J4  =  J*4 

SUE  1  ■' J  1  =  X  <J4-3>  ♦♦£+■  <X  CJ4-3'1  -1 .  DO  1  ♦♦3 
SUES  1 J >  =  X ( J4— 3> *  0> ' J4-3>  +1 . D0> 

SUB  3  <  J  ■'  =  X  04-3  > 

SU B -5  <  J >  =  X  <  J 4 - 1  >  < X  < J 4  '  - 1 .  D  0 >  ♦  ♦  3 

SUE?  <J>  =  X  0J4-1 )  ♦  <X  <J4>  + 1 .  D O  ' 

SUB '3  <S>  -  X  <  J4 ' 

40  CONTINUE 

l” 

C  EVALUATE  THE  MAGNITUDE  Oc  THE  TRANSFER  FUNCTION.  YHT ? 

S  AT  EACH  OF  THE  M  GRID  POINTS 

DO  GO  I = 1 . M 


=  3 .  D  0 *DCOS  1  P I  ♦  W  I )  ■ 
=  C  0  S  T  3  <  I  ♦  ♦  3 


97 


MUM  =  1.20 
DEM  =  1 . DO 
DO  "50  J=1«K 

MUM  =  MUM*  USUEl  +TUEE  <  J  1  ♦COSTE  •  D  I  JE 3  •'  J>  ♦CDT2T 4  I >  • 

DEN  =  DEM*  CSUB6  <  J>  ♦  i'JE?  •"  J '  ♦COSTS  f  I  >  fS-JB$  <  S>  ♦<! OS£T4  <  I )  • 

TO  CONTINUE 

YHT  <I>  =  DCQ^T  ''NUM  -  DEN-' 

CO  CONTINUE 

C  FIND  THE  LARGEST  AND  SMALLEST  VALUES  OF  VHT  IN  THE 
C  FASSBANDS 

Y 1  =  1 . DIE 
YE  =  0.D0  ' 

DO  70  I  =  1 «  M 

I F  Ci'  ■.  I  ■  .  EQ.  0.  D0>  GO  TO  7  0 
VI  =  DM  INI  '.Y 1  ?  VHT  <  I  '>  ’> 

Y£  =  DMAX1  CY2.  VHT  ■  I'>  ,• 

70  CONTINUE 

DEFINE  THE  CONSTANT  MULTIPLIER  A  TO  EE  THE  RECIPROCAL  ^F 

:  tj£  average  of  the  smallest  *nd  largest  yht't.  im  the  fa::- 
C  EAMDS.  THIS  MAKES  DENSE  IF  ALL  THE  FASSBANDS  HAVE  THE  SANE 
I  DEC  I FED  VALUE  OF  1  RND  THE  SANE  TOLEc*NCE;  OTHEFHIIE  THIS 
C  CRLCULRTION  OF  R  MUST  EE  MODIFIED 

c  * 

h  =  d.  Du  cfi+Y*) 

.  r 

I-  CRLCULRTE  THE  WEIGHTED  E=  RGR  RT  THE  GRID  -GIMTS 

F  =  0 .  D  0 
DO  3  0  1  =  1  jM 

V  HT  1  I  ■  =  A*  YHT  '  I  ■ 

E  '  I  •  =  WEIGHT  ■  T’  ♦  DAB!  '  YHT  '  I  '•  -V  1  I  ■ 

F  =  BMAX1  (P.E(D 
50  CONTINUE 

FOUNT  =  FOUNT  + 1 
IF  NOT.  PRINT'1  CET'JRM 

C  I -  REQUIRED <  PRINT  OUT  THE  RESULT!:  CJr4  THE  CORF  IE  GRID 
UIF I TE  <  I OUTD*  9999>  R .  F 

9999  cOc'MAT  <■,'/”♦♦♦♦♦  c  I NRL  RESULT!  ♦♦♦♦♦"  17HTHE  CONCTRNT  R  1 1  • 

+•  DEO.  1 0'-3cHcUMCTIOM  VALUE  OM  COARSE  GRID  I !  •  DEO.  10  ' 

WR'ITEaOUTD.  9993  > 

9993  eORMRT  07  V  "FREQUENCY"  >  7X>  9HDE SIRED  Y'iG  EHRCTURL  V *> 

*  1 1 : i «  5HEFP0R':1 

DO  90  1  =  1-  M 

I.P  I  TE  ( I OUTD .  9997  >  i.l  <  I  >  .  Y  < I  >  «  YHT  <  I )  «  E  <  I  ’■ 

90  CONTINUE 
9 r 97  FOFMRT  < 4D 1 C . 3 1 

P 

C  TO  RE-CRLC IJLRTE  R<  CRLCULRTE  N-'D=YHTHT  ON  THE  rINE  G=IB 
C  THE  SAME  RESERVATION  ABOUT  THU  CRLCULRTION  Oc  R  AR='LIEC 

:  a:  above 
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.'I  =  1.012 
V£  =  0.  DO 
MM  =  M— 1 
DO  130  1=1- MM 

c 

:  SKI®  TRANSITION  BANDS.  DEFINED  AS  THOSE  WHERE  VCD  CHANG 

C 

IF(V(I)  .ME.  Y<I+1"N  GO  TO  130 
DO  120  11=1.11 

IFai.NE.il)  '30  TO  100 

3  find  cases  where  the  last  point  of  the  band  mill  mot  be 

:  COMPUTED  LATE=- 

IFa.EQ.MM>  GO  TO  100 
I F  < Y  a  + 1  :• .  ME .  V  C I  ♦£>  >  GO  TO  1  0 0 
GO  TO  120 

freq  =  w  a>  ew  a  d  -  w  a>  >  ♦ .  ido^floatcii-d 

CCC  =  £. D0*DCQC CPI*FREQ> 

CSS  =  C  C  C ♦ ♦ 2 
NUM  =  1.D0 
DEM  =  1.D0 
DO  110  J=1 
J4  =  J*4 

MUM  =  HUM*  CX C  J4-3>  ♦♦£•*•  CX J4-£ >  -1 .  DC"  **£+X  J4-S 
+•  CX  c  J4-2>  + 1 .  D 0>  ♦CC C>X  <  J4-£>  ♦  1C  l"1 

DEN  =  DEN*  CX  C  J4-1  >  ♦♦£•*•  CX  C  J4>  - 1 .  D0>  **2*X  c  J4- 1  .  ♦ 
+  <  j4>  + 1 .  DO'1  *CCC  +X  C  J4  •  ♦  Z  3S> 

CONTINUE 

VHTHT  C  I,  ID  *  D  CQRT  CNUM.-'DEM) 

I F  C V  CD  .  EQ.  0.  DO1  GO  TO  120 
VI  =  DMIN1  'YD  VHTHT  CD  ID  • 

Y£  =  DMA! ;  1  <Y2 .  VHTHT  'DID'* 

CONTINUE 
CONTINUE 
A  =  2. DO  •  CY1 *Y2> 

S  NOM  REPEAT  THE  CAME  LOOP'S.  FINDING  THE  ERR'C®  -ND  c-®iMT:M 
C  THE  RESULTS 

yR  I  TE  *•  I  DU  TO .  33'rS  1 

333-3  FORMAT  (✓/'“FINAL  RESULTS  ON  FINE  GRID-  10  POINTS  =E=  " 
*  UH  GRID  POINT- 

MRITECIOUTD. 3333> 

F  =  0.  DO 
DO  1  GO  1  =  1. MM 

I F  <y c I > .ME. V  <1  +  1 >  >  GO  TO  1 6 0 
IP  =  1+1 

WRITE  CIOUTD.  3335-  DIP 

3335  FORMAT (”♦♦♦♦♦  BAND  frqn  GRID  POINT". 15. CH  TO* 

+  15. GH  ♦♦♦♦♦> 

DO  150  11=1.11 

IF  a  I.  ME.  ID  GO  TO  140 
IFCI.EQ.MM>  GO  TO  140 
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IF •va  +  l>  .ME. V-  53  TO  14  0 

50  TO  150 

1 4 o  fpeq  =  w  a >  +  *:m  •  i  +  i  >  -m  •  i  .■  ■♦.id o^rloat  •  1 1 — i  • 

Hi  =  n*VHTHTa  .  1 1  > 

ERP’OR  =  ME  1 5HT  <  I  >  ♦DAB  S  <  ft  1  - V  '■  I  )  1 
F  =  DflftXl  <F.  EFFORT 

MP I TE  <  I  CUTD»  9?97>  FREQ .  V  < I  >  -  ft  1 .  ERROR 
150  CONTINUE 
160  CONTINUE 

MR I TE  < I DUTD .  9994>  fl.R 

'?'?  ?4  FORMAT  '"THE  VftLUE  OF  ft  FROM  THE  FINE  GRID  IS"*  DEO.  10 

+4SHTHE  FINftL  VftLUE  OF  THE  ERROR  ON  THE  -INE  GRID  IS.  DEO.  10 
RETURN 

END  v 


C  SUBROUTINE:  HftNDJ 

C  THIS  IS  ftN  INRLENENTftT ION  OF  ft  RANDOMISED  HOOKE  ftND  JEEVES 
:  SEARCH  ALGORITHM.  IT  FINDS  LOCAL  OPTIMA!  THE  COORDINATE:  ARE 
C  RANDOMLY  ORDERED  EACH  TIME  EXPLR  IS  IHVOl ED.  FLOW  CHARTS 
C  FO»  THIS  SUBROUTINE  AND  THE  NEXT  ARE  INCLUDED  IN  THE 
C  REFERENCED  PAPER . 


•SUBROUTINE  HAND  J  ON  *  DELTAS «  S  DELTA  <  R  HO  •  ^  IN. LIMIT. RSI • 
DOUBLE  PRECISION  RSI ' B6>  .  THETA  >:li  <  .  ®wj  •:  l  i 
DOUBLE  PR EC  I  SIGN  W  KB  0 1 . r  020 • . WEIGHT  020 > 

DOUBLE  PRECISION  FR* I . RRH I . DELTA. RHC* DELTAS. "DELTA. FMIN 
COMMON  PftM  -  M  .  Y .  ME  I  GMT.  M .  »>  DUN T 
;  OMMGN ■-'MACH-  I  NO?  I OUTD 
CALL  FUNCT  <ti  .RSI.  FRS  I'* 

DELTA  =  DELTAS 
MR  I  T£.;  I  OUTD. 

?A?A  F-CMAT<"NO.  CALLS  c'AT.  SIZE".  UN.  "DELTA",  r  '. 

♦  "OLD  VALUE"  "NEW  VALUE” > 

10  IF  >: DELTA.  LT.  SDELT'ft"'  50  TO  50 
IF' -PSI.LT.FMIN>  SC  TO  50 
I  f.'KG'JNT.GT.  LIMITS  50  TO  50 
CALL  SET  '."RH I .  R  S I «  N> 

R='H  I  =  FRS'  I 

CALL  EXPLP  <PH  I .  RRH  I .  M.  DELTA') 

I  =  0 

BO  IF-rFPHI.GE.FPSI)  50  Tu  4  0 
IF<F»SI.LT.FMiro  50  TO  40 
IF •rKQUNT.GT. LIMIT)  50  TO  4  0 
I  =  1*1 

HP  I  TE-' I  OUTD.  9999-  V'GUNT.  I .  DELTA.  FRS  I .  FPHI 
999  B  FORMAT  CBM. .16. 5X .  1 6 .  BD 1 6 . 2 > 

C  ALL  CcT  <THETA . PS I . N> 

CALL  SET 'RSI. PHI. N> 

F 5  S I  =  FPH I 
DO  BO  J=  1 .  N 

PHI  r  J)  =  £ .  ♦bh I  J >  -  THETA  < J> 

BO  CONTINUE 

CALL  f;JMC T  >:'M.  PH  I .  C-H I 
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CALL  E.YPLP  -  PHI «  CCHI «  M.  DELTA- 
SO  TO  20 

4  0  I c  a.  ST.  CO  SO  TO  10 

DELTA  =  PHO*DELTA 
50  TO  10 

50  WR I TE  •! I OUTD  ?  YOUNT- FPSI 

999?  FORMAT  FINAL  VALUES-  NO.  OF  RiJNCT  CALLS*" - IS  - 
+  7H  FUNCT*»D20. 10> 

RETURN 

END 


C  SUBROUTINE:  EYPLR 

C  THIS  IS  A  LOCAL  E>;PLOPATION  SUBROUTINE  USED  BY  H-AND-J 


SUBROUTINE  EXPLP  <RHI  ,  FPHI >  N*  DELTA- 

DOUBLE  PPEC I S I ON  PH  I  - 3S- -  FPH I - DELTA , S  AVE - FNEW -  SAVE I  -  SAVE J 
INTEGER  LIST  CS*;- 

COMMON  S  T  E  E R ■ ' PRINT*  TWOPT «  FIXCOF 
LOGICAL  PRINT*  TWOPT? PI XCOF 

C  RANDOMIZE  THE  ORDER  OF  THE  COORDINATES 


CALL  CHUFF  <N» LIST;- 
I COUNT  *  1 

10  IF  Cl  COUNT. ST.. N>  RETURN 

i  =  li  st i  count:- 

C  IF  ONLY  THREE  COEPF ICIENTC  WERE  READ  IN  PER  STAGE 
S  THEN  DO  NOT  CHANGE  THE  SECOND  NUMERATOR  COEFFICIENT 

V SAVE  =  I--  I  4;-  *4 

Ir  -  FI'CSOF.  AND.  <}■  ."SAVE.  EC!.  2'-  >  SO  TO  100 

SAVE  =  phi  a;- 

C'HI  I  :>  =  ='HI  a.:-  +  DELTA 

SALL  rUNCT  < N -  PHI  ~NEW> 

IF-iFNEW.  SE.  FPHI)  SO  TO  20 
CPHI  =  FNEW 
SO  TO  100 

20  PHI  -:'r>  =  =HI  -I':-  -2.  D0*DELTA 
CALL  RUNCT  <N, PHI • FNEW> 

I p  <PHEW . Sc . FPH I >  SO  TO  30 
PR HI  =  FNEW 
SO  TO  100 

SO  PHI  a  -  =  SAVE 

IP  <•  NOT.  TWGPT;-  SO  TO  100 

l“ 

•S  BEGINNING  OF  2-OPT 

r 

JCOU.NT  =  I  COUNT  + 1 
40  IF  <  JCOUNT.  ST.  N':-  GO  TO  100 
J  =  LIST  -  JCDUNT'- 
C 

C  SIMILARLY  FDR  TWOPT-  I-  THREE  SZSF-ICIENTS  DEP  STAGE  WERE 
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C  READ  IN.  DC  HOT  CHANGE  THE  SECOND  NUMERATOR  COEFFICIENT 
KSAVE  =  J-CJ'"4-*4 

IF<PIXCOF.  AND.  CKCRVE.EQ.  ,;q  JO  90 
SAVE  I  =  PHI  a> 

SAVE  J  =  PHI <J> 

PHI  a  ’  =  PHI  < I '■  +  DELTA 
PHI <  J!>  =  PHI <J'  +  DELTA 
CALL  FUNCT  <N .  PH  I .  FNEU> 

I F  CFNEU .  GE .  FPH I  >  30  TO  50 
PPHI  =  FNEM 
30  TO  90 

50  PHI  <1:-  =  PHI  CI>-£.DO*DELTA 
CALL  CUNCT  CM. PHI .  FNEM:- 
I F CFNEU . GE . FPH I >  GO  TO  6 0 
FPH I  =  FNEM 
GO  TO  90 

GO  PH  I  >:  J>  =  PH  I  <  J>  -2 .  D  0* DELTA 
:  ALL  p’JNCT  lN.  PH  I .  FNEW> 

Ir  <FNEM.  GE.  FPH  I'1  GO  TO  70 
PP'H  I  =  FNEM 
GO  TO  90 

70  PH  1(1)  *  PHI  < I >  +2 .  E O'* DELTA 
CALL  FUNCT CN .  PH I , FNEM> 

IF  <FNEW. GE. F»HI>  GO  TO  SO 
FPH I  =  PNEM 
GO  TO  90 

SO  PHI 1  T>  »  CAVE  I 
*  CAVE  J 

?0  J  COUNT  =  JCOUNT-t-1 

GO  TO  40 

C  END  OF  E-OPT 

100  ICOUNT  =  MOUNT*! 

GO  TO  10 
END 


C  SUBROUTINE:  SET 

C  THIS  SETS  VECTOR  A  EQUAL  TO  VECTOR  B 


SUBPOUT I NE  SET  < A .  E«  N> 

DOUBLE  PR  EC  I  SION  A  C  3*S>  >  E  CSS  > 
DO  10  I=l< N 
A  C I)  =  B < I > 

10  CONTINUE 
RETURN 
END 


C  CUBPOUTINE:  CHUCC 
C  THIS  c ANDOMLV  ORDERS  LIST 
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SUBROUTINE  SH'JPR  '.N-  Li  ST  • 
INTEGER  LIST  <'36) 

DO  10  1=1 j N 
LIST  <I>  =  I 
10  CONTINUE 

HO  £0  LL= 1 -  N 
L  =  N-LL+1 

j  =  i  nt  -:rloat  <l>  ♦un  i  •:  o:>  >  ♦  i 

K  =  LIST  >'L> 

LIST  <L>  =  LIST  CJ> 

LIST  <J)  =  K 
£0  CONTINUE 

RETURN  v 

END 


C  RUNCTION:  UNI 

C  AUTHOR:  — I L I TH  M.  GROSS 

C  BELL  LABORATORIES-  MURRAY  HILL*  NEM  JERSEY 

C  07?  7  4 

C  PORTABLE  RANDOM  NUMBER  GENERATOR 


FUNCTION  UNI -MO 
INTEGER  I  BYTE  -;4> 

DATA  ICS  EEC -'O'- ITSEED  0-- IFCN si' 

UN  I  =R  1 UNF  <.  I C  SEED «  ITSEED  - 1  BYTE  >  IFCN> 

RETURN 

END 

FUNCTION  R 1  UNF  ■'  I CSEED  -  ITSEED-  I  BYTE-  I-CN'- 
LOGICAL  “IRST 


INTE 

SEP  C  SEED  ■'  £>  »  TSEEI 

<"££"•  -  xor  . 

£'?•  -  I BYTE 

•  4 r  I  'I.  C  F  '  •' 

DATA 

xcRa>''  i''-:  :op<£>  ' 

£  •  -  XOR  <  2  '• 

S  -  XOR  '  4 

-  ■  S  <  XD»  < :s>  • 

1 

XOR 

’  G  '<  / 1 ''  -  XOR  ( 7-'  '4 ''  •  ',■■■ 

OR <S>  - 

XOR  C?  •  i-' 

XOR < 1 0 ) s 7/  - 

z. 

XjR 

ai  ■  5  '  -  XOR  •:  1  £>  '"4 

-XOR  •IS).'' 

7  '  • '  !OR  ■  1 4 

■S$  s ,  XOR  <  1 

1 

T; 

XOR 

(  1 6?  6/’-  XOR  <  1 7>  '7s 

•  XOR  •.  1  S') 

4  <  ;OR 1  l '? 

•  '5  '  •  XOR  <2 0> 

'  £  ■* 

4 

XOR 

<  £  1  >  '  2  '-XOR  •'££  »  '  7 

-  XOR.  C23')  ' 

G'-XOR <24 

■  •  j',  XOR <25> 

•4 

5 

XOR 

''£6.''  s S't  XOR  <27>  £  ' 

-  XOR  C2S>  ■ 

ix-  :op  <:.£•? 

/•'  0  s 

DATA 

FIRST''.  TRUE.  -  JCS 

EED  '1£345 

-  JTSEED-' 

1  073 

DATA  CSEED'' 0-  0-  0-  0-  0-  O' 

DATA  T  S  E  E  D  0  -  0  -  0  -  0  -  0  -  0  -  0  -  0  -  0  -  U  -  U  *  U  -  0  -  0  -  0  -  0  -  0  -  0  -  0  •  0  -  0  -  0  - 
1  0  -  0 »  0  -  0  -  0  -  0 »  0  -  0  -  0  -  0 
R1UNF=0. 0 

IF  ■'  <.  NOT.  CIRST  '  .  AND.  <IFCN.  3T.  0>  •  GG  TO  50 
IF  <  IFCN.  GT.  O'*  GO  TQ  to 


JC  SEED®  I ABS  I CSEED> 
JTSEED=IAES  ^ ITSEED) 

io  f:rst=.ralse. 

c 
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c:eed<i>  =  JCSEED 
DO  20  1=1.5 

CSEED  a  +  1'-  =C  SEED  <  I  ■  ''64 
cseed  a.- =c  seed  -n -cseed  <  i+i  -  +64 
CSEED ■  6)  =MOD  <CSEED  -,6>  « 4  > 


IF  <  < JCSEED. NE.  0'>  .AMD.  <MOD  '-.CSEED  -i:>  .  2> 
+  C SEED  < 1 > =C SEED < 1 > + 1 

TSEED  < 1 >  * JTSEED 
DO  30  1=1.11 
TSEED  <1  +  1  >  =TSEED  (I)  '2. 

TSEED  •  I  ■  =T3EE0  -  D  -TSEED  <1  +  1  ■  +2 


DO  40  I =1 2. 3c 
TSEED • I > =0 

IF <JTSEED.  NE.  0>  TSEED -  l.»=l 
I-  <IFCN.E3.  0>  RETURN 
CONTINUE 
DO  60  1=1.17 

T  SEED  <  I  >  =  I  AES  •:  TSEED  •;  I  >  -T  SEED  '1  +  15  -  - 
DO  70  1=13.32 

TSEED  -'I.-  =  1  AES  -CTSESD  <D  -TSEED  -II -17'-  - 

CSEED  <6  =1 3+CSEED  -  6>  +55 +C SEED  <5>  +16+C 
C S EED  <5 -  =  1 3+C SEED  -  5'-  +55+C SEED  -  4 :■  ♦  1 6+C 
CSEED  -  4  -  =1  3+C  SEED  1  4  -  +55+CSEED  -'  S'-  +16+C 
C  SEED  -.3.)  =1  3+CSEED  C3  >  +55+C  SEED  <2>  +16+C 
C SEED  1 2  >  = 1 3+CSEED  <2> +55+C SEED < 1 > 

CSEED <1 > =1 3+CSEED <1 > 


ED  <4 
ED  -■  3 


IChPFY=0 
DO  30  I=l«5 

K  =  k'+4 

C  SEED  ( I =C  SEED  <  I  -  + 1  CA*RY 

I  SAPPY* CSEED  -  I.-  ''64 

CSEED -  I  - =CSEED -  I  -  — 6-4  +  ICARRV 
1 2=CSEED  < I >  '3 

I I  =C"SEED  <  I  '■  -3+ 1 2 

J1=4+TSEED  <K'+2."-  +TSEED  <K+1>  +TSEED-T  +1  -  +T  SEED 
Jd=4+ T SEED  1 K  +  5 +T  SEED  < K  +4  -  +T S EED  <1  +4  -  +T  SEED 
I T 1 =23 

IFCI1.GT.  Jl'-  ITl  =  ai  +  Il-U>  '2+  J 1 
IC<I1.LT.  Jl>  I T 1  =  •'  J 1  ♦  J 1  -  J 1  >  2+11 
IT 2=23 

IF  -'12. 3T.  J2>  IT2=  <I2+I2-I2>  ••'2  +  J2 
IF  -  I2.LT.  J2>  IT2=-'J2+  J2-J2  -  ''2+12 
I  SCR  -  I  >  =3+>iOR  <  I T2+1.-  +XOP  •  I T 1  + 1  > 

R 1UNF=  -  » 1 UNF+FLOAT  aSCR  •  I  ■  >  -  64 . 0 
CSEED  <6^  =MOD  ■'CSEED  <6>  + 1  CARRY.  4'- 
J 1  =  TSEED  <  31  -  +TSEED  <32>  +TSEED  -.  32  .- 
IT1=CSEED  <6> -Jl 
IFaTl.LT.  0)  IT  1  =  0— IT1 

IF  <  < I T 1 .  EC .  1>  .  AND.  <CSEED  <6 -  +  J 1 .  EC- .  3'-  >  I T 1  =  3 
R  l'JNF=  -  P 1  UNP+cLOAT <111-  -  ''4.  0 
IF-:  IFCN.EQ.  l  -  RETURN 
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I  BYTE  (.4  <  sirCP  <1>  ♦MQIj  *■  Tl  C-  .£  *  4>  *64 

I  BYTE  '.'3)  *  I  SCR  <’£}  '4+MQD  1  I  7C~  <3 )  »  1G>  ♦IE 

I  BYTE  I  *  CP  C3>  -'16*1  SC*  <4>  *4 

I BYTE < 1 >  *  I :CR <5> ♦ I T 1 ♦G4 

RETURN 

END 

BLOCK  DATA 

COMMON  AABC'X .*  NSET >  NTWOPT*  MBIT  1 «  MBITS.*  MSECT  j  H 
COMMON  •'RAW  '!,J ,  V » WE  IGHT  *  M *  KOUNT 
COMMON. 'STEEP 'PRINT *  TWOPT, FIXCOF 
COMMON 'MACH '  I N  D  .*  I OUTL 

DOUBLE  PREC I  SION  W <3 0) *  V <8  0> *  WE I GHT  <3 0> *  X <36> 
END 


1 
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FWIIR  Input 
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Appendix  E:  A/D  Converter  Software  Listings 


Two  methods  of  software  support  for  the  A/D  conver¬ 
ters  were  implemented.  A  third  method  was  only  partially 
implemented.  Program  ANDIG  uses  a  polling  loop  to  acquire 
data.  It  can  call  one  of  two  subroutines,  AD14B  or  AD15B, 
depending  on  which  I/O  channel  the  user  selects.  These 
two  subroutines  initiate  the  A/D  conversion,  load  the 
data  value  into  the  A  register,  and  then  store  it  in  a 
memory  buffer  location.  This  process  repeats  until  the 
buffer  is  full.  The  subroutine  WAIT1  is  called  to  wait 
in  a  loop  a  specified  amount  of  time  until  the  initiation 
of  the  next  A/D  conversion.  An  alternate  implementation 
of  ANDIG  uses  the  subroutines  AD14  and  ADI 5  (same  entry 
points:  DT14B,  DT15B).  These  subroutines  were  used  to 

measure  the  approximate  conversion  times  of  AD14B  and 
AD15B. 

Program  ANDG2  uses  a  DMA  channel  to  acquire  data. 

It  can  call  one  of  two  subroutines,  DMA14  or  DMA15, 
depending  on  which  I/O  channel  the  user  selects.  These 
subroutines  check  for  the  availability  of  DCPC  channel  1, 
initiate  the  data  transfer,  and  wait  in  a  loop  until  the 
transfer  is  complete.  It  is  possible  to  execute  other 
instructions  rather  than  waiting  in  the  loop. 

The  three  dummy  drivers  listed  are  identical  in 
function  and  differ  only  in  the  names  of  their  respective 
entry  points.  The  listings  provide  sufficient  explanation 
of  their  purpose  and  implementation.  They  were  included 
in  the  latest  generation  of  RTE-III. 
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FTN4, L 

C 

C 

PROGRAM  ANDIG 
C 

C** *********************************** ******************** 


c*  * 

C*  THIS  PROGRAM  TESTS  THE  SUBROUTINES  AD14B  AND  * 
C*  AD15B.  THE  ENTRY  POINTS  ARE  DT14B  AND  DTI  SB ,  * 
C*  RESPECTIVELY.  * 
C*  * 


C* ************* ************* ******* *********************** 

C 

INTEGER  DELAY , DATA , BUFL , DECIDE 
COMMON  DELAY, DATA (8 192) ,  BUFL 
C 

C -  INITIALIZE  VARIABLES 

C 

WRITEd  ,  100  > 

READd,*)  DELAY 
UR ITE (1,200) 

READd,*)  BUFL 
WRITEd, 300) 

READd,*)  LUNUM 
WRITEd, 400) 

READd,*)  DECIDE 
C 

C -  GET  THE  DATA 

C 

IF(  DECIDE.  F.Q.  14)  CALL  DT14B 
IF (DECIDE . EQ . IS)  CALL  DT15B 

C' 

C— -  OUTPUT  DATA 
C 

WRITE (LUNUM ,S00 )  (I,DATA(I) ,1=1 , BUFL ) 

C 

100  FORMAT (/////, "CHOOSE  A  DELAY  VALUE  FROM  0  TO  32767.  ") 
C 

200  FORMAT(//,  "HOW  MANY  DATA  POINTS?  (B192  MAX)  "  ) 

C 

300  FORMAT(//, "WHAT  IS  THE  LU*  OF  THE  OUTPUT  DEVICE?  ") 

C 

400  FORMAT(//, "WHICH  I/O  CHANNEL,  14  OR  15?") 

C 

S00  F0RMAT(22(/) ,128(4(217)  ,/)  ) 

C 

STOP 

END 
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ASMS  ,  L 
* 

* 

NAM  ADI  4B 

* 


*  * 

*  THIS  SUBROUTINE  USES  POLLING  SOFTWARE  TO  * 

*  GET  DATA  FROM  AN  A/D  CONVERTER.  THE  I/O  * 

*  SELECT  CODE  IS  i4.  MAXIMUM  NUMBER  OF  DATA  * 

*  WORDS  IS  9192.  * 

*  * 

*  VARIABLES:  * 

*  DELAY  -  ADDS  DELAY  TIME  BETWEEN  THE  ACQUI-  * 

*  SIT  ION  OF  EACH  TWO  DATA  POINTS  * 

*  DATA  -  DATA  BUFFER  * 

*  BUFL  -  LENGTH  OF  THE  DATA  BUFFER  * 

*  # 


M 

EXT  $LIBR ,$LIBX,WAIT 
ENT  DT14B 

* 

COM  DELAY ,DATA<8192) , BUFL 

* 

DT14B  NOP 


JSB 

HI..IBR 

* 

TURN  OFF  INTERRUPTS  AND  MEMORY  PROTECT 

NOP 

* 

LDA 

DELAY 

* 

3SA 

* 

SKIP  IF  DELAY  >  0  OR  =0 

CMA 

* 

COMPLEMENT  IF  DELAY<0 

STA 

DELAY 

* 

LDX 

BUFL 

* 

SET  COUNTER 

* 

LOOP 

LI)B 

STORE 

* 

SET  BUFFER  POINTER 

ST  C 

1 4  B  ,  C 

* 

SEND  ENCODE  SIGNAL 

SFS 

14B 

* 

IS  DEVICE  BUSY? 

JMP 

*-i 

* 

YES.  TRY  AGAIN. 

LI  A 

14B 

* 

NO.  READ  DATA. 

STA 

STORE, I 

* 

STORE  DATA  IN  BUFFER 

INB 

* 

STB 

STORE 

* 

SET  POINTER  TO  NEXT  BUFFER  LOCATION 

JSB 

WAIT 

* 

DELAY  PERIOD  BEFORE  GETTING  NEW  DATA 

DSX 

* 

DECREMENT  COUNT.  SKIP  AFTER  LAST  DATA  PT 

V 

JMP 

LOOP 

* 

GET  NEXT  DATA 

'P 

CLC 

14B ,  C 

* 

CLEAR  I/O  CONTROL  BIT 

JSB 

$LIBX 

* 

TURN  ON  INTERRUPTS  AND  MEMORY  PROTECT 

<4t 

DEF 

DT14B 

* 

RETURN  TO  CALLING  PROGRAM 

STORE 

* 

DEF 

DATA 

* 

DEFINE.  DATA  BUFFER  ADDRESS 

END 

DT14B 

ASMB,  L 


NAM  ADiSB 

* 

^  ^  ^  ^  ^  ^  ^  ^  viy  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  sW  sw  ^  o>  ^  sW  ^  ^  ^  Oj  ^  ^  sw  ^  ^  ^  ^  if  u«  ^  u«  ^  ^  y  y  ^  y  y 

•v  *  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^n  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 


*  * 

*  THIS  SUBROUTINE  USES  POLLING  SOFTWARE  TO  * 

*  GET  DATA  FROM  AN  A/D  CONVERTER.  THE  I/O  * 

*  SELECT  CODE  IS  iS.  MAXIMUM  NUMBER  OF  DATA  * 

*  WORDS  IS  3192.  * 

*  * 

*  VARIABLES:  * 

*  DELAY  -  ADDS  DELAY  TIME  BETWEEN  THE  ACQUI-  * 

*  SIT  ION  OF  EACH  TWO  DATA  POINTS  * 

*  DATA  -  DATA  BUFFER  * 

*  BUFL  -  LENGTH  OF  THE  DATA  BUFFER  * 

*  * 


***************************************************** 

* 

EXT  -ELIBR  ,*LIBX,WAIT 
ENT  DTiSB 

* 

COM  DELAY, DATA<8192) , BUFL 

* 

DTiSB  NOP 


JSB 

•ELIBR 

* 

TURN  OFF  INTERRUPTS  AND  MEMORY  PROTECT 

NOP 

* 

LDA 

DELAY 

* 

SB  A 

* 

SKIP  IF  DELAY  >  0  OR  =0 

CMA 

* 

COMPLEMENT  IF  DELAY <  0 

STA 

DELAY 

* 

LDX 

BUFL 

* 

SET  COUNTER 

* 

LOOP 

LDB 

STORE 

* 

SET  BUFFER  POINTER 

STC 

i  SB ,  C 

* 

SEND  ENCODE  SIGNAL 

SFS 

iSB 

* 

IS  DEVICE  BUSY? 

JMP 

*-i 

* 

YES.  TRY  AGAIN. 

L I A 

iSB 

* 

NO.  READ  DATA. 

STA 

STORE, I 

* 

STORE  DATA  IN  BUFFER 

INB 

* 

STB 

STORE 

* 

SET  POINTER  TO  NEXT  BUFFER  LOCATION 

JSB 

WAIT 

* 

DELAY  PERIOD  BEFORE  GETTING  NEW  DATA 

DSX 

* 

DECREMENT  COUNT.  SKIP  AFTER  LAST  DATA  i 

* 

JMP 

LOOP 

* 

GET  NEXT  DATA 

<T* 

CLC 

iSB  ,C 

* 

CLEAR  I/O  CONTROL  BIT 

JSB 

SLIBX 

* 

TURN  ON  INTERRUPTS  AND  MEMORY  PROTECT 

£ 

DEF 

DTiSB 

* 

RETURN  TO  CALLING  PROGRAM 

STORE 

DEF 

DATA 

* 

DEFINE  DATA  BUFFER  ADDRESS 

END 

DTiSB 

122 


NAM  WAIT! 


ASMB.L 

* 

* 


¥ 


%  ¥ 

*  THIS  SUBROUTINE  PROVIDES  A  WAITING  LOOP.  THE  * 

¥  VARIABLE  "DELAY"  DETERMINES  THE  LENGTN  OF  ¥ 

¥  TIME  WAITING.  IT  IS  A  NUMBER  FROM  0  TO  ¥ 

¥  32767 .  ¥ 

¥  ¥ 


¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥ 

* 


ENT 

WAIT 

COM 

DELAY 

WAIT 

NOP 

LDY 

DELAY 

¥ 

SET  TIMER 

ISY 

¥ 

INCREMENT  TIMER 

ONCE  BEFORE  DECREMENTING 

¥ 

NOP 

LOOP 

DSY 

¥ 

DECREMENT  COUNT 

.  TIMER=Q? 

JMP 

¥-i 

¥ 

NO.  CONTINUE 

* 

JMP 

WAIT, I 

¥ 

YES.  RETURN  TO 

CALLING  ROUTINE. 

END 

WAIT 
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(fin  rbi 4 


asms, l 


*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦< 


►♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 


♦  THIS  SUBROUTINE  IS  VERY  SIMILAR  TO  AD14E. 

♦  ITS  PURPOSE  IS  TO  PROVIDE  AH  APPROXIMATION 

♦  OF  THE  CONVERSION  TIME  PLUS  THE  TIME  REQUIRED 

♦  TO  STORE  DATA  WHICH  CCUPRS  IN  AD14B.  AD14 

♦  GOES  INTO  AN  INFINITE  LOOP  WHICH  CAN  EE  TERM- 

♦  INATED  BY  PRESSING  THE  "CLEAR  DISPLAY" 

♦  BUTTON  ON  THE  FRONT  PANEL.  WHILE  THIS  ROU- 

♦  TINE  IS  RUNNING?  THE  "STATUS"  SIGNAL  ON  THE 

♦  A -D  CONVERTER  CAN  EE  MONITORED  WITH  AN  OSCIL- 

♦  LOSCOPE. 

♦ 

♦  VARIABLE: 

♦  DELAY  -  ADDS  DELAY  TIME  BETWEEN  THE  ACSUI- 

♦  SITIGN  OF  EACH  TWO  DATA  POINTS 

♦ 

♦ 

EXT  EL  IE®? ELI BX- WAIT 
ENT  DT14E 


DT14E  MOP 
JSB 


SLIER 


♦  TURN  0®®  INTtPcUcTS  AND  MEMORY  ®-:0TEST 


.DA  DELAY 


DELAY 

EUFL 

STORE 


*  SET  THE  SWITCH 

*  REGISTER  TO  1 

* 

♦  SKIP  IF  DELAY'-  0  OF  =0 

♦  COMPLEMENT  IF  DELAY <0 

♦ 

♦  SET  COUNTER 


LOO®  STC 

SFS 


1 4  B .  C 
14B 
♦-1 
14B 
SAVE 
STOc'E 
WA I T 
1 


♦  SEND  ENCODE  SIGNAL 

♦  IS  DEVICE  BUSY' 

♦  YES.  TRY  AGAIN. 

♦  NO.  ®EAD  DATA. 


♦  DELAY  PERIOD  BEFORE  GETTING  NEW  DATA 

♦  LOAD  REG.  E  WITH  SWITCH  ®ES. 

♦  SKIP  WHEN  SWITCH  ®E5.  IS  CLEAR 

♦  GET  NEXT  DATA 


14E.C  ♦  CLEAF  I  ''0  CONTROL  S' IT 

SLIEX  ♦  TUCN  ON  INTERRUPTS  AND  MEMORY  ®®OTECT 

DT14E'  ♦  RETURN  TO  CALLING  PROGRAM 


124 


store  equ  o 

CHVE  EQU  0 
pijPL  EQU  0 


♦  DU  MM'.'  ‘■■’ARI.ABLE 

♦  DUMMY  VARIABLE 

♦  DUMMY  VARIABLE 


END  DT14B 


ASMB.L 


MAM  ADIS 


THIS  SUBROUTINE  IS  VERY'  SIMILAR  TO  AD1SB. 

ITS  PURPOSE  IS  TO  PROVIDE  AM  APPRO  I MAT ISM 
OF  THE  CONVERSION  TIME  PLUS  THE  TIME  REQUIRED 
TO  STORE  DATA  MHICH  OCUcFS  If!  ADI5E.  AD !  T 
GOES  INTO  AM  INFINITE  LOOP  MHICH  7AM  EE  TERM¬ 
INATED  BY  FREISING  THE  "CLEAR  DIIPLA. 

BUTTON  ON  THE  FRONT  PANEL.  .i^ILE  r-II  ROU¬ 
TINE  IS  RUNNING*  THE  "STATUS"  SIGNAL  ?  I-E 
A-'D  CONVERTER  CAN  BE  MCNITEC'ED  WITH  - 
LOSCDPE. 

j  ABLE  < 

DELAY  -  ADDS  DELA .  TIME  BETWEEN  THE  -ICUI- 
SITION  OR  EACH  TWO  DATA  =0:,(': 


DTI  SB 


ENT 

ILIEF ■ 

'■£LIE::> 

WA I T 

ENT 

DTI  SB 

COM 

DELAY 

NCR 

JSE 

ILIEP 

♦ 

TURN  OFF  INTERRUPTS  AND 

NOP 

♦ 

LDA 

1 

» 

SET  THE  SWITCH 

OTA 

I 

♦ 

REGISTER  TO  1 

LDA 

DELAY 

♦ 

SIR 

♦ 

SK I R  IF  DELAY  0  OR  =0 

CM  A 

♦ 

COMPLEMENT  IF  DELAY  0 

STA 

DELAY 

♦ 

LDX 

EUFl 

♦ 

SET  COUNTER 

LDE 

STORE 

♦ 

STC 

1SE-C 

♦ 

SEND  ENCODE  SIGNAL 

SFS 

1SB 

♦ 

IS  DEVICE  BUSY? 

JMP 

♦-1 

♦ 

YES.  TRY  AGAIN. 

LI  A 

l  SB 

♦ 

NO.  READ  DATA. 

STA 

SAVE 

♦ 

STB 

STORE 

♦ 

JSE 

WAIT 

♦ 

DELAY  C'EC'  IOD  EEFO'E  GET' 
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LIB  1  ♦  LOAD  - E 3 .  E  ’.II TH  SWITCH  REG. 

SIB  *  SKIP  WHEN  SWITCH  PEG.  IC  CLEAR 

JMP  LOOP  ♦  GET  MEYT  DATA 

• 

CLC  15B-C  ♦  CLEAR  I  CONTROL  SI” 

JC3  BL I B*<  ♦  TURN  ON  INTERRUPTS  AND  MEMORY  F'PC 

DEF  DT15E  *  RETURN  TO  CALLING  PROGRAM 

♦ 

STORE  EQU  0  ♦  DUMMY  VARIABLE 

SAVE  ECU  0  ♦  DUMMY  VARIABLE 

BURL  EQU  0  ♦  DUMMY  VARIABLE 

♦ 

END  DT15E 
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FTN4  •  L 


PROGRAM  AND62 


THIS  PROGRAM  TESTS  THE  SUBROUTINE  DMA  IN.  THE 
ENTRY  POINTS  APE  CHKGB,  BT14B,  AND  DT15B. 


INTEGER  DATA, BUFL, DEC  I DE • FNSH6 
COMMON  BATA  C3192)  *  EUFL,  FNSH6 

C -  INITIALIZE  VRPIRBLES 

Q 

MR  I TE  a.  2  00) 

"EAD  '•  1 ,  ♦)  BUFL 
BUFL  =  -BUFL 
MR  I TE  a*  BOO) 

READ  Cl ,  ♦  >  L'JMUM 
MRITE  Cl  ,  4 GO'' 

READ  Cl,*"-  DECIDE 

C -  CHECK  IF  DC PC  CHAM.  1  IS  AVAILABLE 

CN  CH6  =  0 
100  CALL  CHK6E 

IFCPNCHE.EO.0-.  GO  TO  100 

C -  GET  THE  DATA 

Ic  '-DECIDE.  EC*.  14)  CALL  DT14B 
I F  CDEC I  DE .  EC .  1 5'>  CALL  DT 1 5B 
C 

C -  MR IT  UNTIL  DRTR  RC CU IS  I T I ON  IS  COMPLETE 

c 

-NSHC  =  0 
150  CRLL  CHKiE 

I F  < FT  1 SH6 .  EG’ .  u)  bO  TO  ICO 

C -  OUTPUT  DATA 

BUFL  =  -BUFL 

MPITE  CLUNUM,5O0>  Cl • DRTR  CD , 1  =  1. BUFL) 


FORMAT  "HOM  MR  NY  DRTR  POINTS' 


MB  19c.  MR;  I) 


FOPMRT  i'/s  "MHRT  IS  THE  LU«  OF  THE  OUTPUT  DEVICES 
FOPMRT  </s*  "I.IHICH  D'Q  CHANNEL,  14  OP  15*“) 


FOPMRT  C22  C  ,  123  C4  1  217)  «  ' 
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stop 

END 


ACME?  L 


MAM  DMA I N 


♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 
♦  ♦ 

♦  THIS  SUBROUTINE  ACQUIRES  DATA  FROM  AN  A/D  CONVERTER  ♦ 

♦  USING  DIRECT  MEMORY  ACCESS.  ONE  OF  TWO  CHANNELS  ♦ 

♦  CAN  BE  SELECTED:  14  OR  15.  ♦ 


►♦♦♦♦♦♦ 


EXT 

'SLIER?  SLIBX 

ENT 

CHK6E?  DT14E 

?  DT15E 

COM 

DATA  <3 192::' 

« 1  , 

BUFL > FNSH6 

CHI 'EE  NOP 

♦ 

ENTRY  POINT 

JSE 

'SLIER 

♦ 

TURN  OFF  INTERRUPTS  AND  MEMORY  PROTECT 

NOP 

♦ 

CCS 

BE 

♦ 

IS  DCPC  CHAN.  1  AVAILABLE' 

JMP 

ATM 

♦ 

NO.  RETURN. 

LDA 

=  E  1 

♦ 

YES. 

STA 

fn  :h€ 

♦ 

SET  FLAG  WORD  TO  1. 

='TN  JSE 

SLIBX 

♦ 

TURN  ON  INTERRUPTS  AND  MEMORY  PROTECT 

DEF 

CHVBE 

♦ 

RETURN  TO  CALLING  PROGRAM 

DT14E  NOP 

♦ 

ENTRY  POINT 

JSE 

■SLIER 

♦ 

NOP 

♦ 

LDA 

BURL 

♦ 

STA 

CW  3 

♦ 

SAVE  -BURL 

LDA 

CWB 

♦ 

SPECIFY  THAT 

ADA 

=  E 1  0  0  0  0  0 

♦ 

OPERATION  IS  INPUT 

STA 

CMS 

♦ 

LDA 

CW114 

♦ 

OTA 

■SB 

♦ 

SEND  CW 1  TO  DCPC  CHAN.  1 

CL'S 

SB 

♦ 

°REPARE  M  REm.  TO  RElEIVE  lWc 

LDA 

CMS 

♦ 

OTA 

BE 

♦ 

SEND  CWB  TO  DCPC  CHAN.  1 

STC 

BE 

♦ 

PREPARE  FOP  CWB 

LDA 

CWB 

♦ 

OTA 

BE 

♦ 

SEND  CWB  TO  DCPC  CHAN.  1 

STC 

6E-  C 

♦ 

ACTIVATE  DCPC  CHAN.  1 

STC 

1 4  E  •  C 

♦ 

START  A/D  CONVERSION 

CL'S 

6E«C 

♦ 

DISABLE  DCPC  CHAN.  1  INTERRUPT 

CL  C 

1 4  E  ?  C 

♦ 

DISABLE  I  0  CHAN.  14  INTERRUPT 

JSE 

SLIBX 

♦ 

DEF 

DT14E 

♦ 

RETURN  TO  CALL I MS  PROGRAM 

i  .i 
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DT15B  MOP 

JCB  iLIEP 
MOP 

LDfl  EUFL 
:TH  CM3 
LDfl  CW£ 

ADA  =£100000 
CTA  CW£ 


LDfl 

OTA 

CLC 

LDfl 

OTR 

STC 

LDfl 

OTR 

CTC 

ETC 

CLC 

CLC 

JCE 

DEF 


CM  115 
6B 
£B 
CW£ 

£E 

£E 

CW3 

££ 


&  B  * C 
15B«C 
-SB?  C 
15B?  C 
‘EL  I BX 
DT15B 


:u  i  1 4 

:vn  15 


qct 

OCT 

BEF 

B 


120014 
120015 
OAT  in 


♦ 


♦  ENTRY  BGIMT 

♦ 

♦ 


♦ 

♦ 

♦ 


♦ 

♦ 

♦ 


♦ 


♦ 

♦ 

♦ 


♦ 

♦  DISABLE  I  Q  CHAN.  15  INTERRUPT 

♦ 

♦ 


♦ 


♦ 

♦ 


END  DT14E 


ft  SMB  »  L 


NAM  DVRS  0 


THIS  IS  A  DUMMY  DRIVER.  ITS  ONLY  FUNCTION  IS  TO 
JUMP  TO  fl  USER  SUPPLIED  SUBROUTINE  WHICH  PERFORMS 
THE  NECESSARY  DRIVER  FUNCTIONS.  THIS  PROCEDURE  ALLOW  I 
TROUBLESHOOTING  OF  THE  DRIVER  AND  CHANGES  TO  EE  MADE 
WITHOUT  REGENERATING  THE  OPERATING  SYSTEM.  A  LARGE 
BUFFER  AREA  IS  INCLUDED  SO  THAT  THE  CORRECTED  VERSION 
CAN  EVENTUALLY  EE  COPIED  INTO  THE  SYSTEM  AREA. 


ENT  DPS  0 » 1 . 5  0 »  C .  5  0 


1.50  NOP 

USB  DR30-I 
JMP  1.50.  I 

♦ 

C.30  NOP 

J SB  DF50.  I 
JMP  C.30*  I 

* 

DR 50  DEF  TMR50 


TMP50  NOP 

JMP  TMP50.I 
ESS  SO OB 
END 

AS ME. L 

« 

NAM  DVR51 


♦  INITIATION  SECTION  ♦ 


♦  CONTINUATION  COMPLETION  SECTION  * 


♦  THE  USER'S  PROGRAM  MUST  CHANGE  TMP50 

♦  STARTING  ADDRESS  OF  HIS  DRIVER  SUBRQU 


THIS  IS  A  DUMMY  DRIVER.  ITS  ONLY  FUNCTION  IS  TO 
JUMP  TO  A  USER  SUPPLIED  SUBROUTINE  WHICH  PERFORMS 
THE  NECESSARY  DRIVER  FUNCTIONS.  THIS  PROCEDURE  ALLOWS 
TROUBLESHOOTING  OF  THE  DRIVER  AND  CHANGES  TO  BE  MADE 
WITHOUT  REGENERATING  THE  OPERATING  SYSTEM.  A  LARGE 
BUFFER  AREA  IS  INCLUDED  SO  THAT  THE  CORRECTED  VERSION 
CAN  EVENTUALLY  BE  COPIED  INTO  THE  SYSTEM  AREA. 


5 1  •  C  . 1 
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* 


1.51 

NOP 

USE 

JMP 

DR51. I 
1.51? I 

♦  INITIATION  SECTION  ♦ 

•3.51 

NOP 

JSB 

JMP 

DR5 1  ?  I 

C .  5 1  ?  I 

♦  CONTINUATION  COMPLETION 

SECTION  ♦ 

DP  51 

♦ 

DEF 

TMF31 

♦  THE  USER'S  PROGRAM  MUST 

♦  STARTING  ADDRESS  OF  HIS 

CHANGE  TMP31 
DRIVER  SUEROU 

TMF31 

NOP 

JMP  TMF3 1  ?  I 
BSS  3  0  OE 
END 


ASMS?  L 

♦ 

NAM  DVR 52 

♦ 

* 

♦  THIS  IS  A  DUMMY  DRIVER.  ITS  ONLY  FUNCTION  13  TO 

♦  JUMP  TO  A  USER  SUPPLIED  SUBROUTINE  WHICH  PERFORMS 

♦  THE  NECESSARY  DRIVER  FUNCTIONS.  THIS  PROCEDURE  ALLOW 

♦  TROUBLE SHOOT  IMG  OF  THE  DRIVER  AND  CHANGES  TO  EE  MADE 

♦  WITHOUT  REGENERATING  THE  OPERATING  SYSTEM.  A  LARGE 

♦  BUFFER  AREA  IS  INCLUDED  SO  THAT  THE  CORRECTED  VERSION 

♦  SAN  EVENTUALLY  EE  COPIED  INTO  THE  SYSTEM  AREA. 


♦ 

ENT  DR3E'«  I.52»C.52 


I. 33  NOP  ♦  INITIATION  SECTION  ♦ 

JS3  DR32j I 
JMP  1.52.  I 

♦ 

C.  32  NOP  ♦  CONTINUATION  COMPLETION  SECTION  ♦ 

JSB  DR52?  I 
JMP  C. 52 j I 

♦ 

DP.52  DEF  TMF'52  ♦  THE  USER'S  PROGRAM  MUST  CHANGE  TMP'53 

♦  ♦  STARTING  ADDRESS  OF  HIS  DRIVER  SUEROU 

♦ 

TMP52  MOP 

JMR  TMP52? I 
ESS  3 DOE 
END 
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Appendix  F:  System  User's  Guide 

This  appendix  contains  supplementary  information 
which  should  assist  the  user  in  getting  acquainted  with 
the  DSP  laboratory.  It  is  not  intended  to  be  comprehen¬ 
sive  because  the  user  can  find  most  information  in 
appropriate  manuals  or  other  documentation. 

1.  The  system  master  security  code  is  AF. 

2.  Off-line  loading  procedures 

a.  System  Disk 

First,  make  sure  the  computer  has  halted,  and 
the  "DRIVE  READY"  light  on  the  disk  drive  is  on. 

Set  the  S  Reg.  to  111202  (the  disk  loader  ROM 
occupies  optional  loader  ROM  socket  #2;  Ref  14:7). 
Press  the  keys  PRESET,  IBL,  and  RUN  in  sequence. 
This  will  load  the  operating  system  into  memory. 

The  system  will  come  up  with  FMGR  running. 

b.  HP  Cassettes 

Select  the  left  or  right  cassette  drive  by 
pressing  GOLD,  one  of  the  FROM:  keys,  and  TRANSMIT 
in  sequence.  Then,  follow  the  same  procedure  as 
described  for  the  disk,  but  set  the  S  Reg.  to 
041300  (optional  loader  ROM  #1). 

c.  Paper  Tape  Reader 

Follow  the  same  procedure  as  described  for  the 
disk,  but  set  the  S  Reg.  to  001700. 

3.  Compile  and  run  a  FORTRAN  IV  program 

To  compile  a  program,  enter  the  following  command: 
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*RU , FTN4 , &f name , 1 

The  asterisk  is  the  RTE  prompt.  The  source  code 
file  name  is  &fname.  HP  convention  is  to  start  all 
source  code  files  with  a  and  all  object  files 

with  a  For  the  above  command,  if  no  file 

exists,  called  %fname,then  RTE  will  create  the  file, 
stop  the  compilation,  and  report  a  DISASTER  97.  In 
this  case,  just  type  the  same  command  in  again,  and 
it  should  run.  If  %fname  existed  before  the  command 
was  entered,  the  new  object  code  will  overwrite 
it  as  the  compilation  proceeds. 

To  load  and  run  the  program,  enter  the  following 

commands  (the  colon  is  the  FMGR  prompt): 

*RU , FMGR 
:LG,0 

: MR , %f name 
: RU , LOADR ,99,1 

If  no  errors  were  discovered  by  the  loader,  it  will 
report  that  the  program  is  ready.  To  run  the 
program,  type: 

: RU , pname 

where  pname  is  the  name  of  the  program.  Usually  the 
FMGR  prompt  will  return  before  all  of  the  output 
has  been  sent  to  the  designated  device.  To  see  the 
remaining  output,  strike  the  RETURN  key.  This  will 
also  terminate  FMGR. 

A  program  can  also  be  run  without  FMGR  so  that  no 
interruption  in  the  output  occurs.  Just  enter  the 
EX  command  after  running  the  loader.  Then  type: 


*RU , pname 

4.  Assemble  and  run  an  assembly  language  program. 

Enter  the  following  commands: 

*RU , FMGR 
:  LS ,  0 

:MS , &f name 
:  LG ,  0 
:  LG ,  3 
:  EX 

*RU , ASMB ,2,1,99 
*RU , FMGR 
:SA,LG, Ifname 

The  program  can  now  be  loaded  and  run  in  the  same 
way  as  described  for  FORTRAN  IV,  except  that  the 
MR  command  is  not  required. 

5.  Report  Files 

The  Source  code  file  name  of  any  program  listed  in 
this  report  is  just  the  program  name  preceeded  by 
a  The  corresponding  object  code  files  are 

preceeded  by  a  Any  corresponding  type  6  file 

has  the  same  name  as  the  program  name. 

6.  TI  Terminal 

The  TI  terminal  is  currently  being  operated  at  110 
baud.  A  300  baud  interface  is  available,  but  the 
terminal  has  trouble  completing  a  carriage  return 
fast  enough.  Another  printing  terminal  has  been 
successfully  operated  on  the  same  I/O  channel  using 
the  300  baud  interface. 
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Appendix  G:  Disk  Storage  Programs 

The  two  programs  implemented  for  storing  data  on 
the  disk  are  listed  in  this  Appendix.  STDSK  demonstrates 
temporary  disk  storage.  FILES  demonstrates  permanent 
disk  storage  in  a  disk  file.  Both  programs  are  self- 
documenting  and  require  no  further  explanation. 


FTN4. L 


PROGRAM  STDSt 


C 

c 


THIS  PROGRAM  USE'  A  DISK  TRACK  ALLOCATION  EXEC  CALL  TO 
ALLOCATE  TRACKS  FOR  DATA  STORAGE.  THE  SUBROUTINE  BINRV1 
STORES  DATA  ON  THE  DISK. 


VARIABLE 

Q 

IBUFR 

r 

Q 

IBUFL 

P 

IDISK 

p 

IS  TRK 

r 

ISC  TR 

r 

IOFCT 

IRDWR  - 

r 

I CODE  - 

THIS  IS  THE  DATA  BUFFER  TO  BE  FEAD  FROM  OR 
WRITTEN  TO 

LENGTH  OF  THE  BUFFER 

LOGICAL  UNIT  NUMBER  OF  THE  DISK  USED  IN  THE 
DATA  TRANSFER 

DISK  STARTING  TRACK  NUMBER 
DISK  STARTING  SECTOR  NUMBER 

D I  SI  SECTOR  OFFSET  • I.E. •  THE  NUMBER  OF  WORDS  TO 

TO  BE  SKIPPED  IN  THE  SECTOR 

DETERMINES  READ  OR  WRITE  OPERATION 

1  =  READ  2  =  WRITE 

SPECIFIES  THE  TYRE  Oc  EXEC  CALL 


DIMENSION  IBUFR<lflO> 


allocate  tracks 


"R  I TE  '•!  1  j  j  U ■' 

?  0  FORMAT  .  "  HOW  MANY  TRACKS?  "  > 

=EAD'1<»>  ITRAK 
IS ODE  =  4 

call  EXEC  <I  CODE*  ITRAK?  ISTRKj  IDISK  «  I  SECT"* 


C -  INFORM  USER  OF  DATA  STORAGE  AREA?  0=  THAT  IT  IS 

I -  NOT  AVAILABLE 


I p  a S TRK . EO . - 1 >  CALL  EX I T 
WRITE  <1  ?  'j-1  ISTRK,  IDISK?  I  SECT 

FOFMAT < " THE  STARTING  TRACK  NUMBER  IS" ? 14.  - 
♦"THE  DISK  LU  NUMBER  IS"?I4«^? 

♦  “THE  NUMBER  OP  G4  WORD  SECTORS'-  TRACK  IS"?  14.  v.> 


C -  SET  VARIABLES  FOR  EINRY 

DO  tOO  1=1? 100 

100  IB:'JPR  a.:-  =  I 

IBUFL  =  100 
I CCTP  =  0 
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iofst  =  o 

IRDMR  =  2 

C -  MR  I TE  THE  BUFFER  VALUES  TO  THE  D I "R 

r 

CALL  BNP  VI  I BUFP «  IBUFLj  IDISK.  ISTRr  .  ISCTF.  I  OF 'ST  *  IP  BMP  1 

C -  SET  THE  BUFFEP  VALUES  TO  OEFO 

DO  200  1=1 « 100 
200  IEUFR •!>  =  0 

C -  READ  THE  BUFFER  VALUES  BACK  IN  FROM  THE  DISK 

C 

I FOUR  =  1 

CALL  BNP VI  <1  BUFP »  IBUFLj  IDIS5-  •  I  STRf  .  ISCTR.  IOFST.  IRDMR;- 

S -  MR  I  TE  OUT  THE  BUFFER  VALUES  TO  VERIFY  THE  C'EAB  UP  I  TE 

S -  OPERATIONS 

MR  I  TE  ■  1  <  :■  00 1  '  I  BUFP  -  I  '•  «  I  =  1  •  1  0  0  • 

3  0  0  FORMA T  •: "  BUFeEP  VALUE  S  ”  .• ,  .  1  0  >  1015*  >  ■ 

STOP 

END 


SUBROUTINE  BURY  l  Cl  BUFP >  IB'JFL*  IDISK.  I  ST®K  .  ISCTR.  IOFST*  I'D-MP- 


THIS  SOUTINE  USES  THE  SYSTEM  UTILITY  “ B I MC'Y "  rOP  READING 
PROM  AMD  MR  ITT  I MG  TO  A  DISK.  THE  ENTRY  RQINT  -UP  'EASING 
IS  "E'EAD".  THE  ENTRY  RQINT  FQR  UR ITT I MG  IS  “BUR  IT”. 


D I MENS 1  ON  I  BUFF:  <  I EUFL 


IF  IRDMP.  EG.  I>  GO  TO  100 
IP  - IFDMF.EC.2--  GO  TO  200 
MR  I  TE  -.1  j  5  0> 

SO  FORMAT <"♦♦♦♦♦  ERROR  IN  INPUT.  IP BUR  DOES  NOT  EQUAL  "• 
+■  "1  OF'  2  ♦♦♦♦♦"> 

RETURN 

100  CALL  BREAD  <  IEUFR.  IBUFLj  IDISK.  ISTP'K  >  ISCTR.  iaF:T> 

RETURN 

200  CALL  BMP' IT -'IBUFR.  IBUFLj  IDISK-  ISTPK-  ISCTR.  IOcST- 
RETURN 
END 
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ZJBROUTIME  EXIT 


HIE  TUB  ROUT  I  ME  II  CALLED  WHEN  THE 
UNAVAILABLE.  THE  USER  IZ  IMFCRMED 
IZ  ZTOF'PEB. 


TRACK! 
OF  THIZ 


ECUS Z TED  AF 
AMD  THE  FRO 


'..IP I TE  1 «  1  000  1 

1000  FORMAT  >'■'.■■  ■<  "TRACK <!>  APE  NOT  AVAILABLE.  "  * 
+  " PROGRAM  TERMINATED.  "  ? 

ZTOP 

END 
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<ri  pi 


:•  i  ».» 


cTN4«  L 


PROGRAM  FILES 


C* 

£♦ 

C* 


THIS  PROGRAM  DEMONSTRATES  THE  USE  OF  FMP  CALLS  FOR  ♦ 
DATA  STORAGE  ON  A  DISK.  A  FILE  IS  CREATED?  AND  DATA  ♦ 
IS  WRITTEN  INTO  AND  READ  FROM  THE  FILE.  ♦ 


INTEGER  DCB1  <144>  ?  NAME  1  <3*>  «SIZE1  <2>  ? CR? TYPE 1 « SECUP < 

+  ERROR  ?  OPTN  ?  ARAY 1  <  1  0>  >  ARAY2  <  1  0>  ?  A  PAY  3  ■'  1  0  > 

i“- 

NAME1  <1.>  -  0 
NAME1  <:■£>  =  0 
NAME  1  CT>  =  0 
UP  ITS  a  ,  15"' 

15  FORMAT  -VENTER  FILE  NAME") 

READ  <1,5 0>  <NAME  1  <  D  ,  1*1,  3> 

50  FORMAT  C3A£> 

3IIEl<i:>  =  3 
TVFE1  =  3: 

3 EC UR  *  0 
CP  -  =0 
C 

: -  CREATE  .A  TYRE  3  FILE  ON  C ARTRIDGE  10 

CALL  CREAT  -'DCBl «  ERROR,  MAMS  1 ,  C IZE1 «  TYPE  1 ,  SEC  UR ,  CR> 

I  =•: ERROP.LT.  O')  30  TO  3 0 0 

'■JR I TE  a  ,  1  0 0.'  (NAME 1  <I>  ,  1*1,  3>  *  CR .  ERF GR 
100  FORMAT  V FILE  '  “  «  3A1,  "  '  MAC  BEEN  CREATED  ZN  CARTRIDGE 
+  II."  MITM  ",13."  E4-U0RD  CEC Tu=C  "V"' 

: -  O-'EN  THE  FILE  C0R  READ  •'VP  I  TE 

OFTM  =  0 

CALL  OREM  <DC£1 ,  ERROR,  MAME1 .  ORTH,  CECUR,  CR"' 

IF (ERROR. LT. 0>  50  TO  900 

C -  GENERATE  DATA  INTO  THREE  ARRAYS 

DO  100  1=1,10 
ARAY1 (1)  =  I 
A  RAYS  ''  I )  =  !♦  I 
AF.AV.3a:'  =  1*1 
00  CONTINUE 


C -  !.IR I  TE  EACH  DATA  ARRAY  AC  ONE  RECORD  IN  THE  FILE 

I**' 

LEN5H  =  10 
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.IP  I TF  •  DC  El -  ECRG=  -  ARAY  1  •  LEM  3H'. 

IP  .  £PPCP.  LT.  0)  GG  TG  •? ij 0 

IA._L  '.IP  ITF  (TCE1  ?  E=FOR-  hRAYLH  LEHIGH ' 

I F  •  ERROR .  L T .  O'-  EG  TO  A n 0 
C  ALL  ‘.'IP  I  TF  *  DCS  1  *  EPF'GR *  ARAY3 .  LEHGH  1 
I f  i  ERROR .  LT .  0.>  GG  TG  A00 

c —  :et  data  arrays  tg  gerg 

DG  300  1=1-10 
ARRV1  CI.>  =  0 
A RAYS '  I  *  —  0 
RPR  Y  3  'I  I  ■>  —  0 
300  COAT  I  HUE  v 

3  -  READ  EACH  °ECGRD  FROM  FILE  INTO  AH  ARRAY 

CALL  fmhDF  CDCB  1  •  ERROR > 

I F  t  ES  =  OR .  L  T .  0"'  GO  TG  A  0  0 

-ALL  P'EADF  CDCE-1  -  cr.R'OR' -  ARAY1  -  LEHGH"1 

IF  ('ERROR.  t_T.  0>  30  TG  A 0 0 

-ALL  P  EADF 1 DCE 1  -  ERROR  - RP  AV2 > LEH'GH > 

IFCERRGR.LT. 0>  GO  TG  A  00 

CALL  RERDP'  DCEl  -  ERROR-  ARAY3.-  LEHGH ;• 

IP  ERROR. LT. 0>  GO  TG  A 00 

I -  -IP I TE  THE  DATA  ARRAY 3  TG  VERIFY  CEAD  WRITE  GFERATIG 

i.ilr  I  TE  Cl  *  4  0  O'-"'  CAPAV 1  C I  *  -  I  =  1 «  1  0  ■'  .  C  AR  ACE'  Cl'*  •  I  =  1 .  1  0  1  - 
+  CAR  AY  3  C I .)  -  1  =  1-  1  0  •* 

4  0  0  FORMAT  c  3  <  1  ij  1 7  > 

-  CLOSE  DATA  PILE 

CALL  C LG I E CDCE 1 , Ec POP > 

Ic  ERROR. LT. 0)  GO  TG  A00 

:tgr 

G  0 0  .IR I  TE  C 1  ?  1  U  0  U  '  ERR  GR: 

1  U  U  0  PQRMAT  1 .  ’  ”  PMF'  ER'RGF'  "-  14  1 


LI* 


■3SS35S 


Appendix 
This 
APLOT.  It 


!:  Auto  Plot  Program 

Appendix  contains  a  listing  of  the  program 
is  self-documenting. 
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f  TN4 , L 

C 

C 

PROGRAM  APLOT 
C 

r'  |L«  ^  ^  4  |«  4 L*  tl,  4 Li  4  j.  ^  sb  ^  4 Li  4 Li  4 Li  ^  4 Li  .  Li  4 Li  ^  ■  Li  \U  4 Li  ^  tji  .  Li  4 Li  ^  s U  4 Li  Oj  4 Li  Uj  (L  iLi  *  Li  -  u  ■!.  4 Li  Ui  \L  4 Li  Ui  \L  •  Li  4 Li  Li  it.  .  Li  4 Li  4 Li  4 Li 


c*  * 

C*  THIS  PROGRAM  DEMONSTRATES  COMPUTER  CONTROL  OF  THE  * 
C*  GRAPHICS  TERMINAL.  THE  EXECUTIVE  CALL  IS  USED  TO  * 
C*  SEND  THE  DECIMAL  VALUE  OF  THE  ASCII  CHARACTERS.  * 
C*  THE  SEQUENCE  OF  CHARACTERS  LISTED  IN  BUFR  WILL  * 
C*  FILL  THE  AUTO  PLOT  MENU ,  ACCEPT  DATA  FROM  THE  COM-  * 
C*  PUTER  AND  PLOT  IT,  AND  TURN  OFF  AUTO  PLOT.  * 
C*  # 


C*  VARIABLES: 

C*  CODE  -  DESIGNATES  A  WRITE  EXEC  CALL 

C*  CTLWD  -  SPECIFIES  THE  SYSTEM  CONSULE  LU* 

C*  BUFR  -  CHARACTER  BUFFER 

C*  BUFL  -  BUFFER  LENGTH ;  =  -2*(#  OF  WORDS  IN  BUFR) 

C* 

^  ^  ^  if  tf  <7k  'T  «r  'F  <r  'p  t  »r  ^  'i'  'A  t  ^  ^  'r  'f  'p  '(*  'r  'i'  'i'  'f  'f  <p  'f  'n  'J'  *1'  'f  '7^  'i'  'l' t  'p 

C 

INTEGER  CODE, CTLWD ,BUFR (43)  ,  BUFL 
C 

C -  INITIALIZE  VARIABLES 

C 

CODE  =  2 
CTLWD  =  1 
BUFL  =  -S6 
BUFR(i)  =  27 
BUFR ( 2 )  =  42 
BUFR ( 3 )  =  97 
BUFR (4)  =  100 
BUFR (5)  =  SO 
BUFR ( 6 )  =  104 
BUFR ( 7 )  =  49 
BUFR ( 8  )  =  105 
BUFR ( 9 )  =  50 
BUFR (10)  =  106 
BUFR (11)  =  49 
BUFR ( 12 )  =  107 
BUFR (13)  =  48 
BUFR (14)  =  103 
BUFR (IS)  =  49 
BUFR (16)  =  48 
BUFR ( 1 7 )  =  48 
BUFR ( 18 )  =  109 
BUFR ( 19 )  =  48 
BUFR ( 20 )  =  110 
BUFR ( 21 )  =  49 
BUFR (22)  =  43 


*  a  a  a  *  •* 


o  u  u  o  u  u  *-t  rj  u  a  u 


BUFR  <  23 )  =  48 
BUFR  ( 24 )  =  ill 
BUFR  < 25 )  =  50 
BUFR  ( 26 )  =  48 
BUFR  <27 )  =  112 
BUFR  ( 28 )  =  49 
BUFR (29)  =  48 
BUFR ( 30  )  =  113 
BUFR (31 )  =  50 
BUFR ( 32 )  =  48 
BUFR ( 33 )  =  114 
BUFR (34)  =  49 
BUFR (35)  =  48 
BUFR ( 36 )  =  115 
BUFR ( 37 )  =  49 
BUFR ( 38 )  =  52 
BUFR (39 )  =  117 
BUFR (40)  =  49 
BUFR ( 41 )  =  118 
BUFR (42 )  =  99 
BUFR ( 43 )  =  65 

-  SET  UP  AUTO  PLOT  MENU  AND  DRAW  GRID 

CALL  EXEC ( CODE , CTLWD , BUFR , BUFL ) 

—  SEND  DATA  TO  BE  PLOTTED 

J  =  1 
K  =  0 

DO  200  1=1,92,7 
UR ITE (1,100)  I,J 
00  FORMAT ( 217 ) 

K  =  K  +  i 
J  =  J+K 
00  CONTINUE 

-  TURN  OFF  AUTO  PLOT 

BUFL  =  -8 
BUFR (4)  =  66 

CALL  EXEC(CODE, CTLWD, BUFR , BUFL) 

STOP 

END 
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Appendix  I:  D/A  Converter  Software  Listings 

The  program  TSTDA  generates  data  for  subroutine 
DA16B  and  calls  it.  DA163  provides  software  support  for 
the  D/A  converter  which  has  an  X-Y  plotter  and  a  conven¬ 
tional  oscilloscope  connected  in  parallel  to  its  outputs. 
Subroutines  MINXY,  MAXXY ,  SCALE,  CBINE ,  PLOTR,  and  SCOPE 
are  called  by  DA16B.  Subroutine  OSCPE  (entry  point  SCPE) 
is  called  by  SCOPE.  Subroutine  PLT  (entry  point  PLOT) 
and  ONE  (entry  point  FIRST)  are  called  by  PLOTR.  The 
purpose  of  each  subroutine  is  explained  in  the  program 
listing . 
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FTN4, L 


program  tstda 

r 

C  ♦  ♦ 

THIS  PROGRAM  GENERATES  DATA  IN  SAMPLE  POINT  FORM  ♦ 

C*  FOP  A  CONTINUOUS  FUNCTION.  THE  DATA  IS  USED  TO  ♦ 

O  TO  DEMONSTRATE  THE  SUBROUTINE  'DA16E  .  ♦ 

♦ 

C ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦ 

COMMON  NPT  S .  I  XAXS  < 356  >  ,  JYAXS  <£56>  » XAX I S  <256>  »  VAX  I S  <236 

•: -  INPUT  THE  NUMBER  OF  DATA  POINTS 

WRITE  a«30> 

50  FORMAT  <22  <  ">  ?  "HGU  MANY  DATA  POINTS?  <236  MAX'”;' 

READ <1 s *>  NCTS 

S -  GENERATE  THE  DATA 

PI  =  4.  O^ATAN  <1 . 0) 

STEP  =  5.  'j*P  I 'FLOAT 'iMPTS  > 
r 

DO  100  IM.NPTS 

SAVE  =  FLOAT' I':' ♦STEP 
••.AXIS';' I)  »  SAVE 

VAX  I  S  <  I  ;•  «  SIN  <SAVE>  +  1  1 . O-'SAVE)  ♦CCS  SAVE) 

100  CONTINUE 
r 

C -  SEND  DATA  TO  D'-A  AND  OUTPUT  DEVICE 

CALL  DAI SB 


SUBROUTINE  DAI  BE 

?♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•♦♦♦♦♦♦ 
'♦  ♦ 

:♦  THIS  ROUTINE  FORMATS  DATA  rOR  OUTPUT  TO  EITHER'  AN  ♦ 

C*  X-Y  PLOTTER  OP  A  CONVENTIONAL  OSCILLOSCOPE.  THE  ♦ 

!♦  INPUTS  AC'E  TWO  REAL  ARRAYS  CORRESPONDING  TO  THE  i  ♦ 

.;♦  AND  V  AXES-  AND  THE  NUMBER  OF  DATA  POINTS.  ♦ 

:♦  ♦ 

c ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦• ♦♦♦♦♦♦ 


COMMON  Ne-TI »  IXAXS  <236  '■  •  JVAXS  •  256)  « XAXIS  ■  136"'  . '  AXIS':  256 ' 
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FIND  THE  MINIMUM  •••'ALUES  IN  :A  'IS  AND  CANIS 

call  minkvo:min»vwin.. 

SET  XMIN  AND  VMIN  TO  MAKE  Pill  ARRAY 
VALUES  POSITIVE 

I F  •".'  CM  I N .  ST .  0 .  0>  XM I N  =  0 .  0 
I F  <YM I N .  ST .0.0 >  YM IN  =  0.0 
XM  I N  =  ABS  OIIN) 

VMIN  =  ABSCVMIN> 

FIND  THE  MANSMUM  VALUES  IN  NANIS  AND  YANIS 
CALL  MANNY CNMAN * YMAN  > 

SCALE  ALL  ARRAY  VALUES 
C  ALL  SCALE  <.NM I N  •  YM  I N .  NMAN  •  YMA'  O 
COMBINE  INTESER  ARRAYS  INTO  ONE  ARRAY 

write <1 1  ioo>  cinaxsci')  »  jy*xs <i>  •  i»i.hpts> 

FORMAT  i  2'T  <1415*  ■'  >  » •  '  1 
CALL  ‘SEINE 

21L2CT  OUTPUT  DEVICE:  PLOTTER  OR  CSC I L LOSS 
‘.‘•RITE  '1  >  130'1 

cORMAT <xx»  ’’WHAT  TYRE  OF  OUTPUT  DEVISE*” •  - 
*  ”l=PLOTTEP  =*OSCILLO SCORE" 

READ  C  I  .  ♦  I=  I  SX 

OUTPUT  DATA 

Is  CIPICK.ES.  1.'  CALL  PLOT® 

Ir  < IP  I C>  . EO . £ ■  CALL  SCORE 

RETURN 

END 


SUBROUTINE  M I  NX  Y  CXMIN.  VMIN'. 

THIS  ROUTINE  DETERMINES  THE  MINIMUM  VAsJJE  -OF  EACH  A 

COMMON  NPTS «  I  NANS  C£5<5>  *  JYAXS  CESS  >  » '  ’A’ ;  I S  ‘.EOS'1  .  YA  IIS  C£3 

XMIN  =  X AXIS  Cl  > 

VMIN  s  YANIS  Cl'* 

DO  £00  I=£«NPTS 

IF  (NANIS  •'!  '  .LT.  '.M IN-'  T’. IN  = 


'•.A  'IS'  I > 


I F  (YAX I S  a  :•  .  LT .  YM I  M  >  VI I M  =  YAX I  S  I ) 
200  CONTINUE 
C 

RETURN 

END 


SUBROUTINE  MAXXY  CXMAX.  YMAX'J 

THIS  ROUTINE  DETERMINES  THE  MAXIMUM  VALUE  FOP  EACH  AXIS 


COMMON  NPT  S » IXAXS  <£56> «  JYAXS <£36> » XAX I S  <£36> > VAX  I S <£36> 


XMAX  =  XAX  I S  ‘T  > 
YMAX  =  VAX  IS1' I  > 


DO  300  1=2* NPTS 

IF  fXAX  IS  CD  .  GT.XMAX.’i  XMAX  =  XAXIS'T> 
I F vYAX 1 1 < I > . GT . YMA  O  YMAX  =  VAX I  3  < I • 
300  CONTINUE 


RETURN 

END 


SUBROUTINE  SCALE  '  . '.MIN.  VNIN*  XMAX.  YMAX' 

COMMON  NFTSj  IXAXS  <i3G>  .  JYAXS  '  £5E  '  •  'AX I 

MAKE  ALL  ARRAY  VALUE'S  RANGE  FROM  0  TO 
INITIALIZE  THE  INTEGER  ARRAYS 


.  YAXIC'-£5G> 
ONE  EYTE 


DO  400  1=1* NETS 

XAX  IS  a>  *  c'X.MIN+XAXI  S  (.1  ■  '  CMIM+XMAX)  *355 
IXAXS  CD  »  I^IX-XA-I  S  1  I 

VAX  i  o  i  I  >  =  CYMIN+YAXISa'  '  ■ ','MIN+YMAX>  *£55 

JYAXS  <  I>  =  I  -  I X  <VAX  I S-  <  I  >  • 

400  SONTINUE 


RETURN 

END 


SUBROUTINE  SCOPE 
INTEGER  READY 

l“ 

MR  I TE  a  -  5  O 0> 

500  FORMAT  <  55  (1H*'  .  "HIT  A  CARRIAGE  RETURN  TO  START  THE 

♦  "SCORE  DISPLAY.  " « "C'R£SS.  CLEAR  DISPLAY  TO  TERMINATE 

*  '  i  55  (  1H*>  < 

READ 1 1 « ♦  '  READY 
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CALL  CCPE 


RETURN 

END 


SUBROUTINE  PLOTR 

C 

INTEGER  READY 

C 

C -  POSITION  THE  PLOTTER  PEN  BY  SENDING  THE  FIRST  POINT 

C' 

CALL  FIRST  ' 

WRITE  1 » 600> 

600  FORMAT  *  55  C 1 H-*  :•  «  " »  "  SET  THE  PLOTTER  FEN  DOWN.",--, 

+  "HIT  A  CARRIAGE  RETURN  WHEN  READY  TO  PLOT."- 

+  -  55  '■  1  H'O  . ''  ) 

R  EAD  <.  1  -  ♦  1  READY 

C 

CALL  ='LOT 

£ 

WR I TE  (.  1  ?  6 1  0  > 

610  FORMAT  (.//%  £5  <1 H* >  -  .  "PLOT  COMPLETE.  LIFT  DEN. 

+  25UH*)'^ 

RETURN 

END 


ACME  *  L 


♦ 

NAM  CMBNE 


♦ 

•  THIS  ROUTINE  COMBINES  TWO  INTEGER  ARRAYS  INTO 

*  JYAX!  IS  THE  MOST  SIGNIFICANT  BYTE-  AMD  V  AS! 

♦  THE  LEAST  SIGNIFICANT  BYTE. 


♦ 

» 


» 

* 


♦ 

ENT  CEINE 


COM 

NPTS « i ha: 

1.; 

£56>  -  JYA  ;!S  '"E56'-1 

SEINE  NOP 

♦ 

ENTRY  POINT 

CLA 

♦ 

CSX 

♦ 

CLEAR  THE  BUFFER  FTP. 

LOOP  LAX 

JYAHS 

♦ 

MOST  SIGNIFICANT  BYTE 

ALF- 

ALF 

♦ 

ROTATE  REG.  A  LEFT  ONE 

BYTE 

and 

=  E 1 7 7 4  0  0 

♦ 

MAKE  SURE  =  ISHT  BYTE  IS 

CLEAR 

LBH 

i  ■  :ahs 

♦ 

LEAST  SIGNIFICANT  BYTE 
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.  ..as / 


:te 

TEMP 

♦ 

pop 

TEMP 

♦ 

COMBINE  TME  TWO  BYTES 

IPX 

IXAXS 

♦ 

INTO  ONE  WORD 

r  *r  • ,,« 

1  j-  *. 

♦ 

SET  FTP.  TO  NEXT  BUFFER  LOCATION 

CXA 

♦ 

C»A 

NPTS 

♦ 

COMPARE  BUFFER  FTP.  TO  NPTS 

JMP 

CBINEj 

I 

RETURN  TO  CALLING  PROGRAM 

.IMP 

LOOP 

♦ 

CONTINUE 

TEMP 

1 

END 

CBINE 

ASMS? 

L 

S 

♦ 

♦ 

♦ 

NAM 

OSCF'E 

♦  ♦♦♦♦ 

♦ 

♦  ♦♦♦•« 

•  ♦♦♦♦♦♦ 

* 

♦ 

THIS 

:  POUT I 

NE  SENDS  DPTh  pP'OM  P  BUFFER  TO  PN  OSCIL-  ♦ 

♦ 

LOSS 

□  PE  THROUGH 

PDA  CONVERTER.  * 

♦ 

* 

♦♦♦♦■« 

•♦♦♦♦♦♦ 

♦ 

EXT 

JL I  E:P « 

■IL I 

ENT 

SC  PE 

♦ 

COM 

NF'TS j  I 

IflXS  <c 

;56.> 

■rroc 

NOP 

♦ 

ENTRY  °OINT 

JOB 

•ELIEF: 

♦ 

TURN  OFc  INTERRUPTS  AND  MEMORY  PROTECT 

NOP 

♦ 

!_DA 

=E  1 

♦ 

OTP 

IB 

♦ 

SET  THE  SWITCH  REGISTER  TO  1 

TP^CS 

CLP 

♦ 

SET 

CAX 

♦ 

THE 

LDY 

NPTS 

♦ 

COUNTERS 

CLP 

16E 

♦ 

START  THE  REFRESH  TIMER 

NE’-.T 

LfiX 

I XAXS 

♦ 

OTh 

l£E 

♦ 

SEND  DATA  POINT  TO  BE  c LOTTED 

ISX 

♦ 

INCREMENT  BUFFER:  RTF. 

DCY 

♦ 

DECREMENT  BUFFER  COUNTER 

JMP 

NEXT 

♦ 

SFS 

I6E 

♦ 

JMP 

♦-1 

♦ 

LIP 

IB 

♦ 

CHECK  THE  SWITCH  REGISTER 

OOP 

♦ 

IS  IT  CLEAR"'' 

JMP 

TRACE 

♦ 

NO.  OUTPUT  ELFFFR:  AGAIN. 

CLC 

1GB 

♦ 

YES'.  CLEAR  I  '□  CONTROL  BIT 

JOE 

IL  I BX 

♦ 

TURN  ON  INTERRUPTS  AND  MEMORY  PROTEST 

DEF 

SCCE 

♦ 

RETURN  TO  CALLING  PROGRAM 

♦ 
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END  SCPE 

RSMB?  L 

♦ 

NRM  PLT 

♦ 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 


♦  ♦ 

♦  THIS  ROUTINE  SEND  BRTR  FROM  R  BUFFER  TO  RN  X-Y  ♦ 

♦  PLOTTER  THROUGH  R  D-R  CONVERTER.  ♦ 

♦  ♦ 


♦  ♦ 


EXT 

■£LIEF. 

•Eli'S:: 

EMT 

PLOT 

COM 

NPTS  ? I 

xrxs  ; 

B5*> 

-LOT 

MOP 

♦ 

EMTPY  POINT 

JBB 

■ELIEP 

♦ 

T!JRM  OFF  IMTERFUPTC  PMD  r,1Ef,1  CP"/  PROTECT 

NOP 

♦ 

TFPCE 

CLP 

♦ 

BET 

CPX 

♦ 

T'PE 

LDY 

N»TB 

♦ 

'BOUNTERC 

♦ 

NEXT 

lpx 

ixpxb 

♦ 

LOPD  BUFFEP  VPLUE 

♦ 

LDE 

-B-327 

£7  ♦ 

BET  UP IT  TIMES 

MiilT 

uTP 

1GE 

♦ 

□TP 

16  E 

♦ 

IME- 

1  S«i.  £ 

♦ 

IMCPEMEMT  COUNT.  - IMIBHSB  MPITING* 

;mD 

UP  IT 

♦ 

MO.  CONTINUE. 

I B  X 

* 

YEB.  UFBPTE  BUFFEP  POINTE* 

EBY 

♦ 

DECPEMEMT  COUNT.  LPIT  BPTP* 

JMP 

MEXT 

♦ 

NO.  GET  NEXT  EPTP  ='CINT 

CLC 

liE 

♦ 

YEB.  CLEPF  I-O  CONTROL  BIT 

JBB 

EL  I  EX 

♦ 

TUPN  GM  IMTEFP'JPTB  ANB  MEMCPY  =CGTECT 

DEF 

DLOT 

♦ 

PETUPM  TO  CALLING  PROGRAM 

EM  It 

FLOT 

R:  MP.iL 

♦ 

♦ 

NRM  ONE 

♦ 


♦  ♦ 

♦  THIS  ROUTINE  SENDS  THE  FIRST  BUFFER  VRLUE  TO  RN  ♦ 

♦  X-V  PLOTTER  THROUGH  R  D-'R  CONVERTER.  THIS  MILL  ♦ 

♦  POSITION  THE  REN  ’0  TMF  •; ::EC  CRN  CET  IT  DOWN.  ♦ 

♦  ♦ 
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♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 


♦ 

EXT 

XL IBP' 

*£L  I BX 

A 

ENT 

FIRST 

♦ 

COM 

NPTSj I 

XflXSO 

L> 

cipst 

NOP 

♦ 

ENTRY  POINT 

USB 

XL  I  BP 

♦ 

TURN  OFF  INTERRUPTS  AND  MEMORY  PROTECT 

NOP 

♦ 

LDA 

IXAXS 

♦ 

LORD  FIRST  BUFFER  VALUE 

OTA 

1 6B 

♦ 

OUTPUT  BUFFER  VALUE 

CLC 

1GB 

♦ 

CLERP  I  -0  CONTROL  BIT 

JSE 

XLIBX 

N  ♦ 

TURN  ON  INTERRUPTS  AND  MEMORY  PROTECT 

DEF 

FIPST 

♦ 

RETURN  TO  CALLING  PROGRAM 

END 

FIPST 
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Appendix  J:  Segmented  Program 

Program  MAINl  will  load  and  execute  one  of  three 
disk-resident  segments.  The  user  is  requested  to  enter 
1,  2,  or  3  corresponding  to  SEGM1 ,  SEGM2,  or  SEGM3 , 
respectively . 
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_ 


r 


FTN4*  L 


_  ♦ 
c* 
r  ♦ 


50 
1  0  0 


PROGRAM  MR  INI 


THIS  PROGRAM  DEMONSTRATES  A  METHOD  "0=  WRITING  SEG¬ 
MENTED  PROGRAMS.  IT  USES  THE  PROGRAM  SEGMENT  LOAD 
EXEC  CALL.  ONE  OR  THREE  SEGMENTS  CAN  EE  CALLED. 


INTEGER  SE5M1X3')  •  SEGM2  0>  •  SEGM3 

DATA  SEGM 1  ■  £HS E ..  2HGM *  1  HI 
DATA  SEGMS  EHSE • 2HGM • 1  HE 
DATA  SEGM 3 -'2HSE  •  2HGM >  1  His 
DATA  I  CODE  3-" 

ASL  USER1  WHICH  SEGMENT 

WRITE  '.i ,  1  00.' 

eOc'1AT  'V»  "WHICH  SEGMENT  1*  2*  OP  3T 
READ 'li*  *  I CHOSE 


.Or 3  AND  EXECUTE  SEGMENT 


EC !.  1  • 

CALL 

ST 1  'C  p 

EC.  2.'  CALL 

E,  :EC 

EC.  l> 

CALL 

EXEC 

rrogram  IF 

EMTR 

EC  a CODE* SEGM 1 ■ 
EC  '.I CODE*  SEGM 2* 


2  0  0 


WRITE  ■  1 «  200.' 

FORMAT '  "ONLY  TH?EE  SEGMENTS  AVAILABLE. 

+  "EXECUTIVE  oROGrAM  TERMINATED. " ■ 

STOP 

END 


-TN4.L 


PROGRAM  SEGM1 <5> 


L  ♦ 

c* 

s  ♦ 


THIS  IS  SEGMENT  1  Oc  PROGRAM  MAIN1 


♦ 

♦ 


MP I TE  < 1 • 10 0 1 

100  CCR,-AT 1  v,  "SEGMENT  I  HAS  BEEN  LOADED  PROM  THE  "  • 
+ " D I S*  AND  EXECUTED" •  ■ 
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RETUPN  TO  MAIN1 


CALL  MR INI 
END 


FTN4»  L 


PROGRAM  3  EG  ME <5> 


SIP  I TE  1  1  <  1  0  CO 

100  -ORMAT  </  ■•’■>  " SEGMENT  £  HR'*  BEEN  LOADED  FROM  THE 
+“DI2K  AND  EXECUTED " «  > 

C 

I -  RET'JRM  TG  MR  I  HI 


CALL  MAIN  1 
END 


- TH4 • L 


PROGRAM  SEGM3 ■ 5> 


THIC  i:  SEGMENT  3  OF  PROGRAM  MAIM  1 . 


..if  I TE  1  •  1  0  0 '' 

100  FORMAT  <v,->  ''SEGMENT  3  HR"  SEEM  LOADED  F=GM  THE 
+  "DI'1-:  AMD  EXECUTED" «  :■ 

l“ 

C -  PET'JFN  TO  MR  ITU 

CALL  HA INI 
END 
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Apperu,  x  K:  Executive  Programs 

Two  programs  were  implemented  to  demonstrate  the 
system's  capability  for  executing  a  large  number  of 
DSP  applications  programs  by  means  of  an  executive 
program.  SCHED  demonstrates  the  use  of  the  Program 
Schedule  EXEC  Call.  The  user  has  the  option  of  specifying 
ore  of  four  types  of  scheduling.  MESG1  demonstrates 
the  use  of  the  System  Message  utility,  MESSS.  MESSS 
can  be  used  to  execute  any  of  the  system  commands.  How¬ 
ever,  only  the  RU  command  is  included  in  MESG1 . 


r 


-TN4.L 

PROGRAM  CCHED 


THE  UOEF:  SPECIFIC-  THE  TYPE 
I  QF  THPEE  PROGRAMS  TO  BE 


c  ♦ 

£:♦ 

THIS  PROGRAM 

SCHEDULE  EYE* 

c* 

0-  SCHEDUL I  Ni 

♦ 

SCHEDULED. 

2+ 

♦ 


♦ 

♦ 


DIMENSION  MAME1  OS'1  *  NAMES  C2'>  .NAMES  O'* 

DATA  NAME 1  SHF®. SHOE* 1H1 
DATA  NAMES  £HDP «  SHOE .  1  He 
DATA  NAMES  SHF®. EHOG* IMS 

■I -  SELECT  T  '»£  OF  SCHEDULE 

SC*  M=ITE  I.  IOC" 

100  PCFMAT  .  SELECT  TYPE  CF  SCHEDULE. " • 

+  "  ?-I  imEDIAtE  SCHEDULE'  HITH  HAIT"- •> 

*  “  l  0*  I MMED I  ATE  SCHEDULE'  NO  MA I T  ”  **  ^  • 

*  23**"JEUE  SCHEDULE*  MITH  UA IT"*  * 


••£4=QU 

EUE 

"  ;* 

:hedule«  no  ha  it" * 

“EAD *  1 « ♦ 

ICODE 

IF ■ I  CODE 

.ED.  ?■ 

GO 

TO 

GOO 

I-  I  CODE 

. ES . 1 0> 

GO 

TO 

3  0  0 

IF  *  MODE 

.  E  S>.  ES'* 

GO 

TO 

300 

. 

I- ■ ICODE 

.EC*.  E4) 

GO 

TO 

3  0  0 

: — 

IN'  'A*_ID 

input 

■P I TS  t . 

1  000  • 

GO  TO  50 

c  — 

SELECT  P 

PQGFAM 

300 

UP  I TE  C 1  * 

400  * 

4  0  0 

cOF'MAT  -•  ✓ 

'  •  ’"SELEC 

T  ® 

ROC 

'p'AM  TO  BE  RUN. 

+ 

"l*PRO 

Gl " 

«'■’<» 

**^sPf?2*ji***  *  **  £ — P‘ C' 3 ^ 3 

PEaD a . ♦ 

'•  IFFDG 

GO  TO  *5 

0  0  *  6  0  0 «  7 

0  0  > 

IP 

>*□5 

INVALID 

INPUT 

450 

HP  I TE  C 1 > 

1  0  0  0  * 

GO  TO  30 

0 

c - 

SCHEDULE 

PROGRAM 

500 

CALL  E  :E 

C  < I  CODE  < 

MAN 

El  • 

L 


•o 
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GO  TO  730 

600  CALL  E'  :EC  a CODE*  NAMES1 
30  TO  730 

700  CALL  EXEC < I CODE « NAME 3 > 

730  MRITEO  >  300- 

300  cOPMAT  K.s/r  "DO  YOU  WANT  TO  SCHEDULE  ANOTHER  PROGRAM? “  >  /« 
+■  “1  =Y£ C  c=NC  “ » /  > 

READ  1 "  ♦  >  IPICK 
IFCIPICK.EQ. 1>  GO  TO  50 

1000  FORMAT  (.■">  ”♦♦♦♦♦  INPUT  EPROP  ♦ 

t-‘ 

STOP 

END  ' 

eTN4,L 


PROGRAM  PR CGI 


,  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦< 

;♦ 

>  THIS  PROGRAM  IS  SCHEDULED  BY  THE  PROGRAM  "SCHED" 


■  ♦♦♦♦♦♦♦♦♦♦♦•< 


URITE  '  1 1.  1 00' 

FORMAT "PROGRAM  1  HAS  BEEN  SCHEDULED  AND  RUN". 


STOP 

END 


-TN4.  l 


PROGRAM  CRGG£ 


.w  •♦♦•*♦♦♦♦♦♦♦♦♦♦♦♦♦♦< 


♦  ♦ 

THIS  R'C'03R AM  IS  SCHEDULED  BY  THE  PROGRAM  "SCHED".  * 

*;♦  ♦ 

c •♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 


URITE-n .  1  GO  ' 

100  FORMAT  f./s*  "PRDGP.AM  2  HAS  BEEN  SCHEDULED  AMD  C1JN"  >  ■'  :• 

~ 

STOP 

END 

FTN4.L 

r- 

£ 

PROGRAM  PROGS 
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O  THIS  PROGRAM  IS  SCHEDULED  BY  THE  PROGRAM  "SCHED".  ♦ 

£♦  ♦ 
C ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦ 

i"* 

WRITE  a*  100> 

100  FORMAT  <ss? "PROGRAM  3  HAS  BEEN  SCHEDULED  HMD  RUN'S //) 


FTN4 '  L 


PROGRAM  MESG1 


c*  • 


c*  • 

C*  THIS  PROGRAM  DEMONSTRATES  THE  USE  OF  THE  SYSTEM  ♦ 

O  MESSAGE  UTILITY  <MESSS> .  THE  UTILITY  ALLOWS  ♦ 

C*  SYSTEM  COMMANDS  TO  EE  EXECUTED  FROM  A  USER  PROGRAM.  ♦ 

IN  THIS  CASE?  THE  USER  HAS  A  CHOICE  OF  RUNNING  ONE  ♦ 
C>  OF  THREE  PROGRAM.  ♦ 

C  ♦  ♦ 


♦ 


I HTEGER  CMND 1 <3> *  CMND  2 <Z> • CMNDS <3> 


£ 

: -  CMND  1  =  "RU  j  PROG  1 '' 

C -  CMND  2  *  "PIJ«  PR0G2" 

S -  CMNDS  =  “R'Jj  PROGS" 

S -  SELECT  PROSPAM  TO  BE  RUN 

50  WRITE  a j 1 Q0> 

100  POP'MAT  <ysi  "SELECT  A  COMMAND" «  "  1=RU«  PROS1 " «  * 

+*  "  S=P’U  .  RRObd  "  i*  ?  "  bsRU  !•  PRC'j  j  "  •  ■* 

READ  1 1<  »■  IFICK 

GO  TO  1 2 0 0 *  S 0 0 j 4  0 0 >  I F  I C K 

S -  RUN  PROGRAM  'PR'OGl' 

200  CMND1<1>  =  2HRIJ 

•:mnD1<2>  =  2H»P 
•SMND1  <3>  =  2HRO 
CMND  1 '4'.  =  2HG1 
I  INTI  =  3 

I  =  ME  SSS ■ CMND 1 » I CNT 1 • 

IFCI.LT.  0>  URITS<1'£000'  <CMMD1  <J>  «  J*l«  3') 

IP  a.GE.  O’)  WRITE  a.  2  000:.  .  CMNDl  <J>  * 

GO  TO  ?00 

S -  RUN  PROGRAM  '■RR0G2' 

c 

300  CMNDS C 1 >  =  cHPU 

CMNDS  '-2  •  =  2H.P 
CMNDS  cS  =  2HRO 
CMNDS  ■:  4  >  =  2HG2 
I2NT2  =  3 

I  =  MESSS  fCMND2»  ICNT2.) 

I F a  . LT .  0 )  WRITE  <’1>£00 0)  <CMND£ <J)  .  J  =  1 .  S > 

I F  ■:  I .  SE .  0.)  WR I TE  <1  •  3 0  0 0  >  0IMND2  <  J>  •  J=  1  •  4  • 

GO  TO  300 

C -  RIJN  PROGRAM  'RPOG2 
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f'j  o 


400 


CMND3 <1>  =  2HPU 
CMND3  '■£'>  =  2H  ?  P 
CMND3  C3-*  =  2HR0 
CMHD3C4-)  =  2HG3 
I3NT3  =  3 

I  =  MESSG  03MND3» ICNT3> 

IF  <1 .  LT.  0)  WRITE  a  » 2000>  CCMMD3  <J>  >  J=  1  •  3> 
IF  <1 . GE. 0>  WRITE <1» 3000  >  03MNB3 <J> ? J=1 -  4  ■ 
C 

C -  fl'SK  Ic  ANOTHER  PROGRAM  13  TO  EE  PIJN 

900  WRITE <1? 4000) 

READ 1  -  IPICK 
IFaPICK.EQ.  It  GO  TO  50 

000  FORMAT  </*',  3A2»x  •-> 

000  c0PMAT ('/> "THE  COMMAND  ' " »  4A£, “ '  HA 3  EE 
4000  FORMAT  - v?  "DO  YO'J  WANT  TO  CUN  ANOTHER  ° 

+  “  1=YES  2=NO"  < 


3  TOP 
END 
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Appendix  L:  Recommended  Software  Structure 


The  software  structure  recommended  for  the  DSP  lab¬ 
oratory  includes  the  Program  Schedule  EXEC  Call  and 
segmented  programs.  Program  EXECT  (the  executive)  has 
five  valid  commands:  TERM  (terminate  the  program), 

PRGM1 ,  PRGM2 ,  PRGM3 ,  and  LIST  (lists  the  valid  commands). 
The  programs  PRGM1 ,  PRGM2 ,  and  PRGM3  are  identical 
except  for  the  names  of  the  segments  they  can  call. 
Therefore,  only  PRGM1  and  its  segments  -  SUBRl ,  SUBR2 , 
and  SUBR3  -  are  listed. 


ctna.l 


PROGRAM  EXECT 


THIS  IS  AN  EXECUTIVE  PROGRAM  WHICH  WILL  SCHEDULE  ♦ 
DISK  RESIDENT  PROGRAMS  FOR  IMMEDIATE  EXECUTION.  ♦ 
THERE  APE  FIVE  VALID  COMMANDS:  PRGM1 ?  PPGM£,  PRGM3,  ♦ 
LIST?  AND  TERM.  THE  FIRST  FOUR  WILL  SCHEDULE  ♦ 
PROGRAMS.  "TERM-'  CAUSES  EXECT  TO  TERMINATE.  ♦ 


C OMMON  I CMND <3> » I  CODE 


1  0  0 

150 
2  0  0 


4  0  0 


*♦♦♦♦♦  c'RObF'AM  •“  OF 
>♦"••11  </>  1 


DIGITAL  SIGNAL 


I CODE  =  9 
WRITE  a  •  100) 

FORMAT  <  1  0  <•-")  '  1  OX, 

+  "PROCESSING  ♦♦•* 

WRITE  a*  £00) 

FORMA T  <y " EN TEC'  COMMAND  " 

I CMND <1)  =  0 
I  CMND =  0 
I  CMND '■GO  =  0 

READ  <1,4  0  0  :•  '■!  I  CMND  <  I  -  *  1  =  1,  So 
FORMAT  <3A£) 

I  r  <  < I  CMND  <  1  >  .  S'? .  2HTE>  .  AND .  <  I  CMND  <  2.» .  EC! .  E'HRM)  )  STOP 
I F  ( I  CMND  <  1)  .SO.  2HP*>  .  AND .  <  I  CMND  •'  2  1  .  EC* .  £HGM> 

+  . AND. < I CMND <3> . EC. 1H1 ■ >  GO  TO  500 

IP<  <  I  CMND  <.1>  .  EG.  2HPR)  .  AND.  <  I  CMND  •  2)  .  EC.  2HGM) 

+  . AND. <1 CMND < 3> .EG. 1H2>>  GO  TO  500 

I F  <  C I  CMND  ’■  1  >  .E0.2HFP)  .  AND.  Cl  CMND  '‘S'  .  EG.  2HGM> 

+  . AND.  Cl CMND  <3> . EG. 1H3) >  GO  TO  500 

I F  (  <  I  CMND  c  1  > .  EG .  *HL  I  .  AND .  <  I C  MND  .  EG .  2HS  T ) '•>  GO  TO  5  0  0 


S -  INTERCEPT  INVALID  COMMANDS 

I* 

CALL  ERR  0=1 
GO  TO  150 

C -  SCHEDULE  c'c'OGRAM 


500  CALL  EXEC  <  I  CODE  •  I  CMND'1 
GO  TO  150 
END 


SUBROUTINE  ERROR 
COMMON  I  CMND  OS"'  *  ICODE 

i* 

WRITE  a  *  1 00 > 

100  -OFMAT-r  ”*♦♦♦♦  IM='JT  E  =  R0c 
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£00 


300 

400 


MR I TE <1 -200' 

FORMAT <' '"SELECT  3PTI0MS  ”  “  1  =  STOP"  "2  =  C3HTINUE" 
+  "3  *  DISPLAY  LIST  OF  COMMANDS" 

READ  <1 - ♦>  IPICK 
I F  a  P I CK . EQ . £>  RETURN 
IF  < IPICK. EQ. 3>  30  TO  300 
STOP 

WRITEa -400> 

FORMAT  "♦♦♦♦♦  LIST  OF  COMMANDS  ♦♦♦♦♦” 


+  " PRGM 1 " - 5X > “ PRGM2 " « 5X > " PPGM3 " - 5X - " L I ST " • £X . “ TERM " 


RETURN 


END 


FTN4>  L 


PROGRAM  PRGM 1 

COMMON  -- 1  FLAG  -" 1  FLAG  1 

INTEGER  SUER  l  <3>  <  SUERS'  <3>  *  SUERS  C3> 

DATA  SUER  1  ■'2H3U «  £HER  ?  1H1 
DATA  S U E R £  •■' £ H S U *  •£ H E P *  1H£-' 

DATA  SIJER3-'£H  I U *  2HBR »  1 H3-" 

DATA  I  CODE  •  I  EG  I  \\/ 0 

IFCIBGIN.NE. 100)  GO  TO  100 
GO  TO  ':50j51»5£>  I  FLAG  1 
100  I EG IN  =  100 

WRITE  a  ,20  0> 

£00  FORMAT'. -••PRGM  1  NAS  BEEN  SCHEDULED  BY  ELECT"' 

S -  LOAD  AND  EXECUTE  FIRST  SEGMENT 

CALL  E:  .EC  I  CODE > SUER  1  > 

C 

C -  LOAD  AND  EXECUTE  SECOND  SEGMENT 

CALL  E"  'EC  ICODE-  SUBR£> 

LOAD  AND  EXECUTE  THIPD  SEGMENT 

CALL  EXEC  <  I  CODE  -  SUERS') 

STOP 
END 

BLOCK  DATA 
COMMON.-'  I  FLAG/- 1  FLAG  1 
END 

FTN4.L 


PROGRAM  SUER  1  0-'-  N*MED  COMMON 
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COMMON-' I  FLAG-'  I  FLAG  1 


write  a.  100> 

100  FORMAT  C'"SUBPl  HAS  BEEN  LOADED  AND  EXECUTED  * 

I™* 

C -  SET  I FLAG 1  FOR  ENTRY  POINT  IN  PRGM1 

C 

I FLAG  1  =  1 
CALL  PRGM1 
END 

-TN4<  L 


PROGRAM  SIJBRE  <5>  »  NAMED  COMMON 

C 

COMMON-"  I  FLAG  "  I  FLAG  1 
WRITE  Cl j 1 00> 

100  FORMAT  ”SUBR£  HAS  BEEN  LOADED  AND  EXECUTED' 

i_! 

: -  SET  I FLAG l  FOR  ENTRY  POINT  IN  PRGM1 

I FLAG 1  =  3 
CALL  PRGM1 
END 

eTN4.L 

C 

r 

PROGRAM  SUERS  <5> •  NAMED  COMMON 
COMMON  ■"  IFLAG--'  I CLAG 1 
WRITE1  1  >  !  0 0 '■ 

100  ~Oc'MAT  <  S'JEC 3  HAS  BEEN  LOADED  AND  EXECUTED' 

C 

C -  SET  ICLAG1  FOP  ENTRY  POINT  I,N  PRGM1 

I FLAG 1  =  3 
CALL  PRGM1 
END 
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